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^ ; Abstract 

"^ ! We show that matter yukawa couphngs of the New Minimal Supersymmetric (SO(IO)) 

I GUT(NMSGUT) are subject to very significant GUT scale threshold corrections. Including 

' p-i 'l these threshold effects relaxes the constraint y^ — y,- ~ y^ — y^ operative in 10 — 120 plet 

P^' generated tree level MSSM matter fermion yukawas yf. We find accurate fits of the MSSM 

p ■ fermion mass-mixing data in terms of NMSGUT superpotential couplings and 5 indepen- 

^ . dent soft Susy breaking parameters Mq, M1/2, Aq, Mfj ^ at Mx- The fits generally have 

I — \'. elevated unification scale Mx near Mpianck, viable values of a^i^Mz), and are consistent 

■ with current limits on B violation, b — )■ 37, muon magnetic moment anomaly and Standard 

>- ' Model p parameter. The associated novel and distinctive soft Susy spectra have light gaug- 

^ ■ inos, a normal s-hierarchy and Bino LSP. The Bino LSP is accompanied by second and first 

Q\ . generation right chiral sfermions light enough to mediate a consistent WIMP dark matter 

^ co-annihilation cosmology and to be discoverable at LHC, while third generation sfermions 

r^ [ are in the LHC undiscoverable range of 3-50 TeV. The fits found require |/i|, \Aq\ ~ 100 TeV 

which imply both deep CCB/UFB minimae and stability of the MSSM standard vacuum 

on cosmological time scales. Our results indicate that a consistent realistic phenomenology 

may be specifiable in terms of SO(IO) (NMS)GUT parameters at Mx alone and that a 

k> i new viable sector of the soft supersymmetry parameter space may exist if flavour violation 

H ' constraints can be satisfied in the 43 dimensional parameter space. 

_ _ I 

1 Introduction 



Supersymmetric Grand Unification based on the SO(IO) gauge group has received well 
deserved attention over the last 3 decades. Models proposed fall into two broad classes 
: those that preserve R-parity down to low energies [H [21 [6l O [121 13 [Z] using as Higgs 
the special representations (126) of SO(IO) that contain R-parity even SM singlets and 
another large class of SO(IO) R-parity violating models[8] that attempt to construct viable 
models using sets of small SO(IO) representations even after sacrificing the vital distinction 
provided by R-parity between matter and Higgs multiplets. Besides structural attractions, 
such as the automatic inclusion of the conjugate neutrino fields necessary for neutrino 
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mass, SO(IO) GUTs offer a number of other natural features. Among these are third 
generation yukawa unification f3] [3], automatic embedding of minimal supersymmetric Left- 
Right models, natural R-parity preservation^ down to the weak scale and consequently 
natural LSP WIMP dark matter, economic and explicitly soluble symmetry breaking at 
the GUT scale [0], explicitly calculable superheavy spectra [HI [IHl [HI [12], interesting gauge 
unification threshold effects [121 [ISl [3 [H HZ] which can lead to a natural elevation of the 
unification scale to near the Planck scale [H], GUT scale threshold corrections to the QCD 
coupling a^i^Mz) of the required [18 j sign and size [19] and a deep interplay between the 
scales of Baryon and Lepton number violation as suggested by the neutrino oscillation 
measurements and the seesaw formulae connecting neutrino masses to the B-L breaking 
scale. 

The fascination of the MSSM RG flow at large tan/3 stems from the tendency of third 
generation yukawa couplings to converge, at the MSSM unification scale [3l[l], in a manner 
reminiscent of gauge unification in the MSSM RG fiow[20| |2T] . For suitably large tan (3 
and for close to central input values of SM fermion couplings at the Susy breaking scale 
Ms ~ Mz, third generation yukawas actually almost coincide at Mx- On the other hand, 
in SO(IO) theories with only the simplest possible fermion mass giving Higgs content (a 
single 10-plet), when all the complications of threshold effects at Mx ~ 10^^G'e\^(not to 
speak of those at seesaw scales Mp ~ 10^ — 10^^ GeV) are ignored, one does expect to 
generate boundary conditions for the gauge and yukawa couplings that are unified gauge 
group wise and (third generation) flavour wise. 

However, fitting the rest of the known fermion data (15 more parameters) definitely 
requires other Fermion Higgs multiplets (more 10-plets, 120,126 s etc). A principled 
position {monoHiggisml) with regard to the choice of Higgs multiplets responsible for 
fermion mass (FM Higgs) is to accept only one of each irrep present in the conjugate of the 
direct product of fermion representations. This principle may be motivated by regarding the 
different Higgs representations as characteristic "FM channels" through which the fermion 
mass (FM) is transmitted in structurally distinguishable ways. For example the Georgi- 
Jarlskog mechanism distinguishes the 45 plet Higgs in SU(5) (126 in SO(IO)) from the 
5 + 5 (10 in SO(IO)) due to their ability to explain the quark-lepton mass relations in 
the second and third generations respectively. Similarly the 126 in SO(IO) is peculiarly 
suitable for implementing the Type I and Type II seesaw mechanisms for neutrino mass. If 
one duplicates the Higgs multiplets transforming as the same gauge group representation, 
for example by taking multiple 10-plets in SO (10), then one abandons the quest for a 
structural explanation of the pattern of fermion masses in favour of "just so" solutions. 
Moreover the many Higgs doublets typically present in the theory mix with the 10-plet 
derived Higgs doublets to make up the light doublets of the effective MSSM. Both these 
features dilute the expectation of exact third generation yukawa unification. 

The technical difficulties of 6 — r unification are attendant upon the fact that the GUT 
scale values of the third generation Yukawa couplings, in sharp contrast to the three gauge 
couplings, are highly infrared sensitive(see [22] for a cogent discussion). The accuracy of 
unification of third generation Yukawa couplings after the MSSM RG flow up to the 1-loop 
GUT scale M^(~ iQ^^-^s GeV) is strongly dependent on the low energy values of the yukawa 
couplings in the MSSM and thus requires large tan/3. It is also strongly dependent on the 
precise values of y^^^^{Mz) and yf^ ^^'^^Mz) {hoih of which can receive significant Susy 
threshold corrections) besides a^{Mz) (which needs significant GUT threshold corrections 



to be consistent with experiment [18] ) . The large threshold corrections [U [23] to the relations 
between Standard Model(SM) and MSSM Yukawa couplings for down type quarks precisely 
in the large tan/3 MSSM variants that are in focus in a SO(IO) context, combined with 
infrared sensitivity of yukawa unification casts doubt on the naturalness of exact t — b — r 
unification. As these yukawa couplings have been measured more precisely over the same 
period as the investigation of third generation yukawa unification in Susy SO(IO) models, 
the strain of these models has increased with time, and various just so mechanisms are 
invoked to preserve the exact third generation unification. Yet it would be strange if 
Susy SO(IO) had no mechanisms to accommodate not only this difficulty but even the full 
complexity of the fermion spectrum. 

Many detailed investigations [H H] of the yukawa unification problem since the early 
days have focussed on models based on small Higgs representations(16, 10, 45) and seek 
to understand only the third generation fermion masses. These studies assume that no 
Higgs Multiplet other than the 10-plet could contribute significantly to the third genera- 
tion masses and use that to justify the assumption that high scale threshold corrections 
to the yukawa couplings are negligible. In the few studies that attempt to tackle the full 
fermion spectrum, recourse is had to a subset of the many possible non-renormalizable 
operators that contribute to fermion yukawa couplings in the effective theory below the 
unification scale. Since the first two generation yukawas are so small relative to the third 
generation it is possible to build models using small non-renormalizable contributions that 
allow one to fit three generation charged fermion data[2H 122] and even keep the operator 
dimension five contributions to baryon violation within or near to experimental limits [25]. 
When other aspects, such as GUT scale spontaneous symmetry breaking, are to be ad- 
dressed multiple Higgs 10,45,16-plets are often introduced [8]. These then require ad-hoc 
discrete symmetries to replace the R/matter parity that is broken in such theories. These 
symmetries also serve to prevent interference between the assigned functions of the multiple 
copies of low dimensional Higgs representation introduced (along with appropriately lim- 
ited non-renormalizable interactions that in effect play the role of the higher dimensional 
Higgs irreps that were disallowed at the outset). With stringent unification constraints at 
GUT scales and large threshold corrections at the Susy breaking scale only narrow regions 
of the Susy-GUT parameter space remain viable. Moreover specification of the soft Susy 
spectrum and couplings to ensure simultaneous cancellations to keep the threshold correc- 
tions to m^^^^^\Ms) small and the b ^ s'j rate consistent with current limits [26] l27t [22] 
is needed. In general such investigations have focussed on evaluating various non minimal 
effects such as non-universal gaugino masses, D-term contributions to scalar non universal- 
ity etc, operating in specific (tight) corners of the soft susy parameter space, to evade the 
confiicting demands of infra red sensitivity, large infrared threshold corrections and accu- 
rate yukawa unification. A much smaller set of papers [H [21 [5l [9l [6l [131 [Zl [HI [H] focusses on 
renormalizable models without ad-hoc symmetries and confronts the full complexity of the 
spontaeous symmetry breaking problem conjointly with fitting the fermion spectra [T^[T7]. 

Renormalizable SO(IO) Susy GUTs[Il [21 [S] employ not only the 10-plet but also the 
other two SO(10)-allowed fermion mass (FM) Higgs representations i.e the 126 and 120. 
Renormalizability of the GUT model is strictly maintained throughout as a criterion for the 
structure of the superpotential. As a result these theories have a very parameter-economical 
structure [H [21 [SI [H]. In 1992, just as yukawa unification came into focus it was proposed [29] 
that use of only the 10 — 126 FM Higgs might be sufficient to account for the entire fermion 



mass spectrum including the neutrino masses and mixings in terms of their (15) yukawa 
couphng parameters. This obviously attractive possibihty was investigated as a generic 
possibihty in SO (10) (i.e without using the formulae in terms of the fundamental GUT 
parameters but only having the symmetries and (naive [50] ) parameter counting implied by 
their FM Higgs yukawa structure). Although initially unsuccessful, increasingly refined 
studies based upon the increasingly well defined neutrino mass data ultimately showed [ 



that the generic 10 — 126 parametrizations permitted accurate fits of all the known 
fermion mass data using 10-plet domination of third generation mass terms combined 
with Georgi-Jarlskog mechanism for the second generation and either the Type I or Type 
II seesaw mechanisms or a combination thereof to fit the neutrino data. Perforce these 
generic models always assumed freedom to dial the relative strength of the two seesaw 
mechanisms. These accurate fits were found even while assuming that the precisely known 
SM yukawa coupling uncertainties persevered unchanged over the RG flow up to the GUT 
scales. In the light of the drastic susy threshold corrections at large tan P this assumption 
was naive and peculiarly uninformed of the attention long paidjH [23] to the role of these 
corrections in the other strand of SO(IO) GUTs. The appearance of [331 ^^ however 
reminded us of this essential feature and we have thenceforth adopted the much more 
realistic estimates of theoretical and experimental uncertainties advocated by [33] in our 
subsequent workjT?]. 

Just as the intriguing results on the successful implementation of the Babu-Mohapatra 
proposal in generic SO(IO) GUTs emerged it was also shown [T31 171 [TBI 150] that : (i) The use 
of MSGUT specific fermion mass formulae implied that neither seesaw mechanism could 
account for neutrino masses while also fitting the charged fermion data, (ii) Conclusions 
derived using generic formulae could not be transferred to concrete and specific GUTs 
because of intrinsic difficulties in untangling highly nonlinear constraints placed by the 
underlying fundamental GUT on the generic parameters, (iii) In a specific scenario [7] 
inclusion of the remaining allowed FM Higgs (i.e the 120-plet) could permit the generation 
of large enough Type I seesaw masses due to suppressed 126 couplings and charged fermion 
mass fits due to the 10 — 120 combination. 

The 10 — 120 FM Higgs combination proposed by us to tackle the charged fermion 
fits (with the 126 couplings too small to affect any but the first generation masses) was 
seemingly shown [35] to face a generic difficulty in providing large enough strange and down 
quark masses. This was only to be expected in any scenario where the Georgi-Jarlskog or 
other mechanism to distinguish second generation down type yukawas is not implemented. 
This difficulty appears in a quite different light when one considers [T7] that Susy threshold 
corrections at large tan /3 tend [U [23] to drastically modify the effective value of the down 
and strange quark yukawas. Thus we proposed[17j that the excellent fits obtainable for the 
other (16) fermion data (besides yd^s) justified searching for the susy threshold corrections 
to implement the required corrections to lower yd^s and this could be achieved not only 
with free soft susy parameters at Ms (as in version 1 of this paper [17]) but even with just 5 
GUT compatible A^ = 1 supergravity(SUGRY), Non Universal Higgs Masses(NUHM), type 
soft parameters (mi , mo, Aq, rri^ g) at mz specified at My[T7] and two more (|/i|, -B ~ m\) 
determined by electro-weak spontaneous symmetry breaking and run back up to Mx to 
provide 7 parameters at Mx coding all the SUGRY-NUHM information. 

These developments have put detailed consideration of the role of 120 — 126 FM Higgs 
representations firmly on the agenda of Susy SO(IO) yukawa unification. In this paper 



we argue that the assumption that third generation fermion yukawas are protected from 
large GUT scale threshold corrections associated with non-10 -plet FM Higgs irreps is 
facile. We show that it is behed in the class of Susy SO(IO) theories|Il El El El El [H] 
which actually face up to the task of accounting for all the available fermion mass data 
in a fully specified model without invoking uncontrollable higher dimensional operators or 
ad-hoc symmetries and rely solely on SO (10) gauge symmetry, supersymmetry, and struc- 
tural (parameter counting) minimality as guiding principles. In particular the NMSGUT 
with a large 120-plet coupling (the 120-plet yukawa is antisymmetric and hence has two 
eigenvalues that are equal in magnitude and one that is zero) requires evaluation of the 
GUT scale threshold corrections to the fermion Yukawa couplings. The large number of 
fields, 120-plet couplings comparable to the 10-plet couplings, and the fact that these 
threshold corrections arise from chiral supermultiplet wave function normalization(so that 
any field that couples to a matter field or a Higgs doublet can run in the self energy loop) 
raises the possibility that these corrections may be far from negligible |36j. In this paper 
we show that this is in fact the case. 

The particular importance of the wave function corrections for the fermion data fitting 
program is that they can relax the stringent constraint yi, — Ut — Us — Vfi that we found 
mills] operative at Mx in SO(IO) models with a 10-120 FM Higgs system. Thus at least 
any model that employs the 120-plet must be sufficiently specified to allow the calculation 
of these very significant corrections which can drastically change the yukawa couplings at 
Mx- Such a specification involves solution of the spontaneous symmetry breaking, calcu- 
lation of the GUT spectrum and couplings and a RG analysis of threshold effects based 
thereon. So far, to our knowledge, such a complete calculation is possible only for the 
SO(IO) NMSGUT [7l[IllE]for which we give the dominant 1-loop contributions ( certain 
sub-dominant mixing effects are postponed to a future calculation). The calculation of the 
(wave function) renormalization that gives rise to the threshold effects modifying the mat- 
ter fermion Yukawa couplings highlights the issue of perturbative consistency. One finds 
that while the corrections to the matter fermion lines of the mass generating Higgs- Yukawa 
vertices are no more than 25%, the wave function corrections(Aj|^ jj) to the (light) MSSM 
Higgs lines can be much larger. The most serendipitous scenario would be if searches that 
restricted these corrections to be less than 1 were able to find fits that are also consistent 
with Baryon decay limits. This is not what we have observed in our searches so far : accu- 
rate fits with I Aji^ jjI < 1 are obstructed in achieving small values (<< 10~^^GeV~^) of the 
Baryon violating dimension 5 operator Wilson coefficients and hence give baryon lifetimes 
~ 10^^ years. On the other hand if one allows large values of Ajj ^ the searches achieve life- 
times that are a million times or more larger. Thus although threshold corrections to gauge 
couplings turned out to be mild due to cancellations [121 H] i^i spite of early alarms [36] the 
wave function renormalization effects at one loop are indeed large. Investigation of higher 
loop contributions is thus be called for to see if the perturbation theory is sensible. In 
this paper we take the view that it is necessary to pursue the investigation of the realistic 
features of the SO(IO) while including the large 1-loop threshold effects even though full 
perturbative convergence may take long to prove ( or indeed may never be possible : as for 
instance in the most precise known theory QED). 

In Section 2 we briefly review the structure of the NMSGUT[7l HH |l7j to estabhsh 
the notation for presentation of our results on threshold effects in Section 3 and Appendix 
A. In Section 4 we present illustrative examples to underline the signiflcance of the GUT 



scale threshold effects and the need to include them. In Section 5 we present examples of 
improved fits(specially with respect to acceptable d = 5 operator Baryon violation rates). 
In Section 6 we discuss the broad features of the emerging phenomenology of the NMSGUT. 
In Section 7 we summarize our conclusions and discuss the various improvements in the 
fitting, RG flows and searches and that are called for. Appendix A contains details of the 
calculation of threshold effects at Mx- Appendix B contains a discussion of the two loop 
RGE flow at large Aq, {M"^ gl values which results in the unconventional susy spectra with 
normal s-hierarchy and gaugino masses not even close to the 1:2:7 ratio found at one loop 
in Susy GUTs with universal gaugino masses. Appendix C contains tables of parameter 
values for additional example flts to provide a wider view of the possibilities and as inputs 
for phenomenological explorations of the viability of our results. . 

2 NMSGUT recapitulated 

The NMSGUT ^14j is a renormalizable globally supersymmetric 5*0(10) GUT whose Higgs 
chiral supermultiplets consist of AM(Adjoint Multiplet) type totally antisymmetric tensors: 
210{^ijki), 126(SijHm), 126(Sjjfc/„)(z,j = 1...10) which break the 50(10) symmetry to 
the MSSM, together with Fermion mass (FM) Higgs 10 (H^) and 120(0^^^). The T26 
plays a dual or AM-FM role since it also enables the generation of realistic charged fermion 
and neutrino masses and mixings (via the Type I and/or Type II Seesaw mechanisms); 
three 16-plets ^a(^ = 1,2,3) contain the matter including the three conjugate neutrinos 
(z/f ). The superpotential (see[6l [9l [TOl [121 13 [HI [U] for comprehensive details ) contains 
the mass parameters 

m : 210^ ; M : 126 ■ T26; Mh : 10^; mo : 120^ (1) 

and trilinear couplings corresponding to the superfleld chiral invariants indicated : 



A 

k 



l3 



210-^ ; r] : 210 ■ 126 ■ 126; p : 120 ■ 120 ■ 210 



10 120 -210 ;7 ©7: 10 -210 (126 ©126) 

120-210- (126©T26) (2) 



In addition one has two symmetric matrices h^s, Iab of Yukawa couplings of the the 
10,126 Higgs multiplets to the 16a-16b matter bilinears and one antisymmetric matrix 
qab for the coupling of the 120 to 16a-16b- One of the complex symmetric matrices can 
be made real and diagonal by a choice of SO(IO) flavour basis. Thus initially complex FM 
Yukawas contain 3 real and 9 complex parameters. Five overall phases (one for each Higgs), 
say those of m, M, A, 7, 7, can be set by flxing phase conventions. One(complex parameter) 
out of the rest of the superpotential parameters i.e m, Mh, M, nio, A, r], p, k, 7, 7, C, C 5 say 
Mh, can be flxed by the flne tuning condition to keep two doublets light so that the effective 
theory is the MSSM. After removing un- physical phases this leaves 23 magnitudes and 15 
phases as parameters : still in the lead out of any theories aspiring to do as much[B]. As 
explained in[6l HI [12] the flne tuning flxes the Higgs fractions i.e the composition of the 
massless electroweak doublets in terms of the (6 pairs of suitable) doublet flelds in the GUT. 
A subtle point here is that even if the other parameters are taken real the flne tuned Mh 
(which does not itself enter into the low energy lagrangian) will be complex. Thus strictly 
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speaking one cannot justify the use of only real superpotential parameters by invoking 
'spontaneity' of CP violation and we will not do so. 

The GUT scale vevs and therefore the mass spectrum are all expressible [01 [121 [ID] in 
terms of a single complex parameter x which is a solution of the cubic equation 



x-^ - 15x^ + 14x - 3 + ^(1 - x)^ = (3) 



where £ = ^^. 

^ rim 

In our programs we find it convenient to scan over the "three for a buck" [15] parameter 
X and determine ^ therefrom . Then the phase of A is adjusted to be that implied by x 
and the relation f = ^^ and is not itself scanned over independently. It is a convenient 

^ rjm ^ •' 

fact that the 492 fields in the Higgs sector fall into precisely 26 different types of SM 
gauge representations which can hence be naturally labelled by the 26 letters of the English 
alphabet |12j. The decomposition of SO(IO) in terms of the labels its "Pati-Salam" maximal 
subgroup S'?7(4) x SU{2)ji x SU{2)l provided^] a translation manual from SO(IO) to 
unitary group labels. The complete GUT scale spectrum and couplings of this theory have 
been given in [121 E] • 

The tree level fermion yukawa couplings and neutrino masses of the effective MSSM 
arising from this GUT below the GUT scale after fine tuning to keep one pair of Higgs 
multiplets light are given in [TH [T7] . 

As mentioned, in the NMSGUT the conjugate(i.e "right handed") neutrino Majorana 
masses are 4 or more orders of magnitude smaller than the GUT scale due to very small 126 
couplings. Therefore for purposes of calculating the threshold corrections to the Yukawa 
couplings at Mx we can consistently treat the conjugate neutrinos as light particles on the 
same footing as the other 15 fermions of each SM family. These fermion mass formulae, 
after correcting for threshold effects, are to be confronted with the RG-extrapolated data 
(from Q = Mz to Q = Mx = 10^^'^^ GeV including neutrino masses and mixing angles). 
The calculation of Ax also fixes the scale m of the high scale symmetry breaking [TJ [161 [H]. 
The stringent simultaneous requirements of of a common unification-seesaw scale, gauge 
unification (including the right high scale gauge RG threshold corrections to shift the GUT 
prediction of a^i^Mz) down to acceptable values[19j), third generation yukawa unifcation, 
as well as fits to all the other fermion masses and mixing matrices, are effective in singling 
out characteristic and suggestive GUT parameters (including Susy breaking parameters at 
Mx). 

3 GUT scale Yukawa threshold corrections 



The technique of[37] for calculating high scale threshold corrections to yukawa couplings, 
generalizes the Weinberg-Hall[38j method for calculating threshold corrections to gauge 
couplings, and has long been available but has not been exploited much; possibly due to 
the assumption that such effects are always negligible. In supersymmetric theories the 
superpotential parameters are renormalized only due to wave function correction and this 
is easy -if tedious- to calculate for the large number of MSSM submultiplets in the 120-plet 
which couple to the light fermions and MSSM Higgs at an SO(IO) yukawa vertex. There 
is also a contribution from heavy gauge field couplings to the light (i.e MSSM) fields at 
the matter yukawa vertices. The calculation involves going to a basis in which the heavy 



field supermultiplet mass matrices are diagonal. This basis is easily computable given the 
complete set of mass matrices and trilinear coupling decompositions given in [9l [121 El- 
For a generic heavy field type $ the mass terms in the superpotential diagonalize as : 

$ = f/*F ; $ = \/*$' ^ $^A/<I> = W^MDiag^' (4) 

Threshold correction to a Yukawa coupling matrix (which occurs in the superpotential 
a.s W = f^YffH±) then have the form 

Yf = Yf + AJ.Yf + ArYj + AH±Yf (5) 

where the ± refers to the Y = ±1 Higgs multiplets appropriate to give mass to Ts = ±| 
fermions and Ajj, Afj± are the 1-loop wave function correction factors. For a generic 
interaction superpotential W = ^ J2ijk ^jifc^*^"''^^, the quantities A, at the renormalization 
scale Q, have the form (A = —■^K in the notation of [37] ) 



AKQ) = ,^^{-2gl,Y.Fii^A,m,,Q)lfj^^ + -J2Y,,iY;,^F^{m,,mi,Q)) (6) 

"^^^ k,A ^ kl 

We have used Q = Mx = 10~^^Mx where Ax is the shift (away from the one loop 
MSSM unification value LoqiqM^ = 16.25) due to loop and threshold effects. Here A is 
a generic gauge field (adjoint) index and i,j,k are generic chiral field indices. The gauge 
couplings and generators gio,!"^ of SO(IO) are related to the usual (5'[/(5)normalization) 
gauge coupling and generators g^ = g, T^ by gio = g/\/2, I^ = \/2T'^. When both the 
fields running in the loop are heavy fields {Fi is a symmetric Passarino-Veltman function) 
Fi should be taken to be 

1 A/f^ M^ 

FMM.. Ms, Q) = ^^^^j_^(Ml/„^ - Mlln^) - 1 (7) 



which reduces to just 



/bf2 

Fn(M^, Q) = Fu{Ma, 0, Q) = In-^ - 1 (8) 



when one field is light {Mb -^ 0). When one of the heavy fields in the loop has MSSM 
doublet type G321 quantum numbers [1, 2, ±1] (so that one eigenvalue is light while the 
other /ive[Tl] are heavy) care should be taken to avoid summing over light-light loops since 
that calculation belongs to the MSSM radiative corrections. 

In the NMSGUT fitting scenario where |/ab| < 10~^ it is an excellent approximation 
to ignore the contributions of the 126 to the high scale threshold corrections. Moreover 
while calculating the wave function renormalization of the Higgs line in the matter fermion- 
antifermion-MSSM Higgs vertex, we shall assume that it is an good approximation to take 
the MSSM higgs to be dominantly made up of the 10-plet (which dominance we know is 
required in order to account for top-bottom-r (near) unification which is an intrinsic part 
of the large tan/3 scenario) so that we can ignore the admixture of the other 5 MSSM type 
Higgs doublets pairs present in the theory. Hence, for the non-dominant i.e non 10-plet 
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derived light Higgs components, the contributions to the wavefunction renormahzation of 
the hght Higgs doublets in the theory will all be suppressed by additional small (in keep- 
ing with the 10-plet dominance of light Higgs composition) factors of |a^P;^ = 2. ...6 (for 
Ho[l,2, 1]) or |q;jP,z = 2. ...6 (for Ho[l,2, —1]). That is to say there will be a suppression 
of the contribution by |ajp or |ajp unless the external doublets come from the 10-plet in 
the fclO.120.210 or 10.210.(7126 -|- 7I26) terms of the superpotential. We can therefore 
sidestep-for the moment- the elaborate calculation required for calculating the wave func- 
tion renormahzation of the external light Higgs when all its 6 possible components (from 
the 10, 210, 126, 126, 120 (two pairs) ) are corrected by wave function renormahzation. 
This is an good ( but not perfect ) approximation in practice as will be seen from the values 
of the Higgs fractions [H [121 [II]- For example in case I-l(Tabl^the values are : 

|a|2 = {0.689,0.0044,0.0038,0.1838,0.0302,0.0886} 

lap = {0.792,0.0068,0.0037,0.1020,0.0052,0.0903} (9) 



where the numbering of the components is[Tl] seen from : 

[1, 2, -l](/ii, h2, h, hi, h^, he) © [1, 2, l](/ii, /12, hs, /14, h^, h^) = (10) 

l-"2'^2 '^2 '"T/f' 2''^2 jy^l-"al)^„i i^al ' T7|"' *^«1' *^ai ) 

The decomposition of SO(IO) invariant terms in the superpotential and gauge terms 
yieldsjHl [121 E] ^ large number(~ 100 ) of vertices. It then requires a tedious but straight- 
forward calculation to determine the threshold corrections explicitly. The explicit expres- 
sions are given in Appendix A. 

Heretofore such threshold corrections have mostly been argued to be negligible(< 1%) 
although at least one paper [[27j faced with the difficulties of literal third generation yukawa 
unification has considered the possibility, without any explicit model which permitted calcu- 
lation, that the third generation yukawa unification relations must necessarily be subject to 
threshold corrections of up to 50%. In which case it was found that the various stratagems 
invoked to permit precise 3 generation Yukawa unification could become redundant. We 
shall see that the calculation of the GUT scale 1-loop Yukawa threshold effects in the 
NMSGUT can actually change the naive(i.e pure 10-plet) unification relations yt = Vb = Vt 
significantly. Furthermore the 10 — 120 plet fermion fits have been shown ( in the absence 
of GUT scale threshold effects) to require a close equality \yh — yr/iys — y^l)\ ~ 1 at Mx 
which is very constricting when searching for fits. The fits we exhibited in [TH [17] were 
all of this type. However in the present case the fits we obtain can deviate significantly 
from y''~y^ ~ 1 which was obeyed by both the fits presented in [TTl where no threshold 
corrections were applied to the yukawas at Mx- The large values of ^h,h can radically 
reduce the magnitudes of SO(IO) yukawas required to reproduce the MSSM couplings at 
Mxi thus threshold effects can help in loosening this constriction of fitting freedom. More- 
over the changes are such as to help in finding fits with a slower B violation rate. We note 
again that one must study the higher loop threshold corrections and the steps necessary 
to define a consistent perturbative expansion (possibly involving some variants of large N 
resummation and use of the exact SO(IO) Susy gauge beta functions[77|) to see if the 1-loop 
results we find are stable. 



4 Numerical fits to the fermion data and threshold 
effects 

To appreciate the significance of tfie tliresliold corrections at Mx for tfie matter fermion 
yukawas it is sufficient to consider tlie values of tlie matrices Ajj for the various MSSM 
fields when the formulae derived in the appendix are evaluated using parameters from the 
examples of fits (found ignoring GUT scale threshold corrections) given in P^ [T7] . The 
NMSGUT is , to our knowledge, the only SO (10) model where the spontaneous symmetry 
breaking and the consequent spectra and Higgs fractions have been explicitly calculated so 
that the Yukawa threshold corrections can be evaluated for specific GUT parameter based 
fits. 

The example fits from [17] are of the 18 known fermion data namely the yukawa cou- 
plings yt,b,T,c,s,^l,u,d,e, the CKM angles and phase 6*^2, 13, 23? '^'^^ the neutrino mass squared 
differences Am2i 32 and Leptonic mixing angles 6^2 13 23 in terms of the NMSGUT hard and 
soft parameters. We search assuming normal neutrino hierarchy, a very light (< 1 meV) 
electron neutrino, and a small neutrino mixing angle 9^^ < 5x1. The fits are found by a 
random search based on the downhill simplex method which requires the definition of a 
X^ function formed from the difference of the GUT implied and target (i.e RG extrapo- 
lated from Mz) values of the fermion parameters normalized by the uncertainties in these 
parameters [33] . An important point is that heretofore fitting in this class of renormaliz- 
able SO(IO) models assumed that the uncertainties involved were merely the (very small) 
error estimates for the SM extrapolated to the GUT scale. The complexities induced by 
the uncertainties due to the strong threshold effects at the Susy and GUT thresholds were 
never given any shrift in this context. Recently the papers [33] IM] motivated us to use the 
threshold effects to evade the difficulty of the failure of the Georgi-Jarlskog mechanism in 
the NMSGUT by lowering the couplings ya^s via threshold corrections due to Susy partners. 
Thus in [T7] we used the more realistic values of[3S] for the error estimates and eschewed 
the spurious precision of most previous efforts in the context of renormalizable MSGUTs. 

The first example of an accurate fit presented in [H] was able to achieve a xx = 0.0538 
for a fit of the 18 fermion data at the scale Mx (accompanied by xz = 0.027 fit at the scale 
Mz for matching the run down charged fermion yukawas to the Standard Model results 
after inclusion of the Susy threshold corrections). When the same data are used with the 
threshold corrections switched on one gets xx = 773.6 and the unification parameters 
are Ax,g,3 are unacceptable. Clearly the threshold corrections can make a great deal of 
difference ! The corresponding changes for the second solution are xx = 286364.7 and 
the unification parameters are Ax,g,3 are unacceptable. It is also worth noting that even 
though the fermion yukawas generated from the NMSGUT formulae change radically when 
one inserts the threshold corrections (using couplings determined by the fits found ignoring 
them), one finds that the accuracy of satisfaction of yt — yr = ys — Vn as measured by 
the ratio Rbr/sfi = {Vb — yT)/ys — Vii remains good with the ratios for these two solutions 
changing as : (1.084,0.97) -^ (1.16,0.96). 



^ After the calculations for this paper were completed the T2K and MINOS results on the first measure- 
ments of 0^3 using reactor muon neutrinos directed at distant neutrino detectors were announced [39j|40]. 
These results indicate that the likely range is ^^3 ~ 10° ± 5°. Since the search found fits with 6*^3 ~ 5° we 
anticipate no problem if the target is raised by 1-5 degrees. 
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NTH -I 




Eigenvalues{Au) 
Eigenvalues{Ag) 
Eigenvalues{Ap) 
Eigenvalues{Ae) 
Eigenvalues^Aq) 
Eigenvalues{AL) 
Ah, Ah 


0.010085 
0.013965 
0.015536 
0.015275 
0.007961 
0.005748 
3.007406 


0.031557 
0.037585 
0.057369 
0.039068 
0.024766 
0.038287 
1.661375 


0.032194 
0.038028 
0.057938 
0.040438 
0.025598 
0.039675 






NTH -2 




Eigenvalues{Au) 
Eigenvalues{Aj) 
Eigenvalues{Ap) 
Eigenvalues{Ae) 
Eigenvalues{AQ) 
Eigenvalues{AL) 
Ah, Ah 


0.069710 
0.066965 
0.078188 
0.086420 
0.063482 
0.077448 
9.309859 


1.007792 
1.024073 
1.049538 
0.999649 
1.058937 
1.075814 
4.858437 


1.033000 
1.053104 

1.074279 
1.027649 
1.085331 
1.104087 



Table 1: NTH12: Eigen values of the wave function correction matrices calculated using coupling 
values found in p/Tj without incorporating threshold corrections. Third generation values are in 
the first column. Note the significant threshold effects on the Higgs lines and the twofold degen- 
eracy among the first two generations as well as suppression of the third generation corrections 
due to cancelation between the large gauge and third generation yukawa contributions. Due to 
reduction of typical SO(IO) 16-plet yukawas both these features are modified in the fits found 
after incorporating threshold corrections at Mx- 

Typical values for the fermion line dressing coefficients due to yukawa couplings of the 
first two generations are of a few percent or smaller as seen in Table [1] where we exhibit 
eigenvalues of the fermion/Higgs line dressing matrices calculated for the two fits given in 
[17] which did not incorporate threshold corrections. Note that values for the Higgs lines 
can add up and achieve large values of up to several 100% ! If we switch off these corrections 
to Ah h then the changes in the fits are much smaller but still by no means negligible. These 
numbers make it clear that the light fermion lines and specially the light Higgs lines suffer 
very significant threshold corrections. Thus realistic GUT theories must face up to the 
task of specifying themselves sufficiently explicitly so that the threshold corrections may 
be calculated, tree level estimates are likely to be only rough pctures or even completely 
misleading. It is in this sense that we speak of Grand Unification "upended" . 

5 Realistic fits with threshold corrections included 

If we use our search programs to find fits after including the threshold effects we can impose 
strict perturbativity in the sense that no threshold correction may exceed 1 i.e 



\Afj^H,H\ < 1 



:iii 



The search programs [17] do find solutions (quite far from the examples of [17] in that some 
couplings, such as r) underwent major changes) which satisfied this constraint and still 
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provided accurate unification and accurate fits of the fermion mass data. However when 
one evaluates the rate of Baryon number violation (in the dominant B — > Meson + u 
channels) one finds (as in the case of the fits in [TTj) that typical solutions predict lifetimes 
of 10^^ years or smaller. This is 6 orders of magnitude below the current experimental 
limits |17] . This problem is an extension of that exhibited by the solutions found without 
threshold effects [T7]. 

Given the vast parameter space it is natural to ask whether there exist solutions where 
there is a suppression of the d = 5 operator mediated Baryon decay. Such fits can be 
found by instructing(via a x^ penalty for rapid baryon number violation) the 'amoeba' 
(i.e is the search engine of the downhill simplex method for nonlinear fitting [52]) to look 
for fits that have sufficiently low B-decay rates. An exhaustive statistical characterization 
of the parameter space and its possibilities requires the marshalling of considerable (su- 
per) computing resources which we shall eventually accomplish. For the moment we avoid 
over-determining an already excruciating search for viable fits by not also demanding such 
strict perturbativity, in the above sense, in addition to all the other demands of a sensible 
supersymmetric phenomenology that we must anyway impose. 

Since the complete programs for calculating B-decay rates (based upon the formulae 
provided in [251 IM]) are large and time consuming it would have slowed down both our 
search engine and the present computations too much to interface the complete Baryon 
violation programs (which include renormalizing some 447 variables from Mx to Mz, as well 
as a time consuming evaluation of the Baryon decay amplitudes) with our search programs 
at this still exploratory stage. Therefore we adopted the expedient of computing the the 
maximal absolute magnitude Max{0^^^) of the LLLL and RRRR coefficients in the d = 
5, AB 7^ effective superpotential for the NMSGUT [TH [T7] . For the solutions found so far 
this quantity was found to be typically of order 10^^^ to 10~^^ GeV^^ corresponding to the 
fast baryon decay rates ~ lO^^^yr"^ obtained. Thus a quick fix to the problem of limiting 
the B-decay rate while searching for accurate fermion fits is to limit (O is the dimensionless 
operator in units od \m/\\) Max{0^'^'') < 10~^GeV~^. This produced fits with proton 
lifetimes above 10'^^ yrs so we also relaxed the limit in some searches to just Max{0^'^^) < 
lO'^^GeV"^. When this condition was imposed simultaneously with the requirement of 
strict perturbativity flTTj) above, we were unable to find any accurate fits so far. On the 
other hand if one removes the condition of strict perturbativity but limits Max{0^^^) < 
10'"^^ Gev~^ or Max{0^^^) < 10'"^"^ GeV~^ then it was possible to find accurate fits that 
gave lifetimes in excess of 10^'^yrs. In addition to the above mentioned penalties we also 
required the fits to satisfy the following consistency/phenomenological constraints : 

• As already explained in detail in [H] the gauge unification RG fiow is constrained so 
that perturbation theory in the gauge coupling at unification remains valid, the uni- 
fication scale is less than Mpianck and the GUT threshold contributions to 0:3 (M^X) 
are in the right range [THl IT^ [T3] : 

- 22.0 < Ag = A{aG\Mx)) < 25 

3.0 >Ax = A{LogioMx) > -0.3 
-.017 < A3 = a3{Mz) < -.004 (12) 

• We constrain the \fi{Mz)\,\{AQ)ii\{i = 1,2,3) parameters to be smaller than 150 
TeV. Typically these parameters emerge in the range ~ 50 — 90 TeV while the 
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gaugino masses Mj are driven to the lower limits imposed (since it is the ratios 
\li{Mz)\/Mi, (Ao)/Mj which cointrol the efficacy of the large tan/3 corrections for our 
purposes. This is the price one must pay to correct the fermion yukawas to achievable 
values in the NMSGUT. Large values of Aq are well known to lead to deep charge and 
colour breaking (CCB) minimae [57] or unbounded from below (UFB) potentials |68]. 
However it is also established [51] that the metastable standard vacua that we are 
considering (with all mass squared parameters of charged or coloured or sneutrino 
scalar fields positive i.e at a local minimum which preserves colour, charge and R- 
parity) can well be stable on times scales (~ 10 Giga-years) of order the age of the 
universe. Thus individual cases need to be re-examined in detail before dismissing 
any otherwise viable fit out of hand. We consider the situation further in the next 
section. 



• In accordance with experimental constraints [51] we also constrain lightest chargino 
(essentially wino W"^ ) masses to be greater than 110 GeV. All the charged sfermions 
as well as the charged Higgs are constrained to lie above 110 GeV and the uncharged 
Higgs(/i°) above 105 GeV. 

• Since the susy threshold corrections to yd,s,h necessary for the survival of the NMSGUT 
as a viable theory of fermion masses depend on logarithms of ratios of soft susy break- 
ing parameters, our scenario is obviously incapable of predicting the all important susy 
breaking scale. On the other hand this is also a theory that counts providing exact 
unbroken R-parity down to the lowest scales (so that the LSP is stable and a good 
Dark Matter candidate) as one of its main virtues. Thus it behooves us to search for 
fits constrained by requiring the mass of the LSP(which is purely Bino due to the 
large value /i ~ 50Tey that emerges) in various ranges best motivated from a LSP 
dark matter scenario point of view, such as > 101 GeV^(range I), 5 — 50 GeV (Range 
II) and 50 — 101 GeV (Range III), to get an initial glimpse of whether and how the 
effective Dark Matter scale is linked to the pattern of superpartner masses. Thus in 
Tables [2] to [T9] we provide examples of fits found with the LSP/Bino mass in each of 
these ranges. In addition, in the Appendix we provide 8 additional examples of fits 
in these categories. We emphasize that these examples are only indicative and that a 
thorough investigation should have a finer grain in the LSP mass ranges chosen and 
must also incorporate loop corrections, specially for sfermion masses besides various 
improvements discussed below. It must also consider the statistics of ffis found by 
much larger scale systematic searches than those we have so far been able to mount : 
due to the scale of computer resources required. We note that the Type II and III fits 
are variations on those in I found by changing just the constraint on the LSP/Bino 
mass. One of the cases namely II-2 did not yield a satisfactorily exact fit in spite of 
extensive running, therefore it is omitted from the Appendix. 

In Tables of Type X-a {X = I, II, III) we give the complete set of NMSGUT parame- 
ters defined at the one loop unification scale M^ together with the values of the soft Susy 
breaking parameters {{mo,mi/2, Aq, B, Mjj ^j}) together with Supersymmetric parameter 
/i. Thus we have tried out only a N=l Supergravity GUT motivated scenario with re- 
laxation permitting Non universal Higgs masses(SUGRY-NUHM) (this is justified from a 
GUT point of view since the light doublets are a mixture of doublets from several sources 
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in different SO(IO) representations). The procedure followed for finding the fits, specially 
the use of Susy threshold correction to correct the yukawas i/d^s to values consistent with 
the 10-120 FM Higgs structure, as well as to raise the effective yb in the MSSM so as 
relieve the well known tension requiring this SM value of this parameter to lie several 
standard deviations above its experimental value in order to achieve b, r unification [M], is 
described in detail in the second paper of this series [17] . Besides these parameter values of 
the SUGRY-NUHM NMSGUT we also give the mass spectrum of superheavy heavy fields 
including the right handed neutrinos and the Type I and Type II neutrino seesaw masses 
as well as the unification parameters Ax,g,3 described in detail in [H]. 

In tables of Type X-b we have given the values of the target (i.e two loop RGE extrapoled 
Susy threshold corrected MSSM yukawas and Susy Weinberg operator coefficients) fermion 
parameters and their uncertainties (estimated a la [33]) together with the achieved values 
and pulls. The reader may check that the fits are all excellent with typical fractional errors 
0(0.1%). We remark that we found it tedious and meaningless to push our program to 
further narrow the fit in view of the large uncertainties and the numerous corrections still 
to be incorporated in our calculations (see below). We also give the eigenvalues of the GUT 
scale yukawa vertex threshold correction operators. We note that Tables X-a show that 
there is a significant lowering of the size of the SO(IO) fermion yukawas so that the universal 
gauge corrections dominate and make the corrections to all three generations almost equal 
specially when the lowering is pronounced. We give also the values of the the "Higgs 
fractions [SI [121 E] ««, ca crucial for determining the fermion mass formulae [121 [ISl 171 [T7]. 
These parameters are determined as a consequence of the GUT scale symmetry breaking 
and the fine tuning to preserve a light pair of MSSM Higgs doublets. They distill the 
infiuence of the SO(IO) GUT (and its spinorial clebsches determined appositely for this 
purpose in [5]) on the low energy fermion physics. The reader may use them together with 
the formulae given in [13] to check the fits even without entering into the details of our 
GUT scale mass spectra. We note that the first components of the a,, at were chosen real 
by convention [I^(see Appendix C). 

In Tables of type X-c values of the SM masses at Mz are compared with those of masses 
from the run down yukawas achieved in the NMSGUT both before and after large tan/3 
driven radiative corrections. Note that the central value of mb{Mz) = 2.9 GeV becomes 
prima facie acceptable in contrast to small Aq scenarios where the need for mb{Mz) > 3.1 
GeV, ie more than one standard deviation from the experimental central value, has been 
a principal source of tension and anxiety for small Ao models [6^ and should perhaps 
motivate acceptance of exploration of the large Aq parameter space which seems almost 
unexplored so far. 

In Tables of type X-d values of the soft supersymmetry breaking parameters which are 
the most crucial and remarkable output of this study since they tie the survival of the 
NMSGUT to a specific type of soft Susy spectrum with large /x, Ao, B and third generation 
sfermion masses generation in the 10-100 TeV range. Remarkably and in sharp contrast 
to received (small Aq^M'Jjjj) wisdom the third s-generation is much heavier than the first 
two generations, which however are themselves not very light except for the right chiral 
smuon/selectron , sstrange/sdown and scharm/sup right handed charged sfermions which 
can actually descend close to their experimental lower limits. In doing so they keep alive 
the effectiveness of the pure Bino LSP (and pure Wino lightest chargino and next to lightest 
Neutralinos) as candidate dark matter by providing co-annihilation channels of the sort a 
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light stau is often enlisted for in standard Susy GUT scenarios. Also remarkable is the 
friability of the familiar 1:2:7 ratio of the gaugino masses Mi, M2, M3 that we discover 
repeatedly in our calculations. This ratio is almost fixed in stone by one loop RGE and GUT 
mandated gaugino mass universality for small Aq (and often provokes baroque elaborations 
of the gauge invariance principle designed to avoid it). Here however the large influence of 
the Aq parameters can change the RGE flow to the point where the gluino is often lighter 
than the Winos which are sometimes as much thrice the mass of the Bino. Note also that 
except for the very heavy third generation the actual sizes of the sfermion trilinear couplings 
are rather modest since they are the product of the Aq parameters and the Yukawas. For 
the third sgeneration the trilinears are roughly the same size as the masses themselves thus 
preserving naturality. We remind the reader that a diagonal two loop RGE flow from Mx 
to Mz was used to determine these soft susy parameters via the Susy threshold effects since 
only the diagonal formulae were easily accessible and seemed justified in view of our limited 
expectations of overall accuracy of sfermion spectra which we have so far evaluated only at 
tree level. In view of the importance of the susy spectra we discuss in detail in Appendix 
B the features of the two loop RGE flow which result in the unconventional spectra noted 
above. 

Finally Tables of type X-e,f give Susy particle determined using two loop RGEs with 
and without generation mixing switched on. They are so similar as to justify the use of the 
diagonal values for estimating the Susy threshold corrections. For the case of the lightest 
sfermions however the corrections are sometimes as large as 10-30%. This again sounds 
a note of caution regarding the exact numerical values of the lighter sfermion masses we 
provide. However even after incorporation of Loop corrections in addition to these effects 
we certainly expect that the broad division into LHC discoverable particles lighter than say 
2.5 TeV (the LHC beam energy available per parton which sets the upper limit of discovery 
potential at the LHC) and those heavier to have some cogency regarding the limits of the 
discoverable. Thus we have provided an additional Table [20] that collects all superparticles 
with masses less than this limit. We also note that the Higgs masses were calculated 
using the 1-loop corrected electroweak symmetry breaking conditions and 1-loop effective 
potential using a subroutine [51] based on[7T]. The wary reader fearful of loop corrections 
destroying the whole scenario may be reassured at least on that count. 
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Parameter 


Value 


Field 


Masses 






\SU(?>),SU{2),Y] 


{Units ofW^GeV) 


Xx 


0.0464 


A[l,l,4 


645.64 


Xz 


0.0148 


5[6,2,5/3] 


0.3633 


/^ii/10-6 


0.0212 


C[8,2, 1] 


35.72,325.28,339.02 


/l22/10-^ 


0.0344 


^[3,2,7/3] 


35.03,349.70,362.74 


/i33 


0.0026 


^[3,2,1/3] 


0.58,26.33,26.33 


/ii/io-^ 


0.0781 -0.1368i 




28.661,393.61,441.49 


/12/IO-6 


-1.9955 - 0.0830i 


i^[l,l,2] 


6.15,6.15 


/l3/10-^ 


0.0580 + 0.0517i 




25.31,325.28 


/22/10-^ 


6.6036 - 4.9627i 


G[1,1,0] 


0.091,0.72,0.72 


/23/10-^ 


2.0080 + 2.3459i 




0.718,30.69,30.92 


/33/IO-3 


-1.0051 + 0.4427i 


/i[l,2,l] 


1.437,20.71,34.27 


^12/10-4 


0.0605 + 0.1232i 




541.46,563.22 


^13/10-^ 


-0.0460 + 1.8407i 


/[3, 1,10/3] 


1.26 


to/10-^ 


6.3251 + 5.7460i 


J[3, 1,4/3] 


1.387,14.31,14.31 


A/10-2 


-0.8601 - 1.4974i 




44.05,383.45 


ri 


-10.3248 + 2.5325i 


^[3,1,8/3] 


50.91,468.83 


P 


0.7042 - 2.2528i 


L[6, 1,2/3] 


24.18,752.08 


k 


0.0151 - 0.0805i 


M[6, 1,8/3] 


761.55 


c 


1.6200 + 0.5400i 


Ar[6, 1,4/3] 


757.79 


c 


1.0084 + 0.4594? 


0[1,3,2] 


1454.74 


m/lO^^GeV 


0.04 


P[3,3,2/3] 


14.50,1130.98 


mo/W^GeV 


-21.047e-^^''3(^) 


g[8,3,0] 


1.041 


1 


3.71 


i?[8, 1,0] 


0.40, 1.55 


7 


-2.8691 


5[1,3,0] 


1.7528 


X 


0.9397 + 0.6629i 


t[3, 1,2/3] 


1.15,19.66,47.99,78.65 


Ax 


1.16 




252.51,337.25,7050.35 


Ag 


4.863 


f/[3,3,4/3] 


1.480 


AaaiMz) 


-0.013 


\/[l,2,3] 


1.046 


{M'^yiO^^GeV} 


0.01,14.65,613.85 


l^[6,3,2/3] 


877.20 


{M^j/10-^^eV} 


0.4997,927.38,38856.88 


X[3,2,5/3] 


0.353,28.201,28.201 


M^{meV) 


2.17,7.63,42.34 


y[6,2,l/3] 


0.44 


{Evals[f]}/10-^ 


0.15,282.68,11844.25 


Z[8,l,2] 


1.54 


Soft parameters 


mi = -88.707 

2 


mo = 4198.698 


Ao = -1.1832 X 10^ 


at Mx 


/i = 9.4206 X 10^ 


B = -5.9399 X 10^ 


tan/3 = 50.0000 




M| = -7.1782 X 10^ 


M2j = -6.7789 X 10^ 


Rb. = 2.6935 


Max{\LABCD\, \Rabcd ) 


6.1149 X 10-23 GeV-i 




s/^ 



Table 2: I-l-a : Column 1 contains values of the NMSGUT-SUGRY-NUHM parameters at Mx 
derived from an accurate fit to all 18 fermion data and compatible with RG constraints. Unificaton 
parameters and mass spectrum of superheavy and superlight fields are also given. The values of 
lJ,{Mx),B{Mx) are determined by RG evolution from Mz to Mx of the values determined by 
the EWRSB conditions. 
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Parameter 


Target = Oi 


Uncert. = 5i 


Achieved = Oi 


Pull = (0, - 0,)/6, 


y«/io-^ 


2.098620 


0.801673 


2.098242 


-0.000472 


Wio-^ 


1.022984 


0.168792 


1.023751 


0.004546 


yt 


0.383855 


0.015354 


0.383875 


0.001336 


i/d/io-^ 


6.826388 


3.979784 


6.843968 


0.004417 


i/s/io-3 


1.294453 


0.610982 


1.301389 


0.011351 


Vb 


0.460136 


0.238811 


0.460919 


0.003278 


l/e/lO-4 


1.201757 


0.180264 


1.202499 


0.004117 


^m/IO-^ 


2.464838 


0.369726 


2.468141 


0.008935 


l/r 


0.527874 


0.100296 


0.523894 


-0.039688 


sin 9 12 


0.2210 


0.001600 


0.2210 


-0.0005 


sin^?3/10"4 


29.1027 


5.000000 


29.1272 


0.0049 


sin^ls/lO-^ 


34.2424 


1.300000 


34.2402 


-0.0017 


5« 


60.0207 


14.000000 


60.0550 


0.0024 


{m\2)llQ-'{eVf 


5.3580 


0.567948 


5.3569 


-0.0019 


{ml,)/lQ-\eVf 


1.7330 


0.346597 


1.7347 


0.0050 


sin^ ^^2 


0.2887 


0.057748 


0.2883 


-0.0080 


sin^ e^^ 


0.4620 


0.138613 


0.4639 


0.0132 


e{^{degrees) 


3.7 


3.7 


4.57 




Eigenvalues{Au) 


0.048796 


0.048800 


0.048805 




Eigenvalues{Ag) 


0.046619 


0.046623 


0.046629 




Eigenvalues{Ap) 


0.052548 


0.052552 


0.052558 




Eigenvalues{Ae) 


0.059078 


0.059082 


0.059087 




Eigenvalues{AQ) 


0.046281 


0.046284 


0.046289 




Eigenvalues{AL) 


0.054386 


0.054390 


0.054394 




Ah, Ah 


71.935369 


62.415435 






ai 


0.8302 + O.OOOOi 


ai 


0.8899 + O.OOOOi 




a2 


0.0644 + 0.0157i 


(52 


0.0516 + 0.0644i 




as 


-0.0425 - 0.0442i 


«3 


-0.0571 - 0.0204i 




a4 


-0.3970 + 0.1618i 


(54 


0.3186 - 0.0237i 




as 


0.1534 + 0.0818i 


"5 


0.0709 + 0.014n 




ae 


0.1116 -0.2759i 


ae 


0.1187 -0.2760i 





Table 3: I-l-b: Fit with xx = \Jj2}liiOi - Oi^/df = 0.0464. Target values, at Mx of the 
fermion yukawa couphngs and mixing parameters, together with the estimated uncertainties, 
achieved values and pulls. The eigenvalues of the wavefunction renormalization increment matrices 
Aj for fermion lines and the factors for Higgs lines are given, assuming the external Higgs is 10- 
plet dominated. Note the close similarity of the eigenvalues which suggests that the small values 
of the SO (10) yukwas utilized when threshold corrections are in play lead to gauge dominated 
corrections which are the same for all three generations. The Higgs fractions ai, di which control 
the MSSM fermion yukawa couplings are also given. Notice the dominance of the first components 
ai, di consistently with the assumption made. Right handed neutrino threshold effects have been 
ignored. We have truncated numbers for display although all calculations are done at double 
precision. 
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Parameter 


SM{Mz) 


m^^'^iMz) 


^MSSM = (^ + Am)^^^(Mz) 


md/10-3 


2.90000 


0.63332 


2.92341 


m,/10-3 


55.00000 


12.04268 


55.60296 


rrib 


3.00000 


3.07778 


3.00662 


me/10-3 


0.48657 


0.47249 


0.48691 


m^ 


0.10272 


0.09694 


0.10270 


ITir 


1.74624 


1.73514 


1.73675 


m„/10-3 


1.27000 


1.09930 


1.27007 


rric 


0.61900 


0.53635 


0.61968 


ITit 


172.50000 


149.05874 


172.47396 



Table 4: I-l-c: Values of standard model fermion masses in GeV at Mz compared with the masses 
obtained from values of GUT derived yukawa couplings run down from M^ to Mz both before 

and after threshold corrections. Fit with xz = \/lfi=iirnf^^^ - m^'^Y /{m^^^'^Y = 0.0148. 



Parameter 


Value 


Parameter 


Value 


Ml 


124.94 


M,, 


5392.74 


M2 


328.39 


M^. 


5392.01 


M3 


569.92 


Mn, 


17882.28 


^\ 


1101.73 


40(0 

^11 


-75168.69 


^I. 


165.22 


.0(0 
^22 


-75083.12 


^3 


11100.23 


40(0 

^33 


-47868.83 


Ml, 


6400.13 


.0(«) 

^11 


-86227.11 


Ml, 


6353.87 


.0(«) 
^22 


-86226.53 


Ml, 


10210.78 


.0(«) 
^33 


-43721.21 


% 


2917.32 


40(d) 


-75501.89 


Mj 

d2 


2916.57 


^22 


-75501.25 


^3 


26552.10 


40(d) 
^33 


-32031.46 


Mq, 


4928.60 


tan/3 


50.00 


Mq, 


4927.98 


KMz) 


76666.76 


Mq, 


22679.30 


B{Mz) 


9.6672 X 10^ 


M^- 


-6.0301 X 10^ 


Mh 


-6.3249 X 10^ 



Table 5: I-l-d: Values (GeV) in of the soft Susy parameters at Mz (evolved from the soft 
SUGRY-NUHM parameters at Mx)- The values of soft Susy parameters at Mz determine the 
Susy threshold corrections to the fermion yukawas. The matching of run down fermion yukawas 
in the MSSM to the SM parameters determines soft SUGRY parameters at Mx- Note the heavier 
third sgeneration. The values of ^{Mz) and the corresponding soft susy parameter B{Mz) = 
m^sin 2/3/2 are determined by imposing electroweak symmetry breaking conditions. mA is the 
mass of the CP odd scalar in the in the Doublet Higgs. The sign of /u is assumed positive. 



Field 


Mass{GeV) 


Ma 


569.92 


M^± 


328.38,76666.85 


M,o 


124.94, 328.38, 76666.81, 76666.82 


M, 


6399.784, 6353.520, 10210.565 


Me 


1102.65, 6400.31, 160.56, 6354.33, 10105.12, 11196.70 


Mu 


4928.29, 5392.61, 4927.63, 5391.92, 17876.33, 22684.87 


M^, 


2917.43, 4928.98, 2916.66, 4928.38, 22663.07, 26566.04 


Ma 


219898.08 


Mh± 


219898.09 


Mho 


219898.06 


Mho 


111.45 



Table 6: I-l-e: Spectra of super symmetric partners calculated ignoring generation mixing effects. 
Inclusion of such effects changes the spectra only marginally. Due to the large values of /i, B, Aq. 
The LSP and light chargino are essentially pure Bino and Wino(W±). The light gauginos and light 
Higgs h^, are accompanied by a light smuon and sometimes selectron. The rest of the sfermions 
have multi-TeV masses. The mini-split supersymmetry spectrum and large fi, Aq parameters help 
avoid problems with FCNC and CCB/UFB instability [53]. The sfermion masses are ordered by 
generation not magnitude. This is useful in understanding the spectrum calculated including 
generation mixing effects. The mass of the Higgs particles(Myi,M/jO,Mj:^±) are calculated by 
incorporating one loop contributions to the Electroweak symmetry breaking i.e to the effective 
potential[7Il[56l[I7]. 



Field 


Mass{GeV) 


Ma 


570.03 


M,± 


328.43, 76664.05 


M,o 


124.95, 328.43, 76664.02, 76664.02 


M, 


6353.56,6399.82,10211.527 


Mg 


133.31, 1102.50, 6353.86, 6400.35, 10109.89, 11195.65 


M, 


4926.11, 4928.47, 5392.09, 5392.79, 17873.85, 22683.48 


M,: 


2916.90, 2917.69, 4926.85, 4929.16, 22661.71, 26564.36 


Ma 


219973.32 


Mh± 


219973.33 


Mho 


219973.30 


MhO 


111.74 



Table 7: I-l-f : Spectra of supersymmetric partners calculated including generation mixing effects. 
Inclusion of such effects changes the spectra only marginally. Due to the large values of fi,B,AQ 
the LSP and light charginos are essentially pure Bino and Wmo{W±). Note that the ordering 
of the eigenvalues in this table follows their magnitudes, comparison with the previous table is 
necessary to identify the sfermions 
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Parameter 


Value 


Field 


Masses 






[SU{3),SU{2),Y] 


{Units oflO^^GeV) 


Xx 


0.0133 


A[l,l,4 


622.89 


Xz 


0.0386 


S[6,2,5/3] 


0.4399 


/ill/10-6 


0.0224 


C[8,2, 1] 


34.59,332.64,344.97 


/i22/10-^ 


0.0283 


£'[3,2,7/3] 


33.75,322.29,340.48 


/i33 


0.0026 


^[3,2,1/3] 


0.70,26.06,26.06 


/11/IO-6 


0.0782 -0.1347i 




27.899,377.48,429.89 


/12/IO-6 


-2.0477 - 0.0905i 


F[l,l,2] 


6.08,6.08 


/13/IO-5 


0.0543 + 0.0544i 




24.31,313.95 


/22/IO-5 


6.4657 -4.723n 


^[1,1,0] 


0.110,0.80,0.80 


/23/10-^ 


1.9840 + 2.3506i 




0.872,33.23,33.50 


/33/10-=^ 


-1.0229 + 0.4338i 


h[l,2,l] 


1.748,20.51,33.32 


^12/10-^ 


0.0560 + 0.1257i 




510.60,538.80 


^13/10-^ 


-0.1809 + 1.6222i 


J[3, 1,10/3] 


1.52 


^23/10-^ 


6.2005 + 5.7438i 


J[3, 1,4/3] 


1.661,14.22,14.22 


A/10-2 


-1.2230 -1.7024i 




42.85,375.68 


V 


-10.0887 + 2.4493i 


^[3,1,8/3] 


49.33,459.01 


P 


0.6745- 2. 1659i 


L[6, 1,2/3] 


23.85,735.63 


k 


0.0165 - 0.0849i 


M[6, 1,8/3] 


746.01 


C 


1.5173 + 0.4434i 


Ar[6, 1,4/3] 


739.68 


c 


1.0209 + 0.4394i 


0[1,3,2] 


1396.55 


m/W^GeV 


0.05 


P[3,3,2/3] 


12.83, 1088.74 


nio/lQ^^GeV 


-20.758e-*^''^(^) 


g[8,3,0] 


1.262 


7 


3.67 


i?[8, 1,0] 


0.49,1.85 


7 


-2.9668 


5[1,3,0] 


2.1051 


X 


0.9290 + 0.6512i 


t[3, 1,2/3] 


1.39,19.29,46.66,79.08 


Ax 


1.15 




249.60,335.17,6705.69 


Ag 


5.955 


^[3,3,4/3] 


1.780 


Aa3{Mz) 


-0.009 


\/[l,2,3] 


1.258 


{M'^yW^GeV} 


0.00,16.54,687.11 


W[6,3,2/3] 


854.17 


{M^j/10-^^eV} 


0.2400,912.69,37917.07 


X[3,2,5/3] 


0.424,27.998,27.998 


M^{meV) 


3.44,8.12,42.61 


F[6,2,l/3] 


0.54 


{Evals[f]}/10-^ 


0.08,287.79,11955.91 


Z[8,l,2] 


1.84 


Soft parameters 


mi = -124.831 

2 


mo = 2371.139 


Ao = -6.1911 X 10^ 


at Mx 


^ = 5.8956 X 10^ 


B = -1.9501 X 10^ 


tan/3 = 50.0000 




M| = -3.0748 X 10^ 


M^j = -2.9538 X 10^ 


Rb. = 2.9200 


Max{\LABCD\, \Rabcd ) 


6.7561 X 10^23 Qgy-i 




SM 



Table 8: Il-l-a: Column 1 contains values of the NMSGUT-SUGRY-NUHM parameters at Mx 
derived from an accurate fit to all 18 fermion data and compatible with RG constraints. See 
caption to Table 2 for explanation. 
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Parameter 


Target = Oi 


Uncert. = 5i 


Achieved = Oi 


Pull = (0, - 0,)/6, 


?/./io-6 


2.141027 


0.817872 


2.141050 


0.000028 


l/c/10-' 


1.043728 


0.172215 


1.043795 


0.000391 


Vt 


0.386577 


0.015463 


0.386580 


0.000221 


l/d/10-^ 


6.891765 


4.017899 


6.919378 


0.006873 


l/s/10-3 


1.311160 


0.618867 


1.311543 


0.000620 


Vb 


0.451549 


0.234354 


0.453520 


0.008411 


1/e/lO-^ 


1.221705 


0.183256 


1.221655 


-0.000272 


1/m/10-^ 


2.580846 


0.387127 


2.580746 


-0.000257 


yr 


0.525812 


0.099904 


0.525048 


-0.007648 


sin ^f 2 


0.2210 


0.001600 


0.2210 


-0.0001 


sin^?3/10-^ 


29.1374 


5.000000 


29.1407 


0.0007 


sin^lg/lO-^ 


34.2832 


1.300000 


34.2832 


0.0000 


5« 


60.0207 


14.000000 


60.0102 


-0.0008 


{m\,)/lQ-\eVf 


5.4100 


0.573459 


5.4100 


0.0001 


{ml,)/lQ-\eVf 


1.7497 


0.349943 


1.7497 


-0.0001 


sin^ ^^2 


0.2887 


0.057750 


0.2887 


0.0000 


sin^ ^^3 


0.4621 


0.138631 


0.4621 


-0.0001 


^f3(rfe^rees) 


3.7 


3.7 


8.18 




Eigenvalues{Ay) 


0.048603 


0.048607 


0.048613 




Eigenvalues{Ag) 


0.046536 


0.046540 


0.046546 




Eigenvalues{Ap) 


0.052666 


0.052670 


0.052675 




Eigenvalues{Ae) 


0.058866 


0.058870 


0.058876 




Eigenvalues{AQ) 


0.046115 


0.046119 


0.046124 




Eigenvalues{AL) 


0.054312 


0.054316 


0.054320 




An, Ah 


70.962744 


62.481215 






tti 


0.8263 + O.OOOOi 


tti 


0.8795 + O.OOOOi 




^2 


0.0646 + 0.0197i 


(52 


0.0509 + 0.0709i 




as 


-0.0437 - 0.0505i 


as 


-0.0594 - 0.0234i 




q;4 


-0.3966 + 0.1515i 


(54 


0.3314 - 0.0234i 




as 


0.1606 + 0.0812i 


"5 


0.0782 + 0.012n 




"6 


0.1088 - 0.2894i 


"6 


0.1174 -0.2904i 





Table 9: Il-l-b: Fit with xx = \Jj2}li{Oi - OiY/df = 0.0133. Target values, at Mx of the 
fermion yukawa couphngs and mixing parameters, together with the estimated uncertainties, 
achieved values and pulls. See caption to Table 3 for explanation. 
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Parameter 


SM{Mz) 


m'^^'^iMz) 


^MSSM _ (^ ^ Am)^^^(Mz) 


md/lO"^ 


2.90000 


0.64294 


2.90511 


m,/10^3 


55.00000 


12.18662 


54.90038 


rrib 


3.00000 


3.05467 


3.02134 


me/10-3 


0.48657 


0.48198 


0.48770 


rrif. 


0.10272 


0.10177 


0.10289 


m-r 


1.74624 


1.74369 


1.74266 


m„/10-3 


1.27000 


1.11799 


1.23689 


rric 


0.61900 


0.54503 


0.60300 


mt 


172.50000 


149.51960 


172.26015 



Table 10: II-l-c: Values of standard model fermion masses in GeV at Mz compared with the 
masses obtained from values of GUT derived yukawa couplings run down from M^ to Mz both 

before and after threshold corrections. Fit with xz = \/J2i=i{''^i 
0.0386. 



,MSSM 



m 



SM\2 



/{m. 



MSSM\2 



Parameter 


Value 


Parameter 


Value 


Ml 


35.00 


M^, 


2490.81 


M2 


113.01 


M-n, 


2490.90 


M3 


83.09 


Mn, 


16928.60 


^\ 


1045.62 


Ml) 
^11 


-39522.44 


^% 


1148.28 


40(0 

^22 


-39473.12 


^3 


14376.75 


40(0 


-25104.91 


Mu 


3699.46 


40W 
^11 


-44729.39 


Mu 


3714.57 


4OH 
^22 


-44729.08 


Ml, 


10786.67 


^33 


-22492.46 


% 


320.40 


.0(d) 

^11 


-39448.33 


M-- 

d2 


323.59 


^22 


-39447.98 


^Ja 


19724.24 


40(d) 
^33 


-16948.89 


Mq, 


2422.85 


tan/3 


50.00 


Mq, 


2423.11 


KMz) 


48104.14 


Mq, 


18383.64 


B{Mz) 


3.0799 X 10^ 


M% 


-2.3660 X 10^ 


Ml 


-2.5187 X 10^ 



Table 11: Il-l-d: Values (GeV) in of the soft Susy parameters at Mz (evolved from the soft 
SUGRY-NUHM parameters at Mx)- See caption to Table 5 for explanation. 
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Field 


^^055(^6^) 


Ma 


83.09 


M,± 


113.01,48104.28 


M^o 


35.00, 113.01, 48104.23, 48104.24 


M, 


3698.868, 3713.977, 10786.469 


Mg 


1046.59, 3699.78, 1147.01, 3715.55, 10777.59, 14383.72 


Mu 


2422.23, 2490.54, 2422.10, 2491.00, 16922.19, 18390.64 


M,- 


321.46, 2423.62, 324.40, 2423.91, 18355.16, 19750.86 


M^ 


124120.40 


Mh± 


124120.43 


Mho 


124120.39 


M,o 


130.18 



Table 12: Il-l-e: Spectra of supersyninietric partners calculated ignoring generation mixing 
effects. See caption to Table 6 for explanation. 



Field 


Mass{GeV) 


Ma 


83.15 


M,± 


113.03,48102.54 


M,a 


35.01, 113.03, 48102.49, 48102.50 


M, 


3698.88,3713.99,10786.642 


Me 


1046.52, 1144.84, 3699.80, 3715.67, 10778.23, 14383.76 


Mu 


2422.27, 2438.70, 2490.58, 2491.10, 16922.33, 18388.56 


M,: 


321.58, 324.46, 2423.67, 2440.53, 18353.16, 19750.82 


Ma 


124163.88 


Mh± 


124163.90 


Mho 


124163.87 


Mho 


130.20 



Table 13: Il-l-f: Spectra of supersymmetric partners calculated including generation mixing 
effects. See caption to Table 7 for explanation. 
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Parameter 


Value 


Field 


Masses 






[5[/(3),5t/(2),r] 


{Units ofW^GeV) 


Xx 


0.0038 


A[l,l,4 


756.12 


Xz 


0.0030 


S[6,2,5/3] 


0.4240 


/ill/10-6 


0.0211 


C[8,2, 1] 


39.13,387.00,401.03 


/i22/10-^ 


0.0329 


£'[3,2,7/3] 


38.83,385.83,402.57 


/i33 


0.0026 


E[3,2,l/3] 


0.68,29.71,29.71 


/11/IO-6 


0.0774 -0.1363i 




31.591,443.85,505.56 


/12/IO-6 


-1.9810 -0.0598i 


i"[l,l,2] 


7.05,7.05 


/13/IO-5 


0.0578 + 0.0556i 




28.49,380.72 


/22/IO-5 


6.4063 - 4.9278i 


^[1,1,0] 


0.105,0.82,0.82 


/23/10-^ 


2.0067 + 2.2917i 




0.836,35.65,35.90 


/33/10-=^ 


-0.9930 + 0.4524i 


/i[l,2,l] 


1.695,23.43,37.99 


^12/10-^ 


0.0610 + 0.1236i 




611.35,639.58 


^13/10-^ 


-0.0066 + 1.7851i 


J[3, 1,10/3] 


1.45 


^23/10-^ 


6.4307 + 5.7410i 


J[3, 1,4/3] 


1.588,16.35,16.35 


A/10-2 


-1.0000 -1.4449i 




49.47,451.05 


V 


-10.5299 + 2.6810i 


^[3,1,8/3] 


57.16,553.55 


P 


0.6944 - 2.2125i 


L[6, 1,2/3] 


27.44,883.40 


k 


0.0166 - 0.0797i 


M[6, 1,8/3] 


896.34 


C 


1.5479 + 0.5373i 


Ar[6, 1,4/3] 


888.19 


c 


1.0081 + 0.4175i 


0[1,3,2] 


1667.08 


m/W^GeV 


0.04 


P[3,3,2/3] 


15.50, 1292.76 


nio/lQ^^GeV 


-23.876e-^^''^(^) 


g[8,3,0] 


1.199 


7 


3.69 


i?[8, 1,0] 


0.47,1.78 


7 


-2.8264 


5[1,3,0] 


2.0164 


X 


0.9272 + 0.6601i 


t[3, 1,2/3] 


1.33,22.24,53.87,86.47 


Ax 


1.22 




286.53,380.42,8307.87 


Ag 


5.014 


^[3,3,4/3] 


1.701 


Aa3{Mz) 


-0.009 


\/[l,2,3] 


1.212 


{M'^yW^GeV} 


0.01,16.26,692.32 


iy[6,3,2/3] 


1005.11 


{M^j/10-^^eV} 


0.2849,723.96,30818.79 


X[3,2,5/3] 


0.409,31.922,31.922 


M^{meV) 


2.09,7.64,42.49 


F[6,2,l/3] 


0.52 


{Evals[f]}/10-^ 


0.11,276.17,11756.40 


Z[8,l,2] 


1.77 


Soft parameters 


mi = -107.371 

2 


mo = 3585.150 


Ao = -1.0442 X 10^ 


at Mx 


/i = 8.4934 X 10^ 


B = -4.7848 X 10^ 


tan/3 = 50.0000 




M| = -5.8077 X 10^ 


M^j = -5.5072 X 10^ 


Rb. = 2.7107 


Max{\LABCD\, \Rabcd ) 


5.8199 X 10^23 Qgy-l 




S^l 



Table 14: Ill-l-a : Column 1 contains values of the NMSGUT-SUGRY-NUHM parameters at 
Mx derived from an accurate fit to all 18 fermion data and compatible with RG constraints. See 
caption to Table 1 for explanation. 
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Parameter 


Target = Oi 


Uncert. = 5i 


Achieved = Oi 


Pull = (0, - 0,)/6, 


?/./io-6 


2.116326 


0.808437 


2.116304 


-0.000027 


l/c/10-' 


1.031618 


0.170217 


1.031400 


-0.001282 


Vt 


0.386527 


0.015461 


0.386531 


0.000285 


l/d/10-^ 


6.872020 


4.006388 


6.880189 


0.002039 


l/s/10-3 


1.303124 


0.615075 


1.301971 


-0.001874 


Vb 


0.468198 


0.242995 


0.468608 


0.001688 


1/e/lO-^ 


1.209876 


0.181481 


1.209817 


-0.000328 


1/m/10-^ 


2.485105 


0.372766 


2.485397 


0.000783 


yr 


0.532529 


0.101181 


0.532449 


-0.000788 


sin ^f 2 


0.2210 


0.001600 


0.2210 


-0.0001 


sin^?3/10-^ 


29.0367 


5.000000 


29.0367 


0.0000 


sin^lg/lO-^ 


34.1648 


1.300000 


34.1652 


0.0003 


5« 


60.0207 


14.000000 


60.0123 


-0.0006 


{m\,)/lQ-\eVf 


5.3985 


0.572236 


5.3985 


0.0001 


{ml,)/lQ-\eVf 


1.7468 


0.349356 


1.7467 


-0.0001 


sin^ ^^2 


0.2886 


0.057728 


0.2886 


0.0001 


sin^ ^^3 


0.4617 


0.138503 


0.4616 


-0.0006 


^f3(rfe^rees) 


3.7 


3.7 


5.38 




Eigenvalues{Ay) 


0.051945 


0.051949 


0.051955 




Eigenvalues{Ag) 


0.049907 


0.049911 


0.049917 




Eigenvalues{Ap) 


0.056963 


0.056967 


0.056973 




Eigenvalues{Ae) 


0.063077 


0.063081 


0.063087 




Eigenvalues{AQ) 


0.048996 


0.049000 


0.049006 




Eigenvalues{AL) 


0.058090 


0.058094 


0.058099 




An, Ah 


72.623699 


62.630532 






tti 


0.8274 + O.OOOOi 


tti 


0.8895 + O.OOOOi 




^2 


0.0650 + 0.0179i 


02 


0.0526 + 0.0664i 




as 


-0.0415 - 0.0457? 


as 


-0.0554 - 0.0217i 




q;4 


-0.4028 + 0.1724i 


(54 


0.3209 - 0.0356i 




as 


0.1511 + 0.0810i 


"5 


0.0705 + 0.0136i 




"6 


0.1099 - 0.2714i 


"6 


0.1150 - 0.2746i 





Table 15: Ill-l-b: Fit with xx = \/Y.i=iiOi - Oi)^/5f = 0.0038. Target values, at Mx of 
the fermion yukawa couplings and mixing parameters, together with the estimated uncertainties, 
achieved values and pulls. See caption to Table 2 for explanation. 
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Parameter 


SM{Mz) 


m'^^'^iMz) 


^MSSM _ (^ ^ Amf^'^iMz) 


md/lO"^ 


2.90000 


0.63252 


2.90621 


m,/10^3 


55.00000 


11.96937 


55.00879 


rrib 


3.00000 


3.09207 


3.00111 


me/10-3 


0.48657 


0.47228 


0.48641 


rrif. 


0.10272 


0.09698 


0.10252 


m-r 


1.74624 


1.74550 


1.74567 


m„/10-3 


1.27000 


1.10652 


1.27009 


rric 


0.61900 


0.53927 


0.61899 


rat 


172.50000 


149.30372 


172.48084 



Table 16: III-l-c: Values of standard model fermion masses in GeV at Mz compared with the 
masses obtained from values of GUT derived yukawa couplings run down from M^ to Mz both 

before and after threshold corrections. Fit with xz = \/J2i=i{''^i 
0.0030. 



,MSSM 



m 



SM\2 



/{m. 



MSSM\2 



Parameter 


Value 


Parameter 


Value 


Ml 


100.07 


M^, 


4539.23 


M2 


270.77 


M-n, 


4538.61 


M3 


428.43 


Mn, 


17182.81 


^\ 


906.34 


Ml) 
^11 


-65947.70 


^% 


142.19 


40(0 

^22 


-65871.94 


^3 


11976.39 


40(0 


-41513.30 


Mu 


5564.41 


40W 
^11 


-75901.65 


Mu 


5528.49 


4OH 
^22 


-75901.14 


Ml, 


10194.21 


^33 


-38216.99 


% 


2291.12 


.0(d) 

^11 


-66155.83 


M-- 

d2 


2290.56 


^22 


-66155.26 


^Ja 


24659.75 


40(d) 
^33 


-27514.18 


Mq, 


4205.32 


tan/3 


50.00 


Mq, 


4204.83 


KMz) 


68576.15 


Mq, 


21289.22 


B{Mz) 


7.4868 X 10^ 


M% 


-4.7966 X 10^ 


Ml 


-5.0786 X 10^ 



Table 17: Ill-l-d: Values (GeV) in of the soft Susy parameters at Mz (evolved from the soft 
SUGRY-NUHM parameters at Mx)- See caption to Table 4 for explanation. 
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Field 


MassiGeV) 


Ma 


428.43 


M,± 


270.76,68576.25 


M^o 


100.07, 270.76, 68576.21, 68576.22 


M, 


5564.015, 5528.094, 10193.999 


Mg 


907.46, 5564.62, 136.00, 5529.05, 10149.71, 12014.32 


Mu 


4204.96, 4539.08, 4204.40, 4538.53, 17176.99, 21294.85 


M;i, 


2291.27, 4205.76, 2290.68, 4205.29, 21272.06, 24674.64 


Ma 


193516.73 


Mh± 


193516.75 


Mho 


193516.72 


M,o 


117.18 



Table 18: Ill-l-e: Spectra of supersymmetric partners calculated ignoring generation mixing 
effects. See caption to Table 6 for explanation. 



Field 


Mass{GeV) 


Ma 


428.53 


M,± 


270.80,68573.70 


M,o 


100.09, 270.80, 68573.66, 68573.67 


M, 


5528.13,5564.05,10194.739 


Me 


105.36, 907.31, 5528.77, 5564.65, 10152.37, 12014.38 


Mu 


4205.11, 4205.54, 4538.67, 4539.24, 17175.20, 21293.36 


M,: 


2290.89, 2291.51, 4205.91, 4206.41, 21270.60, 24673.37 


Ma 


193585.06 


Mh± 


193585.07 


Mho 


193585.04 


Mho 


117.39 



Table 19: Ill-l-f : Spectra of supersymmetric partners calculated including generation mixing 
effects. See caption to Table 7 for explanation. 
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6 Structural features and Phenomenology 

As mentioned the threshold corrections to the Higgs hnes in a hght fermion- hght higgs 
vertex can obtain a thick wave function dressing from GUT scale particles leading to an 
amplification of the effective SO(IO) yukawa couplings and making much weaker SO(IO) 
couplings capable of fitting the fermion data. Moreover the threshold corrections of the 
light fermion yukawa couplings are highly nonlinear, and as such there is no reason to 
expect that the constraint yt, — i/r — ys — y^i found on the basis of the tree level coupling of 
the 10,120-plets to the 16-plets of fermionsjHl H2] continues to be effective. The influence 
of these large corrections can be seen in the value of Rbr/siJ, = -^^5^ given in the tables. We 
see that only in cases 1 — 2, III — 2, /// — 4 is Rbr/sfi still approximately unity while in the 
rest it is typically ~ 3: which is the magnitude of the ratio when all threshold corrections 
to SM values are ignored and the couplings are run up using MSSM RGEs. This shows 
that the tree level constraint can be evaded; and with profit since the magnitudes are no 
longer mismatched. Note that the cases with Rbr/sn ~ 1 typically have A^^-^ smaller by a 
factor of about 5 than the cases where Rtr/sfi ~ 3. 

The ultra small values of the 126 couplings ensure that they make little difference to 
the 2nd and 3rd generation charged fermion yukawas but ensure that M,yc are light and 
thus Type I neutrino masses are viable^. Note that for the same reason in all fits the 
Type II neutrino masses are completely negligible. The other superpotential couplings are 
unremarkable except that rj is somewhat large. However one should recall that it occurs in 
the superpotential divided by 4! = 24. Actual coefficients of the radiative corrections will 
tj^ically be powers of (~ ^). 

As explained in Section 1,2 the unification scale Mx = M'^10^^ , (typically ~ lO^''''^ 
GeV) is identified with the mass of the X[3, 2, ±|] gauge sub-multiplet ( exchange of which 
gives rise to d=6 operators for B decay ) and determines the scale parameter m of the 
superpotential via 

UMgl6.25+Ax 

m = — — = GeV (13) 

Mx 

where 

Mx = ^r\/4|a + a)|2 + 2|p + (I;|2 

and 



is the corrected SO(IO) gauge coupling. We see that the unification scale is generally 
elevated( 10^^-^ — 10^^ GeV) but the SO(IO) coupling at unification is still perturbative 
though sometimes only marginally so. 

The right handed neutrino masses are important for Leptogenesis and for lepton RGE 
flows at intermediate scales. We find M,yc is generically in the range 10^ — 10^^ GeV (with 
normal hierarchy); which is also the preferred range for Leptogenesis. It is determined by 
the necessarily (for viable neutrino spectra) ultra small 126 couplings / ~ 10^^ — 10^'^. 
It is interesting to note that the threshold corrections may also weaken the influence of 
the right handed neutrino thresholds on yukawa unification. The interplay of the GUT 
scale corrections and the M^^c thresholds will be interesting to evaluate, specially since the 
latter are known [65] to lead to tension for yukawa unification. It may be that as in the 
case of the tension regarding the value of mb{Mz) acceptable for yukawa unification [61]. 



which is reheved in our model by the threshold corrections at Ms = Mz, so also GUT scale 
threshold corrections may help with relieving the effect of right handed neutrino thresholds. 

The super heavy masses lie in the range 10^^ — 10^" GeV with a few multiplets some- 
times having an uncomfortably large mass even greater than the Planck mass. Overall the 
unification parameter and the spectra of the Fits with the large Ajjjj ~ 10^ and Rbr/sn ~ 3 
seem more palatable than those for Case 1-2 with Rbr/sfi ~ 1, 'Ah,h ~ 10. 

Let us turn next to the conjectured Soft Spectra determined by requiring EWRSB 
at tan /3 ~ 50 and large threshold corrections to lower t/^ ., by a factor of 5 or so. The 
required values of the |/^|, \Aq\ parameters turn out to be so huge (~ 10^ TeV) that they 
will incite controversy driven by concerns regarding deep CCB minima. For the moment 
we take the pragmatic attitude that we have checked the local stability by ensuing the 
positivity of all scalar mass squared parameters. The stability against tunneling to CCB 
minima on cosmological time scales calls for further investigation after loop corrections to 
scalar masses have been included. However the literature [Ml EH] supports the pragmatic 
attitude we adopt regarding the viability of metastable minima. The seminal and clear 
investigations of [51] regarding met a- stability in the parameter region of ultra large /i, Ao 
are so encouraging that we cannot resist quoting them verbatim. Firstly they note that 
(our interpolations in square brackets) " the height of the barrier separating the [metastable] 
minimum from the CCB minimum is roughly proportional to l/y^in' where ymm is the 
smallest Yukawa coupling associated with the fields that acquire non-zero vev in the CCB 
minimum. The corresponding tunnelling rates are greatly suppressed for small y. " Thus 
since it is only the third generation of matter sfermion fields that have appreciable yukawa 
couplings the violation of the CCB and UFB bounds [671 EH] is not likely to be a matter of 
concern for the first two generations. Moreover they note "i/ie most stringent constraints 
come from the small tan/3 region, where the top Yukawa coupling is larger" . Whereas we 
are in the large tan/3 ~ 50 regime. 

The investigations of [51] focussed on the region |/i|, \Aq\ < 4 TeV. Their findings con- 
firmed that in this region "t/ie larger the trilinear coupling the more dangerous is the corre- 
sponding CCB minimum" . However they had the prescience to realize that understanding 
the behaviour in the regions with much larger |/i|, \Aq\ would illuminate the dynamics of 
CCB and tunnelling and clarify the operation of decoupling arguments which seem vio- 
lated by the above tradecraft maxim but are always crucial to establish an intuitive grip 
on field theoretic dynamics. Thus they note "zt is instructive to examine what happens to 
the tunnelling probability in the limit of very large fi and At (and large enough squark mass 
terms to ensure the existence of the [metastable] minimum) . In that limit, as the CCB 
minimum moves away from the SML [standard] minimum, the barrier separating the two 
becomes thicker, and the false vacuum should become more stable. This is, in fact, what 
happens.... As expected, the tunnelling probability diminishes for very large values of At 
and fi, and nit '^'^^ ^™i ■ ^^ summarize, if the global CCB minimum is nearly degenerate 
with the local SML minimum (thin-wall limit), then the tunnelling probability is extremely 
small. As the trilinear couplings increase, the false vacuum decay rate increases because 
the escape point of the bounce moves out of the flat vicinity of the global minimum into the 
region in which the gradient of the potential is significant. However, a further increase 
in the size of the trilinear couplings, as well as the consequent increase in 
the squark mass terms, makes the barrier thicker and pushes the escape point 
away from the [metastable] minimum. This eventually causes a decrease in 
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Figure 1: ^''Tunneling probability for unphysically large values of At and fi. As the CCB 
minimum moves farther away, it becomes "less dangerous". As before, the stars mark the 
points with S > 400, while the boxes depict those with S < 400.". The points marked with 
stars correspond to MSSM "standard/reahstic vacua" that are long hved on the scale of 
the age of the universe. From ^54j with permission. 



the tunneling rate. In accordance with one's intuition, the low- energy physics 
is unaffected by the physics at the very high energy scales. " These effects are 
clearly visible in the Fig. [T]and indicate that \fi\, \Aq\ > 10 TeV should be utterly safe as 
regards CCB/UFB issues. 

The above arguments justify our retention of the parameter sets which we have found 
and which prima facie suffer from CCB problems. Our very large values of tan/3, |/i|, \Ao\ 
and B militate for metastability on cosmic time scales. However for the dubious we note two 
further points : we have not used the full freedom in to choose trilinears at Mx, contenting 
ourselves with a single Aq{Mx) parameter. The requirement of viable metastability could 
itself be incorporated as a search criteria in future versions of our search routines. In 
sum it would be premature to dismiss the novel scenario emerging from an apparently 
well motivated fundamental theory on the basis of technical grounds of CCB instability 
whose very efficacy in choosing vacua has long been in doubt. We note that recently 
arguments have been advanced that cast doubt on even the very possibility of tunnelling 
from a Poincare invariant vacuum to a lower energy state[SS]. Should it be the case that 
a phenomenology package assumes some version of CCB constraints without calculating 



30 



Particle =^ 


LSP 


Winos 


Gluino 


Higgs 


Sfermions 


Case J| 


(B) 


(lyi.O) 


9 


h^ 


IF 


/-I 


0.125 


0.328 


0.570 


0.111 


1.1(^),0.16(^) 


1-2 


0.105 


0.354 


0.269 


0.113 


1.89(1), 0.3(^),0.82(rf), 0.83(1) 


1-3 


0.147 


0.406 


0.599 


0.115 


1.63(e), 0.39(^) 


I -A 


0.104 


0.302 


0.351 


0.128 


0.12(^),0.89(^),1.7(J), 1.7(1) 


11 -1 


0.035 


0.113 


0.83 


0.130 


1.05(e), 1.14(^), 0.32(rf), 0.32(s), 2.49(m), 2.49(c), 2.49(Qi,2) 


II -3 


.044 


.144 


0.09 


.124 


2.15(^), 2.07(^), 0.46(rf), 0.46(1), 0.20(^), 0.20(g), 2.45(Qi,2) 


II -A 


0.035 


0.110 


0.082 


0.127 


1.67(g), 1.68(^), 0.364(d), 0.36(1), 1.36(^), 1.36(g), 2.11(gi,2) 


III -I 


0.100 


0.271 


0.429 


0.117 


0.105(e), 0.91(^), 2.29(J), 2.29(1) 


III -2 


0.99 


0.342 


0.232 


0.113 


1.87(g), 0.34(^), 0.72(J), 0.71(1) 


III -3 


0.98 


0.273 


0.387 


0.118 


1.06(g), 0.14(^), 2.15(J), 2.15(1) 


III -A 


0.94 


0.284 


0.269 


0.128 


.09(g), 0.11(^),1.17( J), 1.17(1) 



Table 20: Table of nominally discoverable particles. Mass values (in TeV) below 2.5 TeV, rounded 
off to two decimal places, calculated at tree level(except the Higgs which include one loop correc- 
tions) using two loop RGE equations including generation mixing. The principal component of 
the mass eigenstate is indicated in brackets after the mass value. The eigenstates are quite pure. 

metastability then it is the package that must be improved not necessarily our proposed 
parameters that need to be considered as discredited. 

The mass spectra obey a normal s-hierarchy (third sgeneration heavier than degenerate 
first two generations) coupled with a mini split supersymmetry (m? ^ Mi) with pure Bino 
LSP, Wino(W^±) light charginos, and next to lightest neutralino(lV3), and Higgsino heavy 
neutralinos and chargino. What is remarkable and interesting from the point of view of 
the Dark Matter Cosmology is that the quasi inert Bino LSP, which could serve as an ideal 
form of Cold Dark Matter, is generally accompanied by a light Right sfermion of the first 
or second generation (often a smuon). This is in sharp contrast to most Susy GUT spectra 
which predict the stau, stop and sbottom as the lightest sfermions because their masses are 
driven to lower values by the effect of the their large yukawas. Here however the additional 
presence of large negative Higgs mass terms drives the the third sgeneration to large masses 
(see Appendix B). Thus our model is marked out from other GUT models by a distinctive 
low energy spectrum that puts it in a different and novel universality class of models. The 
LHC at 14 TeV will provide about 1-2 TeV per colliding parton and so, as a rule of thumb, 
particles lighter than 2.5 TeV and with reasonably large couplings to SM particles may be 
detectable. Besides the light Bino which is very weakly coupled and the light Higgs /i°, the 
other MSSM gauginos Wa-, g, and from the sfermions some from among fi, e, d, s, u, c are 
lighter than 2 TeV and should be discoverable at LHC. In Table 20 we give the "detectable 
spectra" for each case. 

The cosmology of the Bino LSP Dark matter would be determined by the co-annihilation 
through the sfermion channels (the pseudo-scalar Higgs is unavailable for the purpose being 
very heavy). The combination of a Bino LSP and very light sfermion is thus ideal for Bino- 
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LSP WIMP DM and it is interesting to note that this feature emerged naturally from 
completely independent considerations. Note that-as is often the case - the smuon is the 
lightest sfermion then such annihilation could lead to a excess of charged leptons relative 
to nucleons from DM decay. Such signals have indeed been reported recently by a number 
of experiments [601 ED [Ml [63]. We will return to detailed examination of the rates of the 
co-annihilation and charged particle production therefrom elsewhere. 

In our fits typical gluino masses lie in the 200 — 600 GeV range. It is interesting that 
the ratio of gaugino masses can diverge significantly from the Mi : M2 : M3 :: 1 : 2 : 7 
ratio dictated by the 1-loop RG invariance of Mi/ gf if one begins from universal gaugino 
masses at Mx '■ as is simple and plausible in a GUT context. However we find that 
the gluino can be quite light ( as in Fit 1-2 where it is lighter than even the Wino). 
Further phenomenological analysis to revise the gluino mass bounds in the special soft 
Susy parameter region of high tan/3, and multi TeV n, B, Aq, which has emerged from this 
analysis as a distinct, novel possibly viable region of soft Susy parameter space, is required 
before the viability of such fits can be decided. This unusual feature like many others in 
this scenario, is also due to the effects of the large values of Ao that we have been required 
to consider by the down and strange quark fitting requirements of the NMSGUT. The large 
values of M"^ ^, \Aq\ required by the fermion fitting lead to correspondingly large values 
of required |/i(M^)|, _B(M^) (and therefore also of the additional pseudoscalar mass Ma 
which emerges much to heavy to play any role in Dark Matter cosmology) through the 
EWRSB conditions that tie them together. The large (negative) values of Mjj g have a 
dramatic consequence : the one loop RG equations for the evolution of sfermion masses 
contain terms proportional to Mh^M'Jj times the yukawa couplings squared YIYj. These 
terms dominate the RGBs for the third generation sfermions and drive their masses far 
above the masses of the first two s-generations independently of the value of Aq. These 
important features of the RGE are shown graphically in Appendix B. 

Baryon decay via d = 5 operators is, as usual [25| 153], dominated by the chargino medi- 
ated channels. The heavy sfermions help with suppressing B-decay. The dominant channels 
are Baryon — )■ Meson + neutrino. We emphasize that the flavour violation required by 
d = 5 B violation is supplied entirely by the rundown values of the (off diagonal) Super- 
CKM values determined by the fitting of the fermion yukawas at Mx by the SO(IO) light 
fermion yukawa formulae [9l [121 IBl 13 [H]. Using the formulae given in [14j and adapting 
the formalism of [251 [53], the proton decay decay rates in the six dominant channels for 
the 11 Fits we present in this paper are given in Table [511 These lifetimes are enhanced 
by up to 8 orders of magnitude relative to those calculated for generic fits where [ITj no 
attempt was made to suppress the coefficient of the Baryon violating d = 5 operators or 
to take account of threshold corrections to the yukawa couplings. If one tries to trace how 
this has been accomplished one sees that the minimum value of the masses of standard 
B violating triplets [3,l,±2/3](t © t) (which anyway dominate -being lighter- the novel 
[3, 3, ±2/3] {P(BP), [3, 1, ±8/3] {K © K) triplets) is raised by some two orders of magnitude 
and in addition the SO(IO) yukawa couplings are also reduced significantly relative to those 
that were necessary to reproduce the MSSM fermion yukawas in the theory without thresh- 
old corrections. A point that requires further investigation is the effect of wave function 
renormalizations on the Baryon decay operators. As matters stand we have not applied 
such corrections since the external lines are all light fermion lines and the corrections on 
these lines are typically smaller than the systematic errors in the B- violation calcula- 
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Case 


T,{M+U) 


V{p —)■ TT + Z/) 


BR{p -^ TT+Z/e,/,,^) 


T{p -^ K+u) 


BR{p -^ i^+Z/e,M,r) 


/-I 


2.4 X lO^'^ 


6.2 X 10-38 


{3.6 X 10-^0.082,0.918} 


3.5 X 10-37 


{2.8 X 10-5,0.119,0.881} 


1-2 


6.7 X 10^^ 


2.5 X 10-36 


{3.8 X 10-5,0.06,0.94} 


1.3 X 10-35 


{1.1 X 10-4,0.111,0.889} 


1-3 


5.7 X 10^6 


2.4 X 10-38 


{3.2 X 10-^ 0.100, 0.900} 


1.5 X 10-37 


{2.3 X 10-5,0.139,0.861} 


I -A 


5.7 X 1034 


1.7 X 10-36 


{7.4 X 10-5,0.052,0.948} 


1.6 X 10-35 


{9.1 X 10-5,0.046,0.954} 


11 -I 


1.5 X 10^^ 


9.7 X 10-38 


{1.8 X 10-6,0.114,0.886} 


5.9 X 10-37 


{3.6 X 10-5,0.170,0.830} 


II -3 


1.7 X 10^^ 


7.4 X 10-38 


{1.9 X 10-6,0.153,0.847} 


5.1 X 10-37 


{2.7 X 10-5,0.205,0.795} 


II -A 


6.2 X 10^3 


1.6 X 10-35 


{5.8 X 10-5,0.071,0.929} 


1.5 X 10-34 


{8.0 X 10-5,0.069,0.931} 


III -I 


2.3 X lO^'^ 


6.5 X 10-38 


{5.7 X 10-^0.088,0.912} 


3.7 X 10-37 


{2.5 X 10-5,0.128,0.872} 


III -2 


5.0 X 10^4 


3.3 X 10-36 


{3.3 X 10-5,0.050,0.950} 


1.7 X 10-35 


{9.2 X 10-5,0.093,0.907} 


III -3 


2.2 X lO^'' 


6.2 X 10-38 


{5.9 X 10-^0.103,0.897} 


4.0 X 10-37 


{2.0 X 10-5,0.141,0.859} 


III -A 


6.2 X 10^3 


1.6 X 10-35 


{5.8 X 10-5,0.071,0.929} 


1.5 X 10-34 


{8.0 X 10-5,0.069,0.931} 



Table 21: of d = 5 operator mediated proton lifetimes Tp(yrs), decay rates T[yr 
ratios in the dominant Meson"*" + u channels. 



and Branching 



tion(see tables of type X-b). It is another matter that dimension six operators containing 
external light Higgs vevs may be enhanced by wave function renormalization. However, 
the additional dimensional suppression may well keep these operators subdominant. 

We see that we have been able to suppress the B decay rates to lie comfortably within 
the current limits. Thus the search criteria may even be loosened without conflict with 
experiment. We note that our programs can already calculate the rates in other channels 
driven by Gluino, Neutralino, Higgsino etc exchange. However we defer a presentation 
of the results for the subdominant channels till the various corrections and improvements 
still needed (see below) needed have been implemented. Our aim was to show that the 
NMSGUT is quite compatible with the stability of the proton to the degree it has been 
tested, and even beyond. Firm predictions will ensue only once the susy spectrum is 
anchored in reality by a discovery of a supersymmetric particle. 

The very heavy third sgeneration masses indicate that the rate r(6 — )■ 57), is likely to 
be acceptable and uniform among the flts. The Susy contribution to muon (g-2) Aa^ = 
A{g — 2)^/2 may vary considerably since the mass of the smuon in the loop within which 
the photon couples is quite variable and generally quite low compared to other sfermions. 
Finally change in the p parameter Ap could also in principle be appreciable due to the 6-8 
light superparticle present in most cases. We plugged our susy spectra into the ( tree level) 
Spheno[56j routines to obtain the contributions shown in Table] 



The 6 — )■ S7 branching ratio values are right in the centre of the region (3— 4x 10-4)±15% 
determined by measurements at CLEO, BaBar and Belie f^ HSl SSI ED]- The current 
difference between experiment and theory for the muon magnetic moment anomaly is is 
Aa^ = 255(63) (49) x 10-^^|5T|. The results in Table [22] are thus certainly in the right ball 
park and we may well begin to use the value of Aa^ to discriminate between different models 
provided one is confldent that all instabilities in the parameter determination process have 
been controlled by adequate attention to loop and threshold effects. At the moment however 
we simply note that there is no gross conflict. The predicted change in the p parameter is 
so small as to be insigniflcant compared with the experimental uncertainties ~ .001 [51] . 
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Case 


B.Rih -^ S7) 


Aa^ 


Ap 


/-I 


3.294 X 10"'* 


5.796 X 10-9 


5.985 X 10-6 


1-2 


3.293 X 10"4 


5.471 X 10-9 


2.397 X 10-5 


1-3 


3.294 X 10-4 


2.300 X 10-9 


2.825 X 10-6 


1-4 


3.293 X 10-4 


7.238 X 10-9 


6.064 X 10-^ 


II -I 


3.290 X 10-4 


1.360 X 10-1° 


2.503 X 10-6 


II -3 


3.287 X 10-4 


1.035 X 10-1° 


3.385 X 10-6 


II -A 


3.278 X 10-4 


1.043 X 10-1° 


3.612 X 10-6 


III -I 


3.293 X 10-4 


8.058 X 10-9 


3.718 X 10-6 


III -2 


3.293 X 10-4 


6.824 X 10-9 


2.105 X 10-5 


III -3 


3.295 X 10-4 


8.689 X 10-9 


3.743 X 10-6 


III -A 


3.294 X 10-4 


7.452 X 10-9 


5.989 X 10-^ 



Table 22: Table Low energy constraints from the limits on the branching ratio for b — t- 57, Aa^ 
and Ap. 



The unification scale tends to be raised above M^ in the NMSGUT i.e. Ax > . 
This is especially true once we demand that rf = 5 operators mediating proton decay be 
suppressed. In fact of the Fits we have exhibited here and in [HI HTj the values of Ax 
we encounter are —0.29, 1.82 for the solutions without GUT scale threshold corrections 
(to fermion Yukawas) while with threshold corrections one gets {1.16,2.82,1.28,0.46} for 
Cases I-l to 1-4, {1.15,1.36,0.39} for Cases 11-1,3,4. Ac, and {1.22,2.82,1.21,0.43} for 
Cases III-l to 1-4. Thus we see that the unification scale-defined as the mass of the B- 
violating gauginos of type X[3,2, ±|] is typically raised by one oorder of magintude or 
more. On the other hand the correction to the inverse value of the fine structure constant 
(Ac) at the unification scale varies over a wide range from -20.0 to 8.1 so that the value 
of the unification coupling may as well be weak as not. However it remains true that 
above the new unification scale once we begin to use the SO(IO) RGE beta functions the 
gauge coupling will still explode |57t 158] over an energy scale range of only about 5-10. 
Smaller ac can only postpone this a little. An ideal scenario is then that the theory is 
still weakly coupled at the threshold corrected unification scale Mx > lO^^'^ GeV but that 
thereafter the Susy GUT becomes strongly coupled simultaneously with gravity. In that 
case the Planck scale may be identified as a physical cutoff for the Susy NMSGUT where it 
condenses as strongly coupled Supersymmetric gauge theory described by an appropriate 
sigma model. We envisage |58] the possibility that gravity arises dynamically as an induced 
effect of the quantum fluctuations of the Susy GUT calculated in a coordinate independent 
framework. This may be realized as a path integral over a background metric that begins to 
propagate only at low energies leading to the near canonical N=l Supergravity perturbative 
NMSGUT as the effective theory below Mpianck that we assume in our work. 

Besides the high(GUT)and low (Susy/Electroweak) scale thresholds SO(IO) theories 
are also typically subject to threshold corrections and RG flows associated with the cou- 
plings of the righthanded neutrinos present in the theory. In fact the NMSGUT scenario 
makes essential use of intermediate scale heavy neutrinos(10* — lO*'^ GeV). Such neutrinos 
are however not inert: particularly as far as their effects on RG evolution of the yukawa 
couplings of the light leptons above the mass threshold associated with the right handed 
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neutrino masses are concerned. The techniques for inclusion of these RG flows and thresh- 
old effects are by now standard |l5l |16] and in subsequent full analysis we will include also 
these effects. 

The NMSGUT also provides corrections to the QCD coupling at Mz that are in the 
right range(— .017 < /S.az{Mz) < —.004 to lower it as required [HI [19]. 

As a direct consequence of the dominance of the Type I seesaw mechanism, due to the 
low value of the 126 coupling, the right handed neutrinos emerge in just the range 10^ — 10^^ 
GeV required to implement leptogenesisjSS]. The influence of these threshold on the RG 
evolution has not yet been factored in by us yet and may have important implications for 
the Yukawa uniflcation. This is straightforward to implement and of high priority for the 
next round of improvement of the calculation. 

The large values of the crucial parameters |/i|, |y4o|,iVl4, |M|^jj| ~ lO'^TeV (where Ma 
is the mass of the pseudo scalar Higgs remnant) play a crucial role in structuring the low 
energy phenomenology. Due to the large value of |/i| the LSP is essentially a pure Bino 
and the lighter chargino is a Wino(VF±) while the heavier one H^ is a pure Higgsino, with 
mass set by the large fi parameter ( which also sets the scale for the two heavy neutralinos) 
while the next to lightest neutralino is also an essentially pure WinoCWs) . Since the required 
threshold corrections depend on ratios of scalar masses to gaugino masses there is actually 
a preference for light gaugino masses to enhance the ratios with the heavy masses. Running 
counter to this is only the constraint that the experimental lower mass limit on charged 
gauginos and exotic charged and coloured scalars generally is around 100 GeV. Thus we 
imposed a floor of 110 GeV for all such exotics. It is this chargino limit and the linked 
behaviour of the Bino and Wino masses that prevents the Bino mass from running to very 
low values( in our example flts the lowest Bino mass is 35 GeV). Also due to the link 
between the gaugino masses have We have light gluinos below 500 GeV for the cases I-II 
where we allowed LSP/Binos in the range 5 — 150 GeV. 

The parameters /i, mA at Mz are determined in terms of the run down Higgs mass 
parameters M^ ^ by the Electroweak symmetry breaking conditions which we implemented 
at the one-loop level by including Higgs tadpoles calculated using a subroutine from [5B] 
corresponding to the formulae given by [H]. Ma sets the scale of the mass of the scalar 
Higgs apart from h^ i.e H^, if ^, A all have masses close to Ma- The light scalar Higgs 
h^ typically has a mass in the range of 110 — 125 GeV. Once the values of /i, rriA ( or 
B = m\Sin2l3/2 ) at Mz are known we run them back up to Mx since they do not interfere 
with the running of rest of parameters due to the modular structure of the RGE. Thus we 
give the 7 parameter set (mi ,mo, ^o, ^^jj, 1/^1, -B) which together with the NMSGUT 
parameters completely specifles the theory at all scales and yields a distinctive scenario for 
the low energy supersymmetric accelerator phenomenology as well as for LSP Dark Matter 
cosmology. 

The large but negative values of the Higgs mass-squared parameters found by our search 
programs have a dramatic consequence that sets the type of Susy spectra obtained in 
the NMSGUT in marked contradiction with the generally accepted and used patterns 
of Sfermion masses derived from GUT boundary conditions at high scales. As shown 
in Appendix B these large negative mass-squared parameters drive the third generation 
sfermions to he much heavier than the first two generations even though we assume a 
common mass squared for the sfermions in all three 16 plets at M^. Thus third generation 
sfermion generally lie in the range 5-50 TeV and are effectively decoupled from electroweak 
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scale physics. However the NMSGUT susy spectra are not of the spht supersymmetry type 
since the sfermions of the first two generations typically populate the mass band between 
the light gauginos and the superheavy third sgeneration. Most of the sparticles emerge 
above the direct discovery limits but a couple of right chiral sfermions emerge ( see Table 
[20!) in the discoverable set; of which squarks are generally relatively heavy but right sleptons 
and specially the smuon can descend even to the Electroweak breaking scale. It remains to 
be seen if leptonic flavour violation constraints in will actually permit such light sleptons. 
However the point to emphasize at this stage is that searches based on the assumption of 
a light LSP inevitably lead to a characteristic Susy spectrum : 

M^ < {M^,Mg,Mj^J « Mp^^^ « Mj^^^^ « |/i|, |Ao|,Ma, |M|,_^| (14) 

It may be necessary to further constrain the search (and include 1-loop corrections to 
sfermion masses) to obtain consistency with flavour violation processes involving the first 
two generations (as can be seen from the Table [22] the generally stringent limits due to 
6 — )■ 57 (specially at large tan/3 ) when the third sgeneration is lightest are ineffective in 
our case). Thus whereas it may be premature to point to any one sfermion mass as a 
prediction of the NMSGUT that can be verified at the LHC it seems safe to assert that 
that the NMSGUT does predict a light Bino and Ghargino and more distinctively that the 
first sfermions discovered, at LHG or later, must belong to the first two generations with 
some preference being exhibited for the smuon and or selectron (provided those cases are 
consistent with fiavour violation constraints). 

7 Conclusions and Outlook 

This paper is the third of a series [HI [T7] developing the NMSO(10)GUT as a possibly 
viable and complete theory of particle physics. The emphasis is to develop the theory 
on the same lines and as explicitly as the the Standard Model. We have shown that the 
theory is sufficiently simple as to allow explicit calculation of the spontaneous symmetry 
breaking, mass spectra and eigenstates and allows a computation of the RG fiow in terms 
of the fundamental GUT parameters to the point where one can attempt to actually fit 
the low energy data, i.e the SM parameters together with the neutrino mixing data, in 
its entirety. In carrying out this project the model meets a major obstacle in its inability 
to fit the unmodified down and strange quark yukawas in the MSSM (renormalized up to 
Mx), which it overcomes by staking its viability on the operation of large tan/3 driven 
threshold corrections -providentially known to be operative and important in this context- 
which lower these yukawas from their SM values by a factor of about 5. As a result the 
pattern- although unfortunately not the scale- of the Susy breaking parameters tends to 
become fixed by the need to preserve, or indeed somewhat raise, the b quark yukawa while 
lowering the d, s quark yukawas when crossing the SM-MSSM threshold(s) . Thus the major 
unknown for the model, indeed for the whole field of Susy, remains the mass of at least 
one of the light susy particles. Due to the tight interrelationship of the susy spectra which 
have been generated from just 5 soft susy parameters at Mx and subjected to multiple 
stringent demands such as EW symmetry breaking, yukawa modification, Baryon decay 
suppression, LSP consistency and so on, it seems very plausible that input of any one of 
the susy particle masses would generate a prediction of all the susy masses in the context 
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of a fit tliat was pinned to yield that particular mass value along with satisfying the various 
requirements we have mentioned. In the absence of Susy discovery data and the difficulty 
of extending the mass ranges that can actually be claimed to have been excluded without 
egregious assumptions, we can still call upon the the model to stand by it's claim to be a 
viable theory of Dark matter and yield its stable Bino LSP in the 5 — 150 GeV mass range 
preferred by Cold Dark matter WIMP scenarios. This requirement is sufficient to pin the 
prediction for the Susy spectrum to a specific and completely novel and distinctive type, 
if not yet to specific masses for specific particles. The model thus predicts light gauginos 
and discoverable first or second generation right chiral sfermions below 2 TeV and invisible 
third generation sfermions. Thus we arrive at the attractively dangerous conclusion that : 
// as all other GUT scenarios derive, the third sgeneration emerges lightest the NMSGUT 
may he taken to he falsified. 

This remarkable model has thus added yet more feathers to the already long list of 
its attractions. Besides providing a natural and minimal context for the supersymmetric 
seesaw mechanism and the implementation of R-parity as a part of the gauge structure of 
unification, the theory successfully accounts for the entire available fermion mass data in 
terms of its own parameters in way consistent with all known phenomenology. It also yields 
insight into a necessary structure of it's susy breaking parameters and yields a viable and 
completely natural candidate for the LSP combined with a surprising and novel candidate 
for the scalar NLSP which can lead to an effective DM scenario. Furthermore the theory 
naturally pushes the unification scale towards the Planck scale and allows suppression of 
the dangerous d = 5 B violation operators. The confiation of the Planck scale and the 
unification scale goes a long way towards alleviating a perennial problem of renormalizable 
SO(IO) models [SZIEB] namely divergent couplings in the ultraviolet. The unification scale 
or rather the Landau pole above it becomes a physical cutoff beyond which the theory 
enters a strongly coupled phase together with gravity. 

Till a susy breaking soft mass is pinned by experiment the development of the NMSGUT 
will continue by facing up to technical challenges that we have postponed in the first phase 
of the definition of the model as reported in a series of papers( P, [T21 El [131 El [13 [HI [12] and 
the papers of the current triplet: [HI [T7| and this paper). We conclude by itemizing the 
important issues which may materially alter the numbers obtained by our fitting program 
so far. 



We have, following |44j. chosen Mz as the scale at which we match the SM to the 
MSSM. On the one hand this is well motivated since, as we have seen, the weak 
gauginos, i.e the Bino and Wino tend to emerge as light as we allow them : which 
is as light as experiment permits i.e ~ Mz for the Chargino but much lower for the 
Bino( the connection between the two masses and the limit M^r± however does not let 
the Bino mass descend below about 25 GeV). On the other hand we cannot deny that 
the model has its own little hierarchy problem with the /i, Aq, Mjj parameters all lying 
in the range of 10-100 TeV. Clearly the threshold effects due to the large differences 
of the Susy particle spectra from Mz will cause appreciable threshold corrections to 
the naive RG running which takes a single undifferentiated Susy breaking scale and 
identifies it arbitrarily with Mz- Some development of the techniques for incorporat- 
ing multiple thresholds between the SM and MSSM has already taken place[73] and 
should prove useful. 
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We have calculated only tree tree masses for the susy particles in the theory (the 
Electro weak symmetry breaking and thus the Higgs masses were however calculated 
using 1-loop effective potential [411 EG]- This is not a very good approximation in the 
MSSM with small Aq and it may be worse in scenarios with large Aq parameters. 
The large Aq parameters lead to large trilinear couplings {A = AoYf) only for the 
third sgeneration since the Yukawas 1/ for the other generations are so small and may 
significantly modify the third generation sfermion masses. The complete formulae for 
calculating these modifications-at least in the flavour diagonal case - have long been 
known [H]. They will be incorporated in the next version of our search codes. 

We have not yet incorporated the three thresholds associated with the heavy right 
handed neutrino masses in the theory. Since one progressively introduces the neu- 
trino Dirac couplings as one passes these thresholds going higher in energy, significant 
effects on the yukawa unification may be anticipated. These thresholds will also be 
significant when calculating the amount of lepton flavour violation introduced when 
integrating down from the high scale (driven by the non diagonality of the SO(IO) 
yukawas required to account for the observed quark flavour mixing) . It is for this rea- 
son, and because the mass insertion formalism is ill adapted to securely evaluate novel 
scenarios, that we have not generated tables of Lepton flavour violating mass inser- 
tions for comparison with the existing analyses on lepton flavour constraints [75]. We 
note however that we have calculated some of the common Susy sensitive quantities 
and found them to be compatible with existing limits. Nevertheless the incorporation 
of constraints on our searches based upon electric dipole moments and the strength 
of flavour violation operators involving the sfermions of the first two sgenerations are 
a priority issue for our program. 

It will be clear to the reader familiar with the nuances of the Susy Flavour problem 
that in terms of a Bottom-up approach our results suggest that a radical extension 
of the consistent and viable susy parameter space at large tan/3, Aq, \mu\, |M^| and 
with M? » Mr may be possible and the characteristics of such an extension as 
derived from the NMSO(10)GUT approach could be diametrically opposed to those 
from all previous GUTs. However that reader will also realize that the inversion 
of the standard inverted hierarchy of sfermion masses to a normal hierarchy may 
have drastic implications for the consistency of the theory with flavour violation 
constraints and EDM constraints for the first two generations (the heavy sthird is 
likely to pass such constraints -as already seen in the case of the 6 — )■ 57 constraint 
: which is commonly a stringent one when the third sgeneration is light but is here 
totally insensitive see Table |22|) . Prominent among these well known constraints are 
those on the electric dipole moments of the electron, muon, neutron etc, the severely 
constrained Branching rations of the AF = 1 decays of Mesons such as Bd — )■ fJ^fi, 
and the limits on the supersymmetric contributions to AF = 2 quantities like e^, AK 
etc. These contributions can now be calculated using codes[7H] that do not make any 
simplifying approximations as done in the mass insertion method [7^] which allow 
only order of magnitude estimates of individual terms but cannot be used to evaluate 
total contributions including the effect of cancellations, in unconventional scenarios 
like ours. It may well be that the parameter sets we have presented will fail these tests 
when they put to them. However we emphasize that, from the point of view of the 
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this series of papers which seeks to pin down the viable points of the 45 dimensional 
parameter space SO(IO) NMSGUT, such a challenge would be no different than the 
one that was posed by the onus to show compatibility with B violation rates: which 
the theory overcame by novel use of expedients available to it in a way that pointed 
out fresh approaches to hoary questions. In the same way it may be that when 
the requirement of the consistency of the total value of such parameters predicted 
by the NMSGUT is taken into consideration then again new viable parameter sets 
may emerge. In view of the length of the present paper as well as the considerable 
further effort required to answer definitively to these important challenges the flavour 
violation constraints will be taken up in the sequels. We urge the tolerant reader not 
to prejudge this vital issue but join us in reflecting on the behaviour of the novel type 
of Susy parameter sets suggested by us. 

In this calculation we found parameter sets leading to quasi-stable Baryons by a shot- 
gun carrying brute search. However the characterization of the possible cancellations 
and the least constraining versions compatible with experimental limits remains to be 
done. The search for fits with suppressed Baryon Number violation was carried out 
by simply limiting the size of the maximal element of the LLLL and RRRR operator 
coefficients with no heed paid to the detailed effects of that coefficient on the decay 
rate in specific channels. A more sophisticated (but hugely more computer intensive) 
way of doing this would be to calculate the Baryon violation rates in each channel at 
every iteration and limit the total lifetime. In principle one could hope to implement 
this given enough (super) computing power. 

Due to the large amount of running time required to find an acceptable solution, 
we have only scratched the surface of the enormous parameter space and can by no 
means pronounce on the general structure of the solution space. It will take long runs 
on a super computing cluster to develop a statistical picture of where the solutions 
tend to lie. We are now preparing for both the improvements mentioned above and 
the harnessing of a cluster for the task. 

We have run down only the diagonalized yukawa couplings from M^ to Mz along 
with the susy breaking parameters which are optimized to fit the eigenvalues of the 
SM yukawa couplings to the run down diagonal NMSGUT matter yukawas. A com- 
plete treatment would run down the full set of coupling matrices obtained from the 
GUT and apply the large tan/3 driven Susy threshold corrections to the off diago- 
nal couplings before matching the two sets : or at least the "physical" parameters 
(eigenvalues, mixing angles and phases) coded in the two pairs. The formalism for 
applying off diagonal corrections is still somewhat murky so some theoretical progress 
towards setting out clear algorithms for including all off diagonal 1-loop effects and 
renormalizations is called for. 

We argued that the CCB and UFB minimae(that certainly exist at the large values 
of |/i|, \Aq\ central to the needs of the NMSGUT) need not destabilize the metastable 
standard vacua we have found and cited previous investigations of this issue |54] which 
are specially encouraging in this regard. The fact remains, however, that these decay 
rates must actually be estimated for each set of parameters found. On the other hand 
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once the requisite subroutines for calculating the vacuum tunneling rate are in hand 
one can add them to the search routines to filter the parameter sets. 

Threshold corrections play a central role in our calculation and the large wavefunction 
dressing of the Higgs doublet lines that we find at one loop demand an investigation 
of the two loop effects to determine whether they are yet larger still. If so our 
model, and by implication realistic GUTs generally, may survive only after the Higgs 
wavefunction dressing has somehow been re-summed to all orders. However we note 
that this growth of wave function dressing leads to reduction in the actual size of 
the SO(IO) yukawas actually needed to fit the low energy fermion data by a factor 
of 10-100. Thus the contributions to the fermion lines are effectively degenerate and 
it is only the large number of heavy fields running in dressing of the light Higgs 
boson lines that leads to such a large wave function renormalization of the Higgs 
fields. It is possible that in case of further growth at the multiloop level the theory is 
being driven to a quantum fixed point dominated by the gauge coupling and therefore 
re-summable using the the exact beta functions available for supersymmetric gauge 
theories [77]. In any event we consider that our calculations show that GUTs aiming 
to be realistic have been "up-ended": in the sense that any Grand unified model 
with pretensions to realism must define itself explicitly enough to permit calculation 
of threshold corrections using calculated superheavy field spectra or else risk being 
under suspicion of being an qualitative scenario that may be destroyed as soon as 
quantum effects are included. 

As Karl Popper emphasized so insightfully, the virtue of a comprehensive scientific 
theory that accounts for all known data is that it accepts the challenge and charm 
of living dangerously and in constant confrontation with experimental data that may 
falsify it. It must face every new measurement that it claims to be able to account 
for with its fate hanging in the balance. This indeed is what gives properly scientific 
models their peculiar power and utility: that speculations living at a safe distance 
from the cutting edge of Occam's razor rarely possess. The (N)MSO(10)GUT has 
braved multiple iterations of challenges over the three decades since it was proposed [U 
[2] and, so far, has emerged stronger from every challenge, defining the possible in its 
realm ever more clearly and distinctly. If it dies it will not have lived in vain. 
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Appendix A: Gut Scale threshold corrections to matter 
yukawas 

We give below our results for the threshold corrections to the Yukawa couplings of the 
matter fields due to heavy fields running in a self energy loop on a line leading into the 
Yukawa vertex(here we assume the Higgs line is predominantly the 10-plet derived doublet). 
Let us denote by U'^iV^) the matrices that diagonalize the mass terms for fields of the 
alphabetized label type A : 

$'^M$ = W'^Moiag^' <= $ = W^W ; $ = 1/*$' (Al) 

generically by W"^ and the corresponding index ranging over the multiplicity of that field 
type also by the corresponding lower case roman letter (a). For example for the 6- fold set 

[3, 2, — |](_Ei, E2, £"3, £'4, E^, Eq)Q)[S, 2, |](i?i, E2, E3, E4, £5, Eq) the generic rotation matrix 
is denoted W^{ i.e U^orV^ ) carrying indices e, e' = 1,2,3,4,5,6. Then it is useful to 
define 

a'=dim{W^) 

J-i(W-^a) = Y. \W^,a'\'Fn{m,,,Q) 

a'=l 

h'=2 
a'=dim{W^) 

J-r(l^^,a,m(")) = Yl \W^J'F,2{m^-\ma',Q) 

a'=l 
h'=dim{W") 

h'=2 
h'=dim(W^) a'=dim(W^) 

b'=l a'=l 

h'=6 a'=dim{W^) 

h'=2 a'=\ 

a'=dim{W^) 

+ E \Wt'\'\Ki\'^n{m^''Q) 

a'=l 

a"=dim(W^) 

Ci(W^^,a,a') = Y {Ka"rW<a"Fii{ma'',Q) (A2) 

a"=l 
a"=dim(W") 

C'.iW^'.Kh') = Y {W^,,"rWi^,H"Fn{m,.,Q) 

h"=2 
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b"=dim{W'^) a'=dim(W^) 

b"=l a' = l 

h"=6 a'=dim{W^) 

C',{W'',W^,a,h,h') = E E \W,U'{W,^^,„rW^^,„F,2{ma,m,.,Q) 

h"=2 a'=l 

a'=dim{W^) 

a'=l 

The primes on the function names and summations instruct an omission of any hght 
fields (in practice the hght Higgs [1, 2, ±1] doublets only) from the sum over the heavy fields 
of the given type. If the field is one of the unmixed types (i.e m = A, B, I, M, N, O, S, T, 
U, V, W, Y, Z) then the function carries a superscript (u) thus e.g ^l (Vp, l,Mv) arises 
from a coupling between Fi[l, 1, 2] and V[l, 2, —3]. Such functions arise in the dressing of 
the Higgs lines. The calculation is quite tedious but we applied various consistency checks 
to ensure that we had included contributions from all members of multiplets. Naturally we 
await the contributions of those patient and interested enough to check our results. 



(327r2)As = 2h*hTi{UT, 1) - 4^*^J'i(f/r, 7) - 2iV2h*gCi{UT, 1,7)- 2iV2g*hCi{UT, 7, 1) 

+ h*hJ^,{VT, 1) - 2g*9J'i{VT, 7) - 2g*-9J'i{VT, 6) 

- tiV2h*gCiiVT, 1, 7)) - tiV2g*f£iiVT, 7, 1)) - V2h*gCiiVT, 1, 6) 
+ V2g*hCi{VT,QA) -2t{-g*gCi{VT,6,7)) + 2t{rgCi{VT,7,Q)) 

+ 2Ch*hJ^[{VH, 1) - {l/3)g*gJ'[{VH, 6) + {^)h*gC[{VH, 1, 6)) 

+ {{^)g*hC[{VH, 6, 1)) - g*9J"i{VH, 5) + h*gC[{VH, 1, 5) 

- rhC[{VH, 5, 1) - {{^)-g*-gC[{VH, 5, 6)) + {{^)-g*-gC[{VH, 6, 5))) 

- r^((64/3)J-i(Vb,3) + 8J-i(Vz5,3) 

+ 8J^,{Uk, 2) + 4J-i(Vj, 5) + 8J-i(Vl, 2))) - {2g^) {25 J^^{VG, 6) 

+ 0.5J-i(Vj,4) + 0.5J-i(VV,3) + 2J-i(\/x,3)+J-i(\/e,5)) (A3) 



(327r2)A,- = {2h*hJ^^{UT,l)-4g*gMUT,7) + 2^{V2h*gC^{UT,l,7)) 

+ 2V2g*hC,iUT, 7, 1)) + h*hT^iVT, 1) - 2g*-gJ'iiVT, 7) 

- 2rgJ'iiVT, 6) + tiV2h*gCiiVT, 1, 7)) + tiV2g*hCiiVT, 7, 1)) 

- ^/2h*gCl{VT, 1, 6) + V2g*hCi{VT, 6, 1) + 2t{g*gC,{VT, 6, 7)) 

- 2tig*9CiiVT,7,Q)) + 2Ch*hT[iUH,l) - il/S)r9mUH,Q) 

- fgJ^[(UH,5)+tii^)h*gC[iUHA,6))+tii^)g*hC[iUH,6,l)) 
+ h*gC[{UH, 1, 5) - g*hC[{UH, 5, 1) - i{{^)g*gC[{UH, 5, 6)) 
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+ t{{^)g*gC[{UH, 6, 5))) - g*gm/3)J^^{Uc, 3) + 8J^^{Ve, 6) + 4J^,{Vk, 2) 

+ 8J-i(f/j,5) + 8J-i(\4,2)) - (2/)(0.225J-i(Vg,6) 

+ 0.5J-i(Vj,4) + 0.5J-i(VV,3)+J-i(Vx,3) + 2J-i(\/^,5))) (A4) 



(327r=^)Ap = {2h*hT[iVH, 1) - 2r^J-((VH,5) - Qg*gT[iVH,Q) 

+2h*gC[{VH, 1, 5) - t2V3h*gC[iVH, 1, 6) + i2V3rgC[iVH, 5, 6) 

+ (2r^C;(\/^, 1, 5) - t2V3h*gC[iVH, 1, 6) + i2V3g*gC[iVH, 5, 6))t 

+3h*hT^iVT, 1) - 6^*^J-i(\/t, 7) - QrgJ'iiVT, 6) 

+3y2r^Ci(\/T, 1, 6) - t3iV2h*gCiiVT, 1, 7)) + QtirgCiiVr, 6, 7)) 

+ (3v^r^Ci(VT, 1, 6) - t3iV2h*gCiiVT, 1,7)) + 6tig*gCiiVT, 6, 7)))t 

-r^(4J-i(f/^,4) + 24J-i(f/^,6) + 12J-i(K7,5)) 

-(2(72) (0.625 J-i(Vg, 6) + 1.5J-i(V>, 4) + 0.5T^{Vf, 3) + 3TiiVE, 5))) (A5) 



(327r2)Ag = {2h*hr,{UH, 1) - 2^*^J-((f/^, 5) - ergJ^^iUu, 6) 

+2/i*^CUf//f, 1, 5) - t2V3h*gC[{UH, 1, 6) + t2V3g*gC[{UH, 5, 6) 

+ (2r^C((f/H, 1, 5) - i2V3h*gC[{UH, 1, 6) + t2y/3g*gC[{UH, 5, 6))^ 

+3h*hTiiVT, 1) - 6^*^J-i(Vt, 7) - QrgJ'iiVT, 6) 

+3v^/i*^Ci(\/t, 1, 6) + t3{V2h*gC,{VT, 1, 7)) - 2t3{g*gC,{VT, 6, 7)) 

+ (3y2r^Ci(\/T, 1, 6) + i3iV2h*gC^iVT, 1, 7)) - 2i3{g*gCi{VT, 6, 7)))^ 

-r^(4J-i(f/^,4) + 24J-i(f/z,,3) + 12J-i(\/^,2)) 

-(2/)(025J-i(Vg, 6) + l.bTiiVj, 4) + 0.5J-i(VV, 3) + 3J-i(\/x, 3))) (A6) 



\t 



(327r2)A„ = (r^J-i(f/T,l)-2r^J-i(f/r,6)-y2/i*^Ci(f/T,l,6) 
+y/2g*hCi{UT, 6, 1) + 2/i*/iJ-i(Vt, 1) + h*hJ^[{UH, 1) 
-r^J-((f/H,5) - {l/3)g*gJ'[{UH,6) - h*gC[{UH, 1,5) 
-(2/V3)/i*^CUf/H,l,6) - it/V3)g*gC[{UH,5,6) + (-/i*^C;(f/^, 1, 5) 
-(Vv^)/^*^CUf///, 1, 6) - {i/V3)rgC[{UH, 5, 6))^ + h*hJ'[{VH, 1) 
-r^^i(^H,5) - il/3)rgJ'[iVH,Q) - h*gC[iVH, 1,5) 
-{t/y/3)h*gC[iVH, 1,6)- (2/V3)^*^CUV^h, 5, 6) + i-h*gC[iVH, 1, 5) 
-(Vv^)/^*^C;(\/h, 1,6)- {i/V3)g*gC[{VH, 5, 6))t + g*gi-{S2/3)J^,{Uc, 3) 
-(32/3) J-i(Vc, 3) - 4J-i(f/z,, 3) - 4J-i(f/,,, 6) 
-4J-i(\/p, 2) - 8J-i(V>, 2) - 6J-i(f/p, 2) 
-8J-i(f/L,2)) - (2/)(25J-i(Vg,6) + 0.5J-i(Vj,4) 
+ 1.5J-i(Vx, 3) + 1.5J^i{Ve, 5))) = (327r2)A, (A7) 
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(327r2)Ae = {3h*hJ^i{UT, 1) - 6g*gJ^i{UT, 6) + 3V2h*gCi{UT, 1, 6) 
-3V2g.hC,{UT, 6, 1) + h*hJ^[{VH, 1) - g*9J"i{VH, 5) 
-3-g*9J'[iVH,Q) + h*gC[iVH,l,5) - itV3)h*gC[iVH,l,Q) 
+ {iV3)g*9C[{VH, 5, 6) + {h*9C[iVH, 1, 5) - {tV3)h*gC[iVH, 1, 6) 
+ (2V^)^*^C;(Vh, 5, 6))t + h*hT[{UH, 1) - r^^((f/H, 5) 
-3^*^J-((f/H, 6) - h*gC[{UH, 1, 5) + {iV3)h*gC[{UH, 1, 6) 
+ (2v^)^*^C;(t/H, 5, 6) + i-h*gC[iUH, 1, 5) + iiV3)h*gC[iUH, 1, 6) 
+ (^v^)^*^C;(f/H, 5, 6))t + ^*^(-4J-i(VV, 4) - 12T,iVE, 6) 
-12J-i(Vd,3) - 18J-i(f/p,2)) - (2/)((9/40)J-i(Vg,6) 
+ 1.5J-i(K7, 4) + l.bTiiUx, 3) + 1.5J-i(f/E, 5)) = (327r2)A, (A8) 

^^^!!2t''" = f J'2iVG,UH, 2,2) + ^'T^iVG,UH, 2,3) +^^C',iVG,UH, 2,2,3) 
r 11 1 

+77(C^(Vg, Uh, 2, 2, 3))« + SifT^iVn, Uc, 1, 1) + 7'^2(V^i?, f/c, 1, 2) 
+77C2(V}?, Uc, 1,1,2)+ 77(C2(Vij, f/c, 1, 1, 2))*) + 3(7^2(^7, f/^, 2, 2) 
+ (7)^2(^7, f/^, 2, 1) + 77^2(^7, Ud, 2, 2, 1) + 77(^2(^7, Ud, 2, 2, 1))*) 
+3((7) V2(f/j, Ve, 2, 1) + (7)'^2(f/j, l^E, 2, 2) + 77C2(f/j, I^e, 2, 1, 2) 
+77(C2(f/j, l^E, 2, 1, 2))*) + 3(f J-2(f/i,, Vt, 4, 2) + 7'^2(f/i=;, Vt, 4, 3) 
+77C2(t/i=;, l^T, 4, 3, 2) + 77(C2(f/E, Vt, 4, 3, 2))*) + 3(f J-2(\/x, J/t, 2, 2) 
+7'^2(Vx, Ut, 2, 3) + 77C2(V3f , Ut, 2, 3, 2) + 77(C2(Vx, Ut, 2, 3, 2))*) 
+67'^r(VL, 1, My) + 1272 J-i(Mb, Mm) + 37' J-2(V3,, f/r, 1, 4) 
+67'^2(^E, Uj, 3, 1) + 18f Fi(My, Mw) + 9f J-2(Vx, Up, 1, 1) 
+7=^ J-r(VV, 1, Mv) + 3^^T"tiVH, 4, Mo) + 272 J-^(Vg, f/^, 4, 4) 
+187'Fi(Mb, M^y) + 97'J'2(l^P, f/^, 1, 3) + GfJ^^VL, 1, Mb) 
+3f J-2(Vt, t/s, 4, 3) + 12fFi (My, Mjv) + 3f Fi(My, Mq) 
+2f Fi(My, M^) + f J-^(f/p, Vh, 1,4) + 6f J-2(\/x, f/x, 1, 1) 
+AfT2iVR, Uc, 2, 1) + A^^T2iVR, Uc, 2, 2) + O.bfT^iVc, Uh, 3, 2) 
+0.5^^J'^{Vg, Uh, 3, 3) + 1.57=^^-2(^7, Ud, 3, 1) + l.Sf ^-2(^7, Ud, 3, 2) 
+1.57'J-2(K7, Ue, 3, 2) + 1.5f J-2(K7, Ue, 3, 1) + S^^T^iVc, 1, Mz) 
+8f J-r(Vc, 2, Mz) + 7'^2(f^F, Vh, 2, 3) + fJ^^iUE, Vh, 2, 2) 
+3^'J'2{Vt, Ue, 5, 1) + 3f J-2(\/t, Ue, 5, 2) + 3fj'^iVD, 1, M,) 
+37' J-r(VD, 2, Mj) + 127^ J-r(f/c, 2, Mq) + 12f J-r(f/c;, 1, Mq) 
+1.57'^'" (f^/7, 3, Ms) + 1.5f ^';'(f/^, 2, Ms) + 4.5f J-r(f/i,, 1, M^;) 
+4.57'J'r(f/D, 2, Mt/) + 4.5f J-i"(V^E, 1, Mu) + 4.57Vr(l"E, 2, Mu) 
+ \k\\j^!,{VG, Uh, 1, 5) + 6J-r(Vi, 2, M^) + 3J^2{Vt, Ue, 6, 3) 
+6 J-^f/L, 2, My) + 3J-2(\/x, f/T, 1, 6) + T^iUp, Vh, 4, 4) 
+:F]^{Vp, 4, My) + AT^iVp, Uc, 2, 3) + 8T^{Vc, 3, Mz) 
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+0.5J-^(\/g, Uh, 2, 6) + T^iUp, Vh, 2, 6) + 1.5T2iVE, Uj, 6, 3) 
+3TnVD, 3, Mi) + 1.5J-2(K7, Ud, 3, 3) + ST^iVr, Ue, 5, 6) 
+1.5^'" (Vh, 6, M5) + UJ^nVc, 3, Mq) + A.5T^{Vd, 3, M^/) 
+4.5J-r(f/i,, 6, Mu) + l.5T2iVT, Ue, 7, 4) + 3J-2(\/x, f/x, 2, 2) 
+1.5 J-2(Vx, f/T, 2, 7) + 3T2iVE, Uj, 4, 5) + A.bT^iVp, Ue, 2, 4) 
+4.5J-2(\/x, f/p, 2, 2)) - (2/)(+1.5J-2(\/t, f/x, 1, 3) 
+1.5J-2(\/^, Ut, 5, 1) + 0.1J-^(Vb, Vh, 6, 1) + 0.5J-^(W, Uh, 3, 1))) 



(A9) 



(327r2)A^o 



|[/,^|2 



'111 



(8(f J-2(f/R, \/c, 1, 2) + (7) V2([/,j, Vc, 1, 1) + 1iC2{Ur, Vc, I, 2, 1) 

+11{C2{Ur, Vc, 1, 2, 1))*) + f J-^(f/G, \/h, 2, 2) + (7)'-^2(f^G, ^h, 2, 3) 

+77C^(f/G, V^H, 2, 2, 3) + 77(C^(f/G, ^/f , 2, 2, 3))* + 3(f J-2(f/j, Vp, 2, 1) 

+ (7)'^2(?7j, Vz), 2, 2) + 77C2(f/j, V^d, 2, 1, 2) + 77(C2(f/j, Vp, 2, 1, 2))*) 

+3(f J-2(K7, Ue, 2, 2) + (7)'^2(Vj, Ue, 2, 1) + 77^2(^7, Ue, 2, 2, 1) 

+77(^2(^7, Ue, 2, 2, 1))*) + 3(f J-2(f/x, V^t, 2, 2) 

+ (7)'^2(?7x, Vt, 2, 3) + 77C2(f/x, Vt, 2, 2, 3) + 77(C2(f/x, ^t, 2, 2, 3))*) 

+Kf^2{VE, Ut, 4, 2) + (7'^2(^i^, Ut, 4, 3) + 77^2(^7^, Ut, 4, 3, 2) 

+77(C2(Ve, Ut, 4, 3, 2))*) + 12^^Fi{My, Mm) + 67'J-r(f/L, 1, Mb) 

+37'j-2(f/T, l^i^, 4, 3) + 67' J'2(f/x, V^x, 1, 1) + 18f Fi(Mb, M;^) 

+9f J-2(\/i^, Up, 3, 1) + 37=^Fi(My, Mo) + 2^^F,{Mv, Ma) + l^HiVp, Uh, 1, 4) 

+187^^1 (My, Mw) + fJ'iiUp, 1, My) + 3f Fi(Mv, Mo) 

+2f J-^(\/g, Kh, 4, 4) + 972 J-2(V^p, [/x, 1, 1) + 3f J-2(\/t, f/x, 4, 1) 

+6f 7-2(^7, f/i^, 1, 4) + 6f J-r(^L, 1, ^y) + 12f Fi(Mb, Mm) 

+4f J-2(\/ij, l^c, 2, 2) + 472j-2(Vp, Vb, 2, 1) + O.Sf J-^(\/g, f/^, 3, 2) 

+0.57'J'2(l^/7, f/c, 3, 3) + 1.57^2(1/0, f/j, 2, 3) + 1.67^2(^0, Uj, I, 3) 

+1.57'J'2(K7, Up, 3, 1) + 1.5f J-2(K7, Ue, 3, 2) + 87Vr(f/c, 2, M^) 

+8f J-r(f/c, l,Mz)+ l^KiVp, Uh, 2, 3) 

+f J-^(VV, f//7, 2, 2) + 37'J-2(f/T, l^s, 5, 2) + 3f J-2(f/T, V,,, 5, 1) 

+3f J-r(f/z>, 2, M,) + 37' J-r(f/D, 1, Mi) + 127^ J-«(Vc, 1, Mq) 

+12f J-«(V^c, 2, Mq) + 1.57=^r;'(\/H, 3, Ms) + l.Sf r;'(\/H, 2, Ms, Mh, 6) 

+4.5f J-i"(f/E, 2, Mu + 4.57'J'r(f/E, 1, M^/ + 4.5f J-i"(Vd, 1, Mu) 

+4.57'J-r(V^D, 2, M^) + |A;|2(J-^(Vg, V^h, 1, 5) + QJ^'^iyp, 2, My) 

+3 J-2(f/x, V^T, 1, 6) + J-^(VV, f/^, 4, 4) + J-f (?7p, 4, My) 

+6J-r(f^L, 2, Mb) + 3J-2(f/T, Vp, 6, 3) + 8J-r(f/c;, 3, Mz) 

+F'^{Vp, Uh, 2, 6) + A:F2{Vr, Vc, 2, 3) + 0.5J-^(Vg, 1^7, 3, 6) 

+1.5J-2(f/i,, K7, 6, 3) + 3J-r(f/D, 3, M,) + 1.5J-2(\/z5, f/j, 3, 3) 

+3J-2(\4;, Ut, 6, 5) + l.S^'i (f/^, 6, Ms) + 12J-r(f/c, 3, Mq) 

+4.5 J-^t/o, 3, Mc;) + A.^J^'^iVp, 6, Mt;) + 1.5 J-2(\/s, Ut, 4, 7) 
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+3 J-2(f/x, Vx, 2, 2) + 1.5J-2(f/x, Vt, 2, 7) + 3J-2(f/i,, K/, 4, 5) 

+4.5J-2(t/p, Ve, 2, 4) + 4.5J-2([/x, Vp, 2, 2)) - (2(7=^)(+1.5J-2(f/T, Vx. 1, 3) (15) 

+1.5J-2(f/E, T, 5, 1) + 0.1J-^(G, Uh, 6, 1) + 0.5J-^(?7ir, H, 3, 1))) (AlO) 



Appendix B : Discussion of RGE features 

We used the two loop Renormalization group evolution equations for the softly broken 
MSSM given in |75j to evolve randomly chosen SUGRY-NUHM parameters {wi ?, rrii/2, ^o? B, 
Mjj g} together with the /x parameter and the rest of the superpotential parameters down 
from Mx = 10^^'^^ GeV to Mz- It is notable that the universal gaugino mass M1/2 at 
Mx is found to be negative in all our fits. Furthermore the Higgs mass parameters M^ g 
are very large and negative being typically ~ — 10^ TeV^. The sfermion beta functions at 
one loop contain terms proportional to the Higgs mass parameters squared [7H, [75] . For 
example : 

/^S = (m^, + 2M^)Yt Y. + (m| + 2M|)YtY, + (16) 

which dominate the RGE for the third generation sfermions and drive their masses far 
above the those of the first two generations as one flows from Mx to Mz- This behaviour 
is only slightly modulated by the contributions of Aq and is a one-loop feature immune to 
significant modification by the two loop contributions. The presence of terms A^A added 
to the Higgs contributions has a countervailing effect on the scalar mass evolution since it 
tends to decrease the mass squared in the infrared. However the huge Higgs masses always 
prevail resulting in the third sgeneration always being heavier than the first two. This is 
an invariant feature of our spectra. Heavy third sgeneration is a distinctive feature of our 
fits and counterposes them to all (to our knowledge) previous GUT based predictions which 
have a third sgeneration lighter than the first two. 

These tendencies can be clearly seen in the plots of the RG evolution given in Fig- 
ures 2-5 which refer to the actual two loop RG evolution (Fig. 2) of M^ , and the next 
three figures refer to hypothetical cases with {Ao{Mx) = 0}(Fig.3), {M]j{Mx) = = 
M|(Mx)}(Fig.4),{M^(Mx) = = Ml{Mx) = Ao(Mx)}(Fig.5). 

Similarly we can understand why the one loop prediction (Mi : M2 : M3 as gf '■ g^ '■ g1 
which is ~ 1 : 2 : 7 at Mz-) for the ratio of gaugino masses Mj which follows from the 
1-loop RG invariance of Mi/ gf., -i = 1, 2, 3 is badly violated at two loops when Aq is large. 
Although the 1:2:7 seems set in stone by the known gauge couplings at Mz if GUT 
mandated equahty of gaugino masses at Mx is accepted, the influence of the additional 
terms J74il75j at two loop in the gaugino mass RGE : 

i^a = ^A''Ma + j^,\J:B^:^gl{M. + M,) (17) 



b=l 



+ y: C:{Tt[Y^A,]-M,Tt[YX] 



:it 
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Figure 2: Two loop RG evolution of Ml from M^ to Mz for Case I-l. Red: Ml, Blue: 
M|, Green: M^. Note the strong growth in the the third sgeneration mass at low energies. 
The same behaviour is exhibited by all sfermions. 
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Figure 3: Little effect of Aq : Hypothetical Two loop RG evolution of M"^ from M^ to Mz 
with Aq{Mx) = for Case I-l. Red: Ml, Blue: M?, Green: M|. Note the strong growth 
in the the third sgeneration mass at low energies. The same behaviour is exhibited by all 
sfermions. 
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Figure 4: Effect of large M|^^ 



Hypothetical Two loop RG evolution of M^ from M^ 
to Mz with Ml{Mx) = M]^{Mx) = for Case I-l. Red: Ml, Blue: M?, Green: M|. 
Note the strong decrease in the the third sgeneration mass at low energies while the first 
two generations are unaffected. Putting Aq = has essentially no effect except that the 
increase becomes linear. The same behaviour is exhibited by all sfermions. 
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Figure 5: Little effect of A^ : Hypothetical Two loop RG evolution of M'j from M^ to Mz 
with Mfj{Mx) = M]^{Mx) = = Ao{Mx) for Case I-l. Note the strong decrease in the 
the third sgeneration mass at low energies while the first two generations are unaffected. 
The removal of the curvilinear decrease in favour of a linear one is the only effect of putting 
y4o = in addition. The same behaviour is exhibited by all sfermions. 
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Figure 6: Hypothetical 2-loop RG evolution of gaugino masses with Aq j^ (full lines) 
and with Ao = (dashed lines) for Case II-l. Red:Mi, Blue M2, Green M3. Notice how 
the Gluino mass can even fall below the Wino mass when Ao 7^ 0. 

The terms containing the product of the Yukawa gauge couplings and the corresponding 
soft trilinear couplings that are generated from Aq, imply that in practice the ratios can 
vary widely if Aq is large and gluinos can even be lighter than winos : This is seen clearly 
from the graphs of the RG flow of the gaugino masses with and without Aq (Fig. 6) and 
the graph of the ratios of the gaugino masses with and without 74o(Fig. 7). 
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Figure 7: Hypothetical Two loop RG evolution of ratios of gaugino mass ratos with Aq y^ 
(full lines) and with Aq = (dashed lines) for Case II-l. Red : M3/M1, Blue: M2/M1 , 
Green M3/M2. In the case Aq = 0, the gaugino masses follow the standard evolution to 
the 1:2:7 ratio at low energies. 

Appendix C: Additional Tables of Parameter values 
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Parameter 


Value 


Field 


Masses 






[5[/(3),5f/(2),r] 


{Units ofW^GeV) 


Xx 


0.0232 


A 1,1,4] 


17462.07 


Xz 


0.0712 


5[6,2,5/3] 


0.4028 


/ill/10-6 


-0.4808 


C[8,2, 1] 


733.73,2706.98,4133.40 


/i22/10-^ 


0.3364 


£'[3,2,7/3] 


1514.69,9791.18,11582.78 


/i33 


0.0141 


E[3,2,l/3] 


0.91,529.34,529.34 


/11/IO-6 


-0.0989 - 0.02442 




1724.816,5358.19,7929.76 


/12/IO-6 


1.6150 + 0.0073i 


i^[l,l,2] 


202.02,202.02 


/13/IO-5 


-0.6328 - 0.0360i 




1129.96,8959.17 


/22/IO-5 


-1.9062 + 0.0805i 


^[1,1,0] 


0.088,0.73,3.65 


/23/10-^ 


1.0036 + 0.0037i 




3.646,128.52,128.85 


/33/10-=^ 


-0.4254 + 0.029H 


/i[l,2,l] 


1.727,930.24,1430.99 


^12/10-^ 


-2.8201 - 0.7277i 




10152.47,13862.15 


^13/10-^ 


-57.2234 -43.9416i 


J[3, 1,10/3] 


1.02 


^23/10-^ 


135.6224 + 115.2576i 


J[3, 1,4/3] 


1.224,355.07,355.07 


A/10-2 


0.0664 + 0.0432i 




642.66, 7578.30 


V 


-8.7343 - 0.3535i 


^[3,1,8/3] 


1099.55,11032.05 


P 


1.2543 - 0.6648? 


L[6, 1,2/3] 


878.68, 14142.89 


k 


-0.0708 - 0.0060i 


M[6, 1,8/3] 


14763.05 


C 


0.6792 + 1.3334i 


iV[6, 1,4/3] 


14078.87 


c 


3.5786 - 0.5914i 


0[1,3,2] 


26102.41 


m/W^GeV 


0.05 


P[3,3,2/3] 


1157.34,18368.19 


nio/lQ^^GeV 


-947.952e-*^'^3(^) 


g[8,3,0] 


0.718 


7 


1.01 


i?[8,l,0] 


0.46,1.53 


7 


-1.3309 


5[1,3,0] 


1.6833 


X 


0.8559 + 1.0678i 


t[3, 1,2/3] 


0.87,248.68,833.99,1244.35 


Ax 


2.82 




4579.86,4885.32,76873.18 


Ag 


-19.885 


f/[3,3,4/3] 


1.345 


Aa3{Mz) 


-0.011 


\/[l,2,3] 


1.197 


{M'^yW^GeV} 


0.00,12.97,1184.75 


Vr[6,3,2/3] 


10511.05 


{M^j/10-^^eV} 


0.0012,4.56,416.72 


X[3,2,5/3] 


0.439,554.214,554.214 


M^{meV) 


2.14,7.64,42.49 


F[6,2,l/3] 


0.53 


{Evals[f]}/10-^ 


0.01,49.25,4499.03 


Z[8,l,2] 


1.51 


Soft parameters 


mi = -387.435 

2 


mo = 2609.655 


Ao = -1.8398 X 10^ 


at Mx 


/i = 1.4670 X 10^ 


B = -1.5219 X 10^° 


tan/3 = 50.0000 




M| = -1.4822 X 10^° 


M2j = -1.4528 X 10^° 


Rbr = 0.8725 


Max{\LABCD\, \Rabcd ) 


7.1909 X 10-22 GeV-i 




S/l 



Table 23: I-2-a: Column 1 contains values of the NMSGUT-SUGRY-NUHM parameters at Mx 
derived from an accurate fit to all 18 fermion data and compatible with RG constraints. See 
caption to Table 2 for explanation. 
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Parameter 


Target = Oi 


Uncert. = 5i 


Achieved = Oi 


Pull = (0, - 0,)/6, 


?/./io-6 


2.106596 


0.804720 


2.107032 


0.000542 


l/c/10-' 


1.026873 


0.169434 


1.027891 


0.006009 


Vt 


0.388445 


0.015538 


0.388409 


-0.002299 


l/d/10-^ 


5.979116 


3.485824 


6.027227 


0.013802 


l/s/10-3 


1.131477 


0.534057 


1.125268 


-0.011625 


Vb 


0.544641 


0.282669 


0.545745 


0.003905 


1/e/lO-^ 


1.279312 


0.191897 


1.279764 


0.002359 


1/m/10-^ 


2.593121 


0.388968 


2.591469 


-0.004247 


yr 


0.566253 


0.107588 


0.567375 


0.010426 


sin ^f 2 


0.2210 


0.001600 


0.2210 


-0.0001 


sin^?3/10-^ 


28.5820 


5.000000 


28.5803 


-0.0003 


sin^lg/lO-^ 


33.6305 


1.300000 


33.6299 


-0.0005 


5« 


60.0207 


14.000000 


59.9564 


-0.0046 


{m\,)/lQ-\eVf 


5.3851 


0.570823 


5.3850 


-0.0002 


{ml,)/lQ-\eVf 


1.7467 


0.349345 


1.7466 


-0.0004 


sin^ ^^2 


0.2880 


0.057607 


0.2880 


-0.0008 


sin^ ^^3 


0.4595 


0.137850 


0.4597 


0.0014 


^f3(rfe^rees) 


3.7 


3.7 


3.44 




Eigenvalues{Ay) 


0.119568 


0.119717 


0.125583 




Eigenvalues{Ag) 


0.119570 


0.119714 


0.125569 




Eigenvalues{Ap) 


0.149527 


0.149678 


0.155646 




Eigenvalues{Ae) 


0.149568 


0.149724 


0.155689 




Eigenvalues{AQ) 


0.106161 


0.106308 


0.112283 




Eigenvalues{AL) 


0.136273 


0.136419 


0.142375 




An, Ah 


14.262526 


11.287996 






tti 


0.8067 + 0.0000i 


ai 


0.9095 + O.OOOOi 




^2 


0.0179 -0.0131i 


(52 


0.0280 - 0.0024i 




as 


-0.0364 - 0.0057i 


as 


-0.0315 + 0.0016i 




q;4 


0.0024 + 0.5830i 


(54 


0.0171 - 0.4055i 




as 


0.0486 - 0.0185i 


"5 


0.0246 - 0.0388i 




"6 


-0.0287 - 0.0629i 


"6 


-0.0384 - 0.0525i 





Table 24: I-2-bFit with xx = \/j2}liiOi - Oi)^/6f = 0.0232. Target values, at Mx of the 
fermion yukawa couphngs and mixing parameters, together with the estimated uncertainties, 
achieved values and pulls. See caption to Table 3 for explanation. 
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Parameter 


SM{Mz) 


m'^^'^iMz) 


^MSSM _ (^ ^ Am)^^^(Mz) 


md/lO"^ 


2.90000 


0.52473 


3.07224 


m,/10^3 


55.00000 


9.79652 


57.51101 


rrib 


2.90000 


3.24592 


2.90989 


me/10-3 


0.48657 


0.47320 


0.48644 


rrif. 


0.10272 


0.09578 


0.10240 


m-r 


1.74624 


1.74402 


1.74715 


m„/10-3 


1.27000 


1.10467 


1.27126 


rric 


0.61900 


0.53890 


0.62018 


rat 


172.50000 


148.26724 


172.47971 



Table 25: I-2-c: Values of standard model fermion masses in GeV at Mz compared with the 
masses obtained from values of GUT derived yukawa couplings run down from M^ to Mz both 

before and after threshold corrections. Fit with xz = \/J2i=i{^ 
0.0712. 



,MSSM 



m 



SM\2 



/{m. 



MSSM\2 



Parameter 


Value 


Parameter 


Value 


Ml 


104.66 


M-^, 


3630.75 


M2 


353.60 


M^, 


3627.51 


M3 


268.96 


M-n, 


20705.08 


^H 


1887.68 


.0(0 
^11 


-110809.91 


^h 


265.45 


40(0 

^22 


-110677.57 


^3 


11488.56 


.0(0 
^33 


-67135.87 


Ml, 


6440.25 


^11 


-133970.62 


Ml, 


6303.95 


^22 


-133969.80 


Ml, 


10623.40 


^33 


-66830.19 


^. 


825.29 


.0(d) 
^11 


-110623.57 


M-,, 


819.67 


^22 


-110622.83 


^2, 


42238.45 


40(d) 
^33 


-37954.46 


Mq, 


5443.89 


tan/3 


50.00 


Mq, 


5442.40 


^Ji{Mz) 


112330.33 


Mq, 


33482.05 


B{Mz) 


2.1590 X 10^ 


Ml 


-1.2007 X 10^° 


Ml 


-1.3857 X 10^0 



Table 26: I-2-d: Values (GeV) in of the soft Susy parameters at Mz (evolved from the soft 
SUGRY-NUHM parameters at Mx)- See caption to Table 5 for explanation. 
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Field 


MassiGeV) 


Ma 


268.96 


M,± 


353.59,112330.39 


M^o 


104.66, 353.59, 112330.36, 112330.37 


M, 


6439.906, 6303.596, 10623.190 


Mg 


1888.22, 6440.43, 254.80, 6304.74, 10423.50, 11670.42 


Mu 


3630.56, 5443.61, 3627.28, 5442.15, 20701.92, 33484.59 


M;i, 


825.70, 5444.24, 820.02, 5442.76, 33474.52, 42244.47 


Ma 


328626.38 


Mh± 


328626.39 


Mho 


328626.36 


M,o 


112.60 



Table 27: I-2-e: Spectra of supersymmetric partners calculated ignoring generation mixing effects. 
See caption to Table 6 for explanation. 



Field 


Mass{GeV) 


Ma 


269.14 


M,± 


353.67,112325.66 


M,a 


104.69, 353.67, 112325.64, 112325.65 


M, 


6303.72,6440.01,10625.538 


Mg 


296.70, 1887.75, 6294.24, 6440.47, 10487.30, 11631.85 


M^ 


3627.98, 3631.26, 5423.60, 5444.05, 20695.60, 33484.45 


M,: 


822.45, 828.15, 5424.23, 5444.67, 33474.38, 42241.87 


Ma 


328751.52 


Mh± 


328751.53 


Mho 


328751.50 


Mho 


112.62 



Table 28: I-2-f: Spectra of supersymmetric partners calculated including generation mixing 
effects. See caption to Table 7 for explanation. 
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Parameter 


Value 


Field 


Masses 






[5f/(3),5f/(2),r] 


{Units ofW^GeV) 


Xx 


0.0483 


A 1,1, 4] 


809.66 


Xz 


0.0065 


5[6,2,5/3] 


0.2796 


/ill/10-6 


0.0210 


C[8,2,l] 


43.65,387.28,404.49 


/i22/10-^ 


0.0353 


^[3,2,7/3] 


42.78,445.41,456.44 


/i33 


0.0026 


^[3,2,1/3] 


0.45,32.26,32.26 


/11/IO-6 


0.0781 -0.1376i 




34.752,488.08,543.11 


/12/IO-6 


-2.0155 - 0.0665i 


F[l,l,2] 


7.60,7.60 


/13/IO-5 


0.0607 + 0.0478i 




31.02,407.44 


/22/IO-5 


6.5549 - 4.8873i 


^[1,1,0] 


0.070,0.55,0.69 


/23/10-^ 


1.9973 + 2.3436i 




0.692,30.02,30.20 


/33/10-=^ 


-1.0174 + 0.4408i 


/i[l,2,l] 


1.176,25.35,41.75 


^12/10-^ 


0.0608 + 0.1221i 




681.54,704.75 


^13/10-^ 


-0.0465 + 1.7512i 


/[3, 1,10/3] 


0.96 


^23/10-^ 


6.3552 + 5.8153i 


J[3, 1,4/3] 


1.069,17.56,17.56 


A/10-2 


-0.4664 - 0.9791i 




54.25,474.36 


V 


-10.4880 + 2.4932i 


^[3,1,8/3] 


62.72,581.65 


P 


0.6621 - 2.2887i 


L[6, 1,2/3] 


29.75,931.63 


k 


0.0151 -0.0820i 


M[6, 1,8/3] 


942.87 


C 


1.6335 + 0.5734i 


Ar[6, 1,4/3] 


939.66 


c 


0.9799 + 0.4779i 


0[1,3,2] 


1814.58 


m/W^GeV 


0.03 


P[3,3,2/3] 


18.17,1405.68 


nio/lQ^^GeV 


-25.800e-^^''^(^) 


g[8,3,0] 


0.795 


7 


3.91 


i?[8, 1,0] 


0.31,1.19 


7 


-2.8770 


5[1,3,0] 


1.3505 


X 


0.9444 + 0.6732i 


t[3, 1,2/3] 


0.88,24.06,59.01,100.44 


Ax 


1.28 




307.25,415.94,8824.64 


Ag 


1.218 


^[3,3,4/3] 


1.139 


Aa3{Mz) 


-0.013 


\/[l,2,3] 


0.808 


{M'^yW^GeV} 


0.01,14.15,596.40 


Vr[6,3,2/3] 


1079.74 


{M^j/10-^^eV} 


0.4641,865.17,36452.94 


X[3,2,5/3] 


0.273,34.505,34.505 


M^{meV) 


2.14,7.63,42.49 


r[6,2,l/3] 


0.34 


{Evals[f]}/10-^ 


0.15,283.31,11937.12 


Z[8,l,2] 


1.19 


Soft parameters 


mi = -193.241 

2 


mo = 4908.909 


Ao = -1.6043 X 10^ 


at Mx 


/i = 1.2834 X 10^ 


B = -1.1313 X 10^° 


tan/3 = 51.5000 




M| = -1.2330 X 10^° 


M2j = -1.1748 X 10^° 


Rbr = 2.6024 


Max{\LABCD\, \Rabcd ) 


4.9596 X 10-23 GeV-i 




S^l 



Table 29: I-3-a Fit : Column 1 contains values of the NMSGUT-SUGRY-NUHM parameters at 
Mx derived from an accurate fit to all 18 fermion data and compatible with RG constraints. See 
caption to Table 2 for explanation. 
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Parameter 


Target = Oi 


Uncert. = 6i 


Achieved = Oi 


Pull = (0, - 0,)/6, 


?/./io-6 


2.084307 


0.796205 


2.084498 


0.000240 


l/c/10-' 


1.016016 


0.167643 


1.017365 


0.008048 


Vt 


0.383569 


0.015343 


0.383558 


-0.000709 


l/d/10-^ 


7.238813 


4.220228 


7.247117 


0.001968 


l/s/10-3 


1.372490 


0.647815 


1.376178 


0.005694 


Vb 


0.505108 


0.262151 


0.505457 


0.001331 


1/e/lO-^ 


1.283357 


0.192503 


1.285487 


0.011068 


1/m/10-^ 


2.637591 


0.395639 


2.639897 


0.005827 


yr 


0.575157 


0.109280 


0.570576 


-0.041920 


sin ^f 2 


0.2210 


0.001600 


0.2210 


0.0002 


sin^?3/10-^ 


28.8487 


5.000000 


28.8604 


0.0024 


sin^lg/lO-^ 


33.9439 


1.300000 


33.9429 


-0.0008 


5« 


60.0208 


14.000000 


60.0314 


0.0008 


{m\,)/lQ-\eVf 


5.3665 


0.568846 


5.3642 


-0.0040 


{ml,)/lQ-\eVf 


1.7445 


0.348909 


1.7468 


0.0065 


sin^ ^^2 


0.2875 


0.057499 


0.2869 


-0.0109 


sin^ ^^3 


0.4576 


0.137268 


0.4591 


0.0114 


^f3(rfe^rees) 


3.7 


3.7 


4.55 




Eigenvalues{Ay) 


0.053857 


0.053861 


0.053868 




Eigenvalues{Ag) 


0.051644 


0.051649 


0.051655 




Eigenvalues{Ap) 


0.058766 


0.058770 


0.058776 




Eigenvalues{Ae) 


0.065404 


0.065409 


0.065415 




Eigenvalues{AQ) 


0.050688 


0.050692 


0.050697 




Eigenvalues{AL) 


0.060022 


0.060026 


0.060032 




An, Ah 


77.639981 


63.294285 






tti 


0.7919 - O.OOOOi 


«! 


0.8754 - O.OOOOi 




^2 


0.0651 + 0.0139i 


(52 


0.0558 + 0.0625i 




as 


-0.0416 - 0.039H 


as 


-0.0561 -0.0179i 




q;4 


-0.4597 + 0.2078i 


(54 


0.3572 - 0.0342i 




as 


0.1484 + 0.0808i 


"5 


0.0685 + 0.0189i 




"6 


0.1097 -0.2650i 


"6 


0.1198 -0.2740i 





Table 30: I-3-b: Fit with xx = \/j2}liiOi - Oi)^/5f = 0.0483. Target values, at Mx of the 
fermion yukawa couplings and mixing parameters, together with the estimated uncertainties, 
achieved values and pulls. See caption to Table 3 for explanation. 
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Parameter 


SM{Mz) 


m'^^'^iMz) 


^MSSM _ (^ ^ Amf^'^iMz) 


md/lO"^ 


2.90000 


0.62818 


2.90275 


m,/10^3 


55.00000 


11.88721 


54.95005 


rrib 


2.90000 


3.07412 


2.89882 


me/10-3 


0.48657 


0.47199 


0.48604 


rrif. 


0.10272 


0.09706 


0.10272 


m-r 


1.74624 


1.73492 


1.73563 


m„/10-3 


1.27000 


1.09437 


1.26906 


rric 


0.61900 


0.53416 


0.61943 


rat 


172.50000 


148.47236 


172.66514 



Table 31: I-3-c: Values of standard model fermion masses in GeV at Mz compared with the 
masses obtained from values of GUT derived yukawa couplings run down from M^ to Mz both 

before and after threshold corrections. Fit with xz = \/J2i=i{^ 
0.0065. 



,MSSM 



m 



SM\2 



/{m. 



MSSM\2 



Parameter 


Value 


Parameter 


Value 


Ml 


147.27 


M-^, 


6318.67 


M2 


406.23 


M^, 


6317.40 


M3 


598.97 


M-n, 


21142.35 


^H 


1614.43 


.0(0 
^11 


-98630.83 


^h 


320.14 


4 0(0 

^22 


-98506.40 


^3 


16050.48 


.0(0 
^33 


-59069.03 


Ml, 


7827.10 


^11 


-117145.75 


Ml, 


7747.15 


^22 


-117144.96 


Ml, 


13923.16 


^33 


-59199.13 


^. 


3334.46 


.0(d) 
^11 


-98885.51 


M-,, 


3333.11 


^22 


-98884.62 


^2, 


36975.96 


40(d) 
^33 


-37795.22 


Mq, 


6120.46 


tan/3 


51.50 


Mq, 


6119.44 


^Ji{Mz) 


101035.83 


Mq, 


30209.06 


B{Mz) 


1.7444 X 10^ 


Ml 


-1.0097 X 10^° 


Ml 


-1.1107 X 10^° 



Table 32: I-3-d: Values (GeV) in of the soft Susy parameters at Mz (evolved from the soft 
SUGRY-NUHM parameters at Mx)- See caption to Table 5 for explanation. 
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Field 


MassiGeV) 


Ma 


598.97 


M,± 


406.23, 101035.90 


M^o 


147.27,406.23, 101035.87, 101035.87 


M, 


7826.822, 7746.869, 13922.999 


Mg 


1615.07, 7827.26, 316.41, 7747.59, 13877.17, 16090.39 


Mu 


6120.21, 6318.57, 6119.06, 6317.42, 21138.66, 30212.30 


M^, 


3334.56, 6120.76, 3333.19, 6119.76, 30199.68, 36983.67 


Ma 


299785.40 


Mh± 


299785.42 


Mho 


299785.39 


M,o 


114.58 



Table 33: I-3-e: Spectra of supersymmetric partners calculated ignoring generation mixing effects. 
See caption to Table 6 for explanation. 



Field 


Mass{GeV) 


Ma 


599.14 


M,± 


406.47, 100968.77 


M,a 


147.34,406.47, 100968.75, 100968.75 


M, 


7744.59, 7824.46, 14034.822 


Me 


393.14, 1632.84, 7745.08, 7824.89, 13992.95, 16284.21 


Mu 


6108.76, 6117.73, 6311.63, 6312.80, 21005.60, 30236.11 


Ma 


3331.17, 3332.57, 6109.44, 6118.29, 30223.69, 37095.48 


Ma 


300060.73 


Mh± 


300060.74 


Mho 


300060.72 


Mho 


114.85 



Table 34: I-3-f: Spectra of supersymmetric partners calculated including generation mixing 
effects. See caption to Table 7 for explanation. 



59 



Parameter 


Value 


Field 


Masses 






[SU{'i),SU{2),Y] 


{Units ofW^GeV) 


Xx 


0.0771 


Al,l,4 


168.86 


Xz 


0.0225 


5[6,2,5/3] 


0.3197 


/ill/10-6 


0.0156 


C[8,2, 1] 


5.54,205.24,206.53 


/i22/10-^ 


-0.0131 


Dp, 2, 7/3] 


4.90,110.54,138.09 


/i33 


0.0033 


^[3,2,1/3] 


0.38,3.89,6.18 


/11/IO-6 


0.0198 + 0.1542i 




6.184,165.12,192.29 


/12/IO-6 


-1.2260 + 0.3835i 


F[l,l,2] 


1.15,1.15 


/13/IO-5 


0.3221 + 0.1153i 




3.20,85.63 


/22/IO-5 


5.6438 - 6.8723i 


G[1,1,0] 


0.089,0.28,0.28 


/23/10-^ 


2.0515 + 3. 1055i 




0.722,16.74,17.12 


/33/10-=^ 


-1.0661 + 0.7067i 


/i[l,2,l] 


1.325,3.74,5.08 


^12/10-^ 


0.0506 - 0.0960i 




145.15,190.66 


^13/10-^ 


-1.5789 - 2.0450i 


/[3, 1,10/3] 


1.29 


^23/10-^ 


2.9105 + 6.4235i 


J[3, 1,4/3] 


1.261,3.10,3.10 


A/10-2 


-6.6788 - 1.6361i 




8.01,130.00 


V 


-14.9636 + 0.9886i 


i^[3, 1,8/3] 


8.83,154.03 


P 


0.4151 - 1.6189i 


L[6, 1,2/3] 


4.56,252.21 


k 


0.0254 - 0.0733i 


M[6, 1,8/3] 


258.11 


C 


0.9369 + 0.5584i 


Ar[6, 1,4/3] 


249.35 


c 


1.5736 -0.4378i 


0[1,3,2] 


456.61 


m/W^GeV 


0.03 


P[3,3,2/3] 


1.32,401.98 


nio/lQ^^GeV 


-3.749e-*^'^»(^) 


g[8,3,0] 


1.138 


7 


4.49 


i?[8, 1, 0] 


0.34,1.34 


7 


-2.4523 


5[1,3,0] 


1.5523 


X 


0.8697 + 0.4711i 


t[3, 1,2/3] 


0.94,3.03,8.58,10.22 


Ax 


0.45 




162.05,171.60,3725.03 


Ag 


19.546 


f/[3,3,4/3] 


1.367 


Aa3{Mz) 


-0.008 


\/[l,2,3] 


0.875 


{M'^yW^GeV} 


0.01,6.59,279.67 


1^(6,3,2/3] 


382.85 


{M^j/10-^^eV} 


5.5373,3939.17,167187.45 


X[3,2,5/3] 


0.293,6.732,6.732 


M^{meV) 


2.30,7.53,41.51 


F[6,2,l/3] 


0.36 


{Evals[f]}/10-^ 


0.46,326.47,13856.04 


Z[8,l,2] 


1.34 


Soft parameters 


mi = -221.928 

2 


mo = 5063.109 


Ao = -1.4400 X 10^ 


at Mx 


/i = 1.2715 X 10^ 


B = -9.2002 X 10^ 


tan/3 = 51.0000 




M| = -1.4452 X 10^° 


M2j = -1.3796 X 10^° 


Rb^ = 3.6620 


Max{\LABCD\, \Rabcd ) 


4.5150 X 10-22 GeV-i 




S/J. 



Table 35: I-4-a Fit : Column 1 contains values of the NMSGUT-SUGRY-NUHM parameters at 
Mx derived from an accurate fit to all 18 fermion data and compatible with KG constraints. See 
caption to Table 2 for explanation. 
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Parameter 


Target = Oi 


Uncert. = 5i 


Achieved = Oi 


Pull = (0, - 0,)/6, 


?/./io-6 


2.099660 


0.802070 


2.099948 


0.000359 


l/c/10-' 


1.023488 


0.168875 


1.022975 


-0.003035 


Vt 


0.368075 


0.014723 


0.368096 


0.001417 


l/d/10-^ 


6.248810 


3.643056 


6.265778 


0.004658 


l/s/10-3 


1.185571 


0.559590 


1.201681 


0.028788 


Vb 


0.424469 


0.220299 


0.436945 


0.056636 


1/e/lO-^ 


1.173101 


0.175965 


1.171800 


-0.007392 


1/m/10-^ 


2.478866 


0.371830 


2.481095 


0.005994 


yr 


0.527577 


0.100240 


0.523402 


-0.041647 


sin ^f 2 


0.2210 


0.001600 


0.2210 


-0.0001 


sin^?3/10-^ 


29.4359 


5.000000 


29.4173 


-0.0037 


sin^lg/lO-^ 


34.6341 


1.300000 


34.6323 


-0.0013 


5« 


60.0209 


14.000000 


60.0406 


0.0014 


{m\,)/lQ-\eVf 


5.1465 


0.545532 


5.1456 


-0.0016 


{ml,)/lQ-\eVf 


1.6659 


0.333190 


1.6666 


0.0019 


sin^ ^^2 


0.2885 


0.057707 


0.2883 


-0.0043 


sin^ ^^3 


0.4613 


0.138393 


0.4616 


0.0023 


^f3(rfe^rees) 


3.7 


3.7 


4.29 




Eigenvalues{Ay) 


0.011860 


0.011862 


0.011863 




Eigenvalues{Ag) 


0.008409 


0.008412 


0.008413 




Eigenvalues{Ap) 


0.002387 


0.002389 


0.002390 




Eigenvalues{Ae) 


0.012738 


0.012740 


0.012741 




Eigenvalues{AQ) 


0.014518 


0.014520 


0.014520 




Eigenvalues{AL) 


0.011946 


0.011948 


0.011948 




An, Ah 


59.178650 


51.533539 






tti 


0.7545 - O.OOOOi 


tti 


0.8064 - O.OOOOi 




^2 


0.0491 + 0.0409i 


(52 


-0.0007 + 0.0768i 




as 


0.0199 - 0.0535i 


as 


-0.0051 - 0.0324i 




q;4 


-0.5232 - 0.0719i 


(54 


0.3297 + 0.1532i 




as 


0.2001 + 0.0065i 


"5 


0.1175 -0.0307i 




"6 


0.0512 -0.3191i 


"6 


0.0653 - 0.4377i 





Table 36: I-4-b: Fit with xx = \/j2}liiOi - Oi)^/5f = 0.0771. Target values, at Mx of the 
fermion yukawa couplings and mixing parameters, together with the estimated uncertainties, 
achieved values and pulls. See caption to Table 3 for explanation. 
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Parameter 


SM{Mz) 


m'^^'^iMz) 


^MSSM _ (^ ^ Am)^^^(Mz) 


md/lO"^ 


2.90000 


0.57869 


2.88556 


m,/10^3 


55.00000 


11.09336 


55.29916 


rrib 


2.90000 


2.96813 


2.93512 


me/10-3 


0.48657 


0.45919 


0.48201 


rrif. 


0.10272 


0.09719 


0.10198 


m-r 


1.74624 


1.72539 


1.72437 


m„/10-3 


1.27000 


1.11423 


1.27070 


rric 


0.61900 


0.54254 


0.61873 


ITit 


172.50000 


147.09707 


172.15913 



Table 37: I-4-c: Values of standard model fermion masses in GeV at Mz compared with the 
masses obtained from values of GUT derived yukawa couplings run down from M^ to Mz both 

before and after threshold corrections. Fit with xz = \/J2i=i{^ 
0.0225. 



,MSSM 



m 



SM\2 



/{m. 



MSSM\2 



Parameter 


Value 


Parameter 


Value 


Ml 


104.16 


M-^, 


6018.96 


M2 


302.05 


M^, 


6018.69 


M3 


350.55 


M-n, 


32268.70 


^H 


115.35 


.0(0 
^11 


-92737.91 


^h 


123.64 


40(0 


-92630.01 


^3 


27768.92 


.0(0 
^33 


-58880.24 


Ml, 


8173.98 


^11 


-105676.88 


Ml, 


8174.04 


^22 


-105676.19 


Ml, 


21289.08 


^33 


-55613.06 


^. 


1701.52 


.0(d) 
^11 


-92749.44 


^^ 


1702.61 


^22 


-92748.73 


^^3 


40148.17 


40(d) 
^33 


-41722.79 


Mq, 


5628.19 


tan/3 


51.00 


Mq, 


5628.21 


^Ji{Mz) 


106300.71 


Mq, 


36452.35 


B{Mz) 


1.8340 X 10^ 


Ml 


-1.1556 X 10^° 


Ml 


-1.2234 X 10^° 



Table 38: I-4-d: Values (GeV) in of the soft Susy parameters at Mz (evolved from the soft 
SUGRY-NUHM parameters at Mx)- See caption to Table 5 for explanation. 
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Field 


^^055(^6^) 


Ma 


350.55 


M,± 


302.05, 106300.78 


M^o 


104.16, 302.05, 106300.75, 106300.76 


M, 


8173.708, 8173.768, 21288.979 


Mg 


123.89, 8174.12, 114.16, 8174.45, 21282.54, 27774.02 


Mu 


5627.92, 6018.85, 5627.87, 6018.64, 32265.13, 36456.04 


M,- 


1701.72, 5628.53, 1702.77, 5628.56, 36439.70, 40159.71 


M^ 


305894.52 


Mh± 


305894.54 


Mho 


305894.52 


M,o 


128.29 



Table 39: I-4-e: Spectra of supersymmetric partners calculated ignoring generation mixing effects. 
See caption to Table 6 for explanation. 



Field 


Mass{GeV) 


Ma 


350.89 


M,± 


302.21, 106276.46 


M,o 


104.22, 302.21, 106276.43, 106276.44 


M, 


8173.67,8173.76,21295.664 


Me 


89.25, 117.38, 8174.09, 8174.44, 21289.49, 27784.15 


Mu 


5627.05, 5648.66, 6018.00, 6018.22, 32278.73, 36462.50 


M,: 


1699.10, 1700.14, 5627.66, 5649.32, 36446.20, 40166.59 


Ma 


306124.85 


Mh± 


306124.86 


Mho 


306124.84 


Mho 


128.38 



Table 40: I-4-f: Spectra of supersymmetric partners calculated including generation mixing 
effects. See caption to Table 7 for explanation. 
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Parameter 


Value 


Field 


Masses 






[SU{'i),SU{2),Y] 


{Units ofW^GeV) 


Xx 


0.0050 


A[l,l,4 


948.38 


Xz 


0.0038 


S[6,2,5/3] 


0.2867 


/ill/10-6 


0.0200 


C[8,2, 1] 


54.37,436.41,458.35 


/i22/10-^ 


0.0348 


D[3,2,7/3] 


52.04,521.13,532.76 


/i33 


0.0027 


^[3,2,1/3] 


0.47,38.02,38.02 


/11/IO-6 


0.0796 -0.1364i 




43.611,563.34,624.40 


/12/IO-6 


-2.0359 - 0.0639i 


F[l,l,2] 


9.05,9.05 


/13/IO-5 


0.0591 + 0.0506i 




36.80,477.20 


/22/IO-5 


6.5418 - 4.7649i 


^[1,1,0] 


0.072,0.56,0.77 


/23/10-^ 


1.9886 + 2.3416i 




0.767,32.75,32.92 


/33/10-=^ 


-1.0141 + 0.4463i 


/i[l,2,l] 


1.220,30.91,50.43 


^12/10-^ 


0.0574 + 0.1266i 




793.49,821.05 


^13/10-^ 


-0.0387 + 1.6989i 


J[3, 1,10/3] 


0.98 


^23/10-^ 


6.4300 + 5.9592i 


J[3, 1,4/3] 


1.093,20.78,20.78 


A/10-2 


-0.3569 - 0.8704i 




64.58,550.84 


V 


-10.3588 + 2.2825i 


^[3,1,8/3] 


74.44, 677.26 


P 


0.6169 - 2.2892i 


L[6, 1,2/3] 


36.04,1081.79 


k 


0.0102 - 0.0863i 


M[6, 1,8/3] 


1094.90 


C 


1.5786 + 0.4937i 


Ar[6, 1,4/3] 


1091.44 


c 


1.0902 + 0.4822i 


0[1,3,2] 


2107.18 


m/W^GeV 


0.03 


P[3,3,2/3] 


21.00,1627.16 


nio/lQ^^GeV 


-31.348e-*^''^(^) 


g[8,3,0] 


0.808 


7 


3.89 


i?[8, 1,0] 


0.32,1.22 


7 


-2.8834 


^[1,3,0] 


1.3822 


X 


0.9453 + 0.6806i 


t[3, 1,2/3] 


0.90,29.04,70.15,120.87 


Ax 


1.36 




348.42,478.70,10025.59 


Ag 


-0.114 


f/[3,3,4/3] 


1.164 


Aa3{Mz) 


-0.010 


\/[l,2,3] 


0.830 


{M'^yW^GeV} 


0.00,15.70,660.27 


Vr[6,3,2/3] 


1241.12 


{M^j/10-^^eV} 


0.2212,725.65,30521.19 


X[3,2,5/3] 


0.280,40.637,40.637 


M^{meV) 


2.90,7.98,43.17 


F[6,2,l/3] 


0.35 


{Evals[f]}/10-^ 


0.09,283.52,11924.94 


Z[8,l,2] 


1.21 


Soft parameters 


mi = -170.786 

2 


mo = 1582.439 


Ao = -7.8096 X 10^ 


at Mx 


/i = 6.8314 X 10^ 


B = -3.0341 X 10^ 


tan/3 = 51.5000 




M| = -3.4294 X 10^ 


M^j = -3.4202 X 10^ 


Rb. = 2.4959 


Max{\LABCD\, \Rabcd ) 


4.4479 X 10^23 GeV-i 




s^ 



Table 41: II-3-a : Column 1 contains values of the NMSGUT-SUGRY-NUHM parameters at 
Mx derived from an accurate fit to all 18 fermion data and compatible with KG constraints. See 
caption to Table 2 for explanation. 



64 



Parameter 


Target = Oi 


Uncert. = 5i 


Achieved = Oi 


Pull = (0, - 0,)/6, 


?/./io-6 


2.100774 


0.802496 


2.100767 


-0.000009 


l/c/10-' 


1.023520 


0.168881 


1.023538 


0.000104 


Vt 


0.393877 


0.015755 


0.393880 


0.000182 


l/d/10-^ 


7.305382 


4.259037 


7.317435 


0.002830 


l/s/10-3 


1.384369 


0.653422 


1.382982 


-0.002123 


Vb 


0.516686 


0.268160 


0.516774 


0.000329 


1/e/lO-^ 


1.317579 


0.197637 


1.317708 


0.000656 


1/m/10-^ 


2.780859 


0.417129 


2.781086 


0.000546 


yr 


0.583110 


0.110791 


0.582736 


-0.003378 


sin ^f 2 


0.2210 


0.001600 


0.2210 


-0.0001 


sin^?3/10-^ 


28.7161 


5.000000 


28.7167 


0.0001 


sin^lg/lO-^ 


33.7880 


1.300000 


33.7881 


0.0001 


5« 


60.0207 


14.000000 


60.0113 


-0.0007 


{m\,)/lQ-\eVf 


5.5328 


0.586474 


5.5328 


0.0000 


{ml,)/lQ-\eVf 


1.8000 


0.359999 


1.8000 


0.0000 


sin^ ^^2 


0.2873 


0.057461 


0.2873 


0.0000 


sin^ ^^3 


0.4569 


0.137073 


0.4569 


0.0000 


^f3(rfe^rees) 


3.7 


3.7 


6.88 




Eigenvalues{Ay) 


0.057990 


0.057995 


0.058002 




Eigenvalues{Ag) 


0.055908 


0.055913 


0.055920 




Eigenvalues{Ap) 


0.064334 


0.064339 


0.064345 




Eigenvalues{Ae) 


0.070580 


0.070585 


0.070592 




Eigenvalues{AQ) 


0.054281 


0.054286 


0.054292 




Eigenvalues{AL) 


0.064789 


0.064794 


0.064799 




An, Ah 


78.095163 


63.896767 






tti 


0.7879 - O.OOOOi 


«! 


0.8695 - O.OOOOi 




^2 


0.0649 + 0.0142i 


(52 


0.0583 + 0.0634i 




as 


-0.0413 - 0.0387i 


as 


-0.0548 - 0.0158i 




q;4 


-0.4579 + 0.2315i 


(54 


0.3691 - 0.0278i 




as 


0.1505 + 0.0742i 


"5 


0.0719 + 0.0175i 




"6 


0.0984 - 0.2655i 


"6 


0.1081 -0.2814i 





Table 42: II-3-b: Fit with xx = \/E}li{Oi - Oi)^/5f = 0.0050. Target values, at Mx of 
the fermion yukawa couplings and mixing parameters, together with the estimated uncertainties, 
achieved values and pulls. See caption to Table 3 for explanation. 
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Parameter 


SM{Mz) 


m'^^'^iMz) 


^MSSM _ (^ ^ Am)^^^(Mz) 


md/lO"^ 


2.90000 


0.62793 


2.90404 


m,/10^3 


55.00000 


11.86773 


54.92429 


rrib 


2.90000 


3.08544 


2.90008 


me/10-3 


0.48657 


0.48023 


0.48631 


rrif. 


0.10272 


0.10130 


0.10265 


m-r 


1.74624 


1.74579 


1.74474 


m„/10-3 


1.27000 


1.09336 


1.26989 


rric 


0.61900 


0.53271 


0.61905 


mt 


172.50000 


149.44924 


171.97519 



Table 43: II-3-c: Values of standard model fermion masses in GeV at Mz compared with the 
masses obtained from values of GUT derived yukawa couplings run down from M^ to Mz both 

before and after threshold corrections. Fit with xz = \/J2i=i{''^i 
0.0038. 



,MSSM 



m 



SM\2 



/{m. 



MSSM\2 



Parameter 


Value 


Parameter 


Value 


Ml 


43.51 


M^, 


205.86 


M2 


144.40 


M-n, 


200.37 


M3 


89.92 


Mn, 


15100.44 


^\ 


2153.08 


Ml) 
^11 


-47638.60 


^% 


2072.60 


40(0 

^22 


-47572.19 


^3 


12763.49 


40(0 


-28011.67 


Mu 


2923.85 


40W 
^11 


-56421.60 


Mu 


2894.78 


4OH 
^22 


-56421.22 


Ml, 


9407.34 


^33 


-27834.22 


% 


455.68 


.0(d) 

^11 


-47530.07 


M-- 

d2 


455.15 


^22 


-47529.64 


^Ja 


20922.14 


40(d) 
^33 


-17493.87 


Mq, 


2445.28 


tan/3 


51.50 


Mq, 


2445.00 


KMz) 


52905.86 


Mq, 


18336.34 


B{Mz) 


4.0145 X 10^ 


M% 


-2.6746 X 10^ 


Ml 


-3.0610 X 10^ 



Table 44: II-3-d: Values (GeV) in of the soft Susy parameters at Mz (evolved from the soft 
SUGRY-NUHM parameters at Mx)- See caption to Table 5 for explanation. 
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Field 


MassiGeV) 


Ma 


89.92 


M,± 


144.39,52905.99 


M^o 


43.51, 144.39, 52905.94, 52905.95 


M, 


2923.095, 2894.022, 9407.108 


Me 


2153.55, 2924.25, 2068.45, 2898.51, 9391.01, 12775.68 


Mu 


202.54, 2444.65, 196.56, 2444.41, 15095.44, 18341.55 


M,- 


456.42, 2446.04, 455.69, 2445.80, 18317.31, 20938.90 


Ma 


143814.30 


Mh± 


143814.33 


Mho 


143814.29 


M,o 


123.80 



Table 45: II-3-e: Spectra of supersyninietric partners calculated ignoring generation mixing 
effects. See caption to Table 6 for explanation. 



Field 


Mass{GeV) 


Ma 


90.00 


Mx± 


144.43,52903.84 


M^o 


43.52, 144.43, 52903.79, 52903.80 


M, 


2894.06,2923.13,9407.538 


Me 


2065.45, 2153.46, 2899.51, 2924.33, 9392.27, 12775.95 


Mu 


198.31, 204.25, 2444.78, 2451.67, 15094.69, 18340.05 


Mi 


456.03, 456.94, 2446.17, 2453.03, 18315.87, 20938.32 


Ma 


143872.09 


Mh± 


143872.11 


Mho 


143872.08 


Mho 


123.91 



Table 46: II-3-f: Spectra of supersymmetric partners calculated including generation mixing 
effects. See caption to Table 6 for explanation. 
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Parameter 


Value 


Field 


Masses 






[5f/(3),5f/(2),r] 


{Units ofW^GeV) 


Xx 


0.0061 


A 1,1, 4] 


137.29 


Xz 


0.0179 


5(6,2,5/3] 


0.3082 


/ill/10-6 


0.0277 


C[8,2,l] 


4.55,166.93,168.04 


/i22/10-^ 


-0.0081 


D[3,2,7/3] 


3.93,95.60,117.42 


/i33 


0.0034 


E[3,2,l/3] 


0.37,3.21,5.35 


/11/IO-6 


0.0242 + 0.1635i 




5.346, 138.97, 160.78 


/12/IO-6 


-1.6733 + 0.4769i 


F[l,l,2] 


0.96,0.96 


/13/IO-5 


0.3397 + 0.118n 




2.66,69.79 


/22/IO-5 


5.3022 - 7.2813i 


^[1,1,0] 


0.087,0.26,0.26 


/23/10-^ 


2.0313 + 3. 122n 




0.698,14.94,15.31 


/33/10-=^ 


-1.0321 + 0.6903i 


/i[l,2,l] 


1.233,3.11,4.16 


^12/10-^ 


0.0375 -0.1054i 




126.31,161.82 


^13/10-^ 


-1.5603 - 1.8606i 


/[3, 1,10/3] 


1.27 


^23/10-^ 


3.3223 + 6.3071i 


J[3, 1,4/3] 


1.249,2.62,2.62 


A/10-2 


-7.6198 - 1.9488i 




6.73,105.91 


V 


-14.4375 + 1.9014i 


^[3,1,8/3] 


7.33,124.51 


P 


0.3809 -1.5817i 


L[6, 1,2/3] 


3.76,205.78 


k 


0.0229 - 0.0773i 


M[6, 1,8/3] 


210.07 


C 


0.9769 + 0.4323i 


iV[6, 1,4/3] 


203.97 


c 


1.4935 -0.4959i 


0[1,3,2] 


386.98 


m/W^GeV 


0.03 


P[3,3,2/3] 


1.11,339.57 


nio/lQ^^GeV 


-3.093e-*^'^3(^) 


g[8,3,0] 


1.117 


7 


4.63 


i?[8, 1,0] 


0.33,1.33 


7 


-2.3784 


5[1,3,0] 


1.5325 


X 


0.8833 + 0.4731i 


t[3, 1,2/3] 


0.94,2.51,7.17,9.19 


Ax 


0.39 




134.01,142.61,3033.21 


Ag 


20.063 


f/[3,3,4/3] 


1.351 


Aa3{Mz) 


-0.012 


\/[l,2,3] 


0.851 


{M'^yW^GeV} 


0.01,7.43,250.87 


Vr[6,3,2/3] 


320.01 


{M^j/10-^^eV} 


5.8487,6146.10,207404.09 


X[3,2,5/3] 


0.284,5.796,5.796 


M^{meV) 


4.27,8.56,43.09 


F[6,2,l/3] 


0.35 


{Evals[f]}/10-^ 


0.38,400.40,13511.64 


Z[8,l,2] 


1.33 


Soft parameters 


mi = -118.386 

2 


mo = 1773.201 


Ao = -5.9163 X 10^ 


at Mx 


/i = 5.4671 X 10^ 


B = -1.7502 X 10^ 


tan/3 = 51.0000 




M| = -2.5185 X 10^ 


M2j = -2.4818 X 10^ 


Rb. = 3.5550 


Max{\LABCD\, \Rabcd ) 


5.4994 X 10-22 GeV-i 




S/l 



Table 47: II-4-a: Column 1 contains values of the NMSGUT-SUGRY-NUHM parameters at 
Mx derived from an accurate fit to all 18 fermion data and compatible with KG constraints. See 
caption to Table 2 for explanation. 
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Parameter 


Target = Oi 


Uncert. = 5i 


Achieved = Oi 


Pull = (0, - 0,)/6, 


?/./io-6 


2.169056 


0.828579 


2.169068 


0.000015 


l/c/10-' 


1.057321 


0.174458 


1.057338 


0.000102 


Vt 


0.390644 


0.015626 


0.390643 


-0.000072 


l/d/10-^ 


6.611935 


3.854758 


6.614046 


0.000548 


l/s/10-3 


1.255361 


0.592530 


1.256890 


0.002581 


Vb 


0.452379 


0.234785 


0.453431 


0.004481 


1/e/lO-^ 


1.247646 


0.187147 


1.247511 


-0.000719 


1/m/10-^ 


2.635432 


0.395315 


2.635821 


0.000983 


yr 


0.542972 


0.103165 


0.542667 


-0.002953 


sin ^f 2 


0.2210 


0.001600 


0.2210 


0.0001 


sin^?3/10-^ 


29.1115 


5.000000 


29.1102 


-0.0003 


sin^lg/lO-^ 


34.2526 


1.300000 


34.2523 


-0.0002 


5« 


60.0206 


14.000000 


60.0202 


0.0000 


{m\,)/lQ-\eVf 


5.5012 


0.583126 


5.5012 


0.0000 


{ml,)/lQ-\eVf 


1.7832 


0.356636 


1.7832 


0.0000 


sin^ ^^2 


0.2882 


0.057639 


0.2882 


0.0000 


sin^ ^^3 


0.4601 


0.138033 


0.4601 


-0.0001 


^f3(rfe^rees) 


3.7 


3.7 


6.89 




Eigenvalues{Ay) 


0.008011 


0.008014 


0.008014 




Eigenvalues{Ag) 


0.004482 


0.004484 


0.004485 




Eigenvalues{Ap) 


0.002764 


0.002765 


0.002767 




Eigenvalues{Ae) 


0.007822 


0.007824 


0.007825 




Eigenvalues{AQ) 


0.011230 


0.011232 


0.011232 




Eigenvalues{AL) 


0.007511 


0.007513 


0.007513 




An, Ah 


58.805443 


52.078952 






tti 


0.7626 + O.OOOOi 


ai 


0.8079 + O.OOOOi 




^2 


0.0468 + 0.0455i 


(52 


-0.0039 + 0.0795i 




as 


0.0221 - 0.0516i 


as 


-0.0018 - 0.0322i 




q;4 


-0.4860 - 0.102H 


(54 


0.2917 + 0.1584i 




as 


0.2131 + 0.0060i 


as 


0.1288 -0.0330i 




"6 


0.0562 - 0.3402i 


aa 


0.0759 - 0.454H 





Table 48: II-4-b: Fit with xx = \/E}li{Oi - Oi)^/5f = 0.0061. Target values, at Mx of 
the fermion yukawa couplings and mixing parameters, together with the estimated uncertainties, 
achieved values and pulls. See caption to Table 3 for explanation. 
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Parameter 


SM{Mz) 


m'^^'^iMz) 


^MSSM _ (^ ^ Amf^'^iMz) 


md/lO"^ 


2.90000 


0.60092 


2.89886 


m,/10^3 


55.00000 


11.41936 


55.02309 


rrib 


2.90000 


2.98616 


2.90724 


me/10-3 


0.48657 


0.48126 


0.48642 


rrif. 


0.10272 


0.10163 


0.10271 


m-r 


1.74624 


1.74575 


1.74469 


m„/10-3 


1.27000 


1.12845 


1.25425 


rric 


0.61900 


0.55007 


0.61140 


mt 


172.50000 


150.12439 


172.45759 



Table 49: II-4-c: Values of standard model fermion masses in GeV at Mz compared with the 
masses obtained from values of GUT derived yukawa couplings run down from M^ to Mz both 

before and after threshold corrections. Fit with xz = \/J2i=i{''^i 
0.0179. 



,MSSM 



m 



SM\2 



/{m. 



MSSM\2 



Parameter 


Value 


Parameter 


Value 


Ml 


35.00 


M^, 


1355.00 


M2 


110.00 


M-n, 


1354.97 


M3 


82.43 


Mn, 


15017.39 


^\ 


1668.31 


Ml) 
^11 


-37671.35 


^% 


1686.32 


40(0 

^22 


-37622.42 


^3 


12622.76 


40(0 


-23300.14 


Mu 


2873.66 


aO(u) 
^11 


-42600.14 


Mu 


2878.88 


4OH 
^22 


-42599.84 


Ml, 


9306.38 


^33 


-21238.50 


% 


362.45 


.0(d) 
^11 


-37604.10 


M-- 

d2 


363.78 


^22 


-37603.79 


^Ja 


17343.17 


40(d) 
^33 


-16155.25 


Mq, 


2110.14 


tan/3 


51.00 


Mq, 


2110.24 


KMz) 


44372.34 


Mq, 


16267.76 


B{Mz) 


2.6446 X 10^ 


M% 


-1.9734 X 10^ 


Ml 


-2.1327 X 10^ 



Table 50: II-4-d : Values (GeV) in of the soft Susy parameters at Mz (evolved from the soft 
SUGRY-NUHM parameters at Mx)- See caption to Table 5 for explanation. 
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Field 


^^055(^6^) 


Ma 


82.43 


M,± 


109.99,44372.49 


M^o 


35.00, 109.99, 44372.43, 44372.44 


M, 


2872.898,2878.121,9306.146 


Mg 


1668.92, 2874.07, 1683.95, 2881.03, 9294.77, 12631.49 


Mu 


1354.50, 2109.42, 1354.39, 2109.57, 15008.78, 16276.95 


M,- 


363.38, 2111.02, 364.49, 2111.16, 16229.75, 17378.88 


M^ 


116157.09 


Mh± 


116157.12 


Mho 


116157.08 


M,o 


127.41 



Table 51: II-4-e: Spectra of supersyninietric partners calculated ignoring generation mixing 
effects. See caption to Table 6 for explanation. 



Field 


Mass{GeV) 


Ma 


82.49 


M,± 


110.02,44370.87 


M^o 


35.01, 110.02, 44370.81, 44370.82 


M, 


2872.91,2878.14,9306.355 


Me 


1668.80, 1681.53, 2874.11, 2881.63, 9295.61, 12631.53 


Mu 


1354.49, 1354.60, 2109.47, 2123.10, 15008.73, 16275.15 


M,: 


363.58, 364.64, 2111.08, 2124.64, 16228.06, 17378.66 


Ma 


116199.01 


Mh± 


116199.03 


Mho 


116199.00 


Mho 


127.43 



Table 52: II-4-f: Spectra of supersymmetric partners calculated including generation mixing 
effects. See caption to Table 7 for explanation. 
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Parameter 


Value 


Field 


Masses 






[5f/(3),5f/(2),r] 


{Units ofW^GeV) 


Xx 


0.0143 


A 1,1,4] 


16007.39 


Xz 


0.0066 


5[6,2,5/3] 


0.3294 


/ill/10-6 


-0.4946 


C[8,2, 1] 


671.99,2566.01,3875.58 


/i22/10-^ 


0.3413 


^[3,2,7/3] 


1412.23,9043.81,10626.05 


/i33 


0.0143 


E[3,2,l/3] 


0.74,488.82,488.82 


/11/IO-6 


-0.1019 -0.0216i 




1600.508,4925.64,7365.53 


/12/IO-6 


1.6275 - 0.0739i 


F[l,l,2] 


186.39,186.39 


/13/IO-5 


-0.6414 - 0.0365i 




1045.25,8230.98 


/22/IO-5 


-1.8755 + 0.1518i 


^[1,1,0] 


0.072,0.60,3.18 


/23/10-^ 


1.0186 -0.0082i 




3.179,111.33,111.60 


/33/10-=^ 


-0.4301 + 0.017H 


/i[l,2,l] 


1.448,865.06,1334.27 


^12/10-^ 


-2.8487 - 0.6757i 




9342.01,12728.76 


^13/10-^ 


-57.2472 - 43.8940i 


J[3, 1,10/3] 


0.84 


^23/10-^ 


139.1403 + 111.6681? 


J[3, 1,4/3] 


1.001,327.80,327.80 


A/10-2 


0.0590 + 0.0380i 




585.72,6966.80 


V 


-8.6784 - 0.3384i 


^[3,1,8/3] 


1019.21,10139.63 


P 


1.2622 -0.6215i 


L[6, 1,2/3] 


816.96, 12975.96 


k 


-0.0720 - 0.0067i 


M[6, 1,8/3] 


13542.30 


C 


0.6375 + 1.3223i 


iV[6, 1,4/3] 


12916.72 


c 


3.7447 - 0.5884i 


0[1,3,2] 


23941.32 


m/W^GeV 


0.04 


P[3,3,2/3] 


1079.14,16853.68 


nio/lQ^^GeV 


-882.138e-*^'^^(^) 


g[8,3,0] 


0.588 


7 


0.99 


i?[8,l,0] 


0.38,1.25 


7 


-1.3849 


5[1,3,0] 


1.3766 


X 


0.8562 + 1.0667i 


t[3, 1,2/3] 


0.71,234.75,767.39,1155.96 


Ax 


2.82 




4209.99,4493.21,70106.74 


Ag 


-20.827 


f/[3,3,4/3] 


1.100 


Aa3{Mz) 


-0.011 


\/[l,2,3] 


0.979 


{M'^yW^GeV} 


0.00,12.52,1043.05 


Vr[6,3,2/3] 


9640.83 


{M^j/10-^^eV} 


0.0015,5.50,458.44 


X[3,2,5/3] 


0.359,511.842,511.842 


M^{meV) 


2.40,7.72,42.52 


F[6,2,l/3] 


0.43 


{Evals[f]}/10-^ 


0.01,54.51,4543.19 


Z[8,l,2] 


1.24 


Soft parameters 


mi = -397.337 

2 


mo = 2547.245 


Ao = -1.8262 X 10^ 


at Mx 


/i = 1.4564 X 10^ 


B = -1.4999 X 10^° 


tan/3 = 50.0000 




M| = -1.4614 X 10^° 


M2j = -1.4329 X 10^° 


Rb^ = 0.8775 


Max{\LABCD\, \Rabcd ) 


7.8076 X 10-22 GeV-i 




S/l 



Table 53: III-2-a : Column 1 contains values of the NMSGUT-SUGRY-NUHM parameters at 
Mx derived from an accurate fit to all 18 fermion data and compatible with KG constraints. See 
caption to Table 2 for explanation. 
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Parameter 


Target = Oi 


Uncert. = 5i 


Achieved = Oi 


Pull = (0, - 0,)/6, 


?/./io-6 


2.114293 


0.807660 


2.114551 


0.000320 


l/c/10-' 


1.030627 


0.170053 


1.030427 


-0.001178 


Vt 


0.388741 


0.015550 


0.388744 


0.000205 


l/d/10-^ 


5.990999 


3.492753 


6.024454 


0.009578 


l/s/10-3 


1.132757 


0.534661 


1.128624 


-0.007730 


Vb 


0.546568 


0.283669 


0.545182 


-0.004887 


1/e/lO-^ 


1.282003 


0.192301 


1.282314 


0.001614 


1/m/10-^ 


2.593678 


0.389052 


2.594359 


0.001750 


yr 


0.567317 


0.107790 


0.566957 


-0.003338 


sin ^f 2 


0.2210 


0.001600 


0.2210 


-0.0004 


sin^?3/10-^ 


28.5695 


5.000000 


28.5621 


-0.0015 


sin^lg/lO-^ 


33.6158 


1.300000 


33.6166 


0.0006 


5« 


60.0207 


14.000000 


59.9838 


-0.0026 


{m\,)/lQ-\eVf 


5.3891 


0.571243 


5.3891 


-0.0001 


{ml,)/lQ-\eVf 


1.7482 


0.349632 


1.7482 


0.0001 


sin^ ^^2 


0.2880 


0.057602 


0.2880 


-0.0003 


sin^ ^^3 


0.4594 


0.137826 


0.4595 


0.0009 


^f3(rfe^rees) 


3.7 


3.7 


4.86 




Eigenvalues{Ay) 


0.117591 


0.117741 


0.123579 




Eigenvalues{Ag) 


0.117423 


0.117569 


0.123393 




Eigenvalues{Ap) 


0.146546 


0.146698 


0.152632 




Eigenvalues{Ae) 


0.147096 


0.147254 


0.153189 




Eigenvalues{AQ) 


0.104428 


0.104576 


0.110528 




Eigenvalues{AL) 


0.133889 


0.134036 


0.139953 




An, Ah 


14.257424 


11.302210 






tti 


0.7967 + 0.0000i 


tti 


0.8979 + O.OOOOi 




^2 


0.0175 -0.0129i 


(52 


0.0271 - 0.0023i 




as 


-0.0376 - 0.0062i 


as 


-0.0325 + 0.0015i 




q;4 


0.0013 + 0.5966i 


(54 


0.0147 -0.4306i 




as 


0.0485 - 0.0181i 


"5 


0.0250 - 0.0387i 




"6 


-0.0286 - 0.0628i 


«6 


-0.0382 - 0.0523i 





Table 54: III-2-b: Fit with xx = \/Y.i=iiOi - Oi)^/5f = 0.0143. Target values, at Mx of 
the fermion yukawa couplings and mixing parameters, together with the estimated uncertainties, 
achieved values and pulls. See caption to Table 3 for explanation. 
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Parameter 


SM{Mz) 


m'^^'^iMz) 


^MSSM _ (^ ^ Amf^'^iMz) 


md/lO"^ 


2.90000 


0.52468 


2.91331 


m,/10^3 


55.00000 


9.82943 


54.75545 


rrib 


2.90000 


3.24471 


2.89957 


me/10-3 


0.48657 


0.47433 


0.48705 


rrif. 


0.10272 


0.09592 


0.10272 


m-r 


1.74624 


1.74370 


1.74662 


m„/10-3 


1.27000 


1.10824 


1.26917 


rric 


0.61900 


0.54005 


0.61848 


rat 


172.50000 


148.32422 


172.51404 



Table 55: III-2-c: Values of standard model fermion masses in GeV at Mz compared with the 
masses obtained from values of GUT derived yukawa couplings run down from M^ to Mz both 

before and after threshold corrections. Fit with xz = \/J2i=i{''^i 
0.0066. 



,MSSM 



m 



SM\2 



/{m. 



MSSM\2 



Parameter 


Value 


Parameter 


Value 


Ml 


98.92 


M-^, 


3553.79 


M2 


341.43 


M^, 


3550.53 


M3 


232.18 


M-n, 


20687.29 


^H 


1868.65 


.0(0 
^11 


-110036.65 


^h 


204.32 


40(0 


-109904.91 


^3 


11397.10 


.0(0 
^33 


-66704.24 


Ml, 


6369.22 


^11 


-132899.38 


Ml, 


6233.07 


^22 


-132898.56 


Ml, 


10527.00 


^33 


-66235.28 


^. 


712.29 


.0(d) 
^11 


-109811.44 


M-,, 


705.83 


^22 


-109810.71 


^2, 


41913.64 


40(d) 
^33 


-37711.03 


Mq, 


5387.22 


tan/3 


50.00 


Mq, 


5385.73 


^Ji{Mz) 


111532.64 


Mq, 


33269.50 


B{Mz) 


2.1254 X 10^ 


Ml 


-1.1842 X 10^° 


Ml 


-1.3659 X 10^° 



Table 56: III-2-d: Values (GeV) in of the soft Susy parameters at Mz (evolved from the soft 
SUGRY-NUHM parameters at Mx)- See caption to Table 5 for explanation. 
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Field 


MassiGeV) 


Ma 


232.18 


M,± 


341.42,111532.70 


M^o 


98.92, 341.42, 111532.68, 111532.68 


M, 


6368.872, 6232.714, 10526.788 


Mg 


1869.20, 6369.40, 190.06, 6233.87, 10327.33, 11578.52 


Mu 


3553.60, 5386.93, 3550.29, 5385.48, 20684.11, 33272.06 


M;i, 


712.77, 5387.56, 706.24, 5386.09, 33261.88, 41919.73 


Ma 


326055.13 


Mh± 


326055.14 


Mho 


326055.11 


M,o 


112.76 



Table 57: III-2-e: Spectra of supersymmetric partners calculated ignoring generation mixing 
effects. See caption to Table 6 for explanation. 



Field 


Mass{GeV) 


Ma 


232.36 


M,± 


341.50,111528.02 


M,a 


98.94, 341.50, 111528.00, 111528.00 


M, 


6232.83,6368.97,10529.123 


Me 


342.67, 1868.73, 6223.30, 6369.45, 10391.14, 11540.07 


Mu 


3550.99, 3554.30, 5367.30, 5387.36, 20677.92, 33271.88 


Ma 


709.00, 715.54, 5367.93, 5387.99, 33261.71, 41917.16 


Ma 


326179.32 


Mh± 


326179.33 


Mho 


326179.31 


Mho 


112.78 



Table 58: III-2-f: Spectra of supersymmetric partners calculated including generation mixing 
effects. See caption to Table 7 for explanation. 
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Parameter 


Value 


Field 


Masses 






[5f/(3),^f/(2),r] 


{Units ofW^GeV) 


Xx 


0.0031 


A 1,1,4] 


689.26 


Xz 


0.0012 


5[6,2,5/3] 


0.2457 


/ill/10-6 


0.0205 


C[8,2,l] 


38.27,330.90,346.79 


/i22/10-^ 


0.0364 


D[3,2,7/3] 


36.96,399.26,408.77 


/i33 


0.0027 


E[3,2,l/3] 


0.40,27.84,27.84 


/11/IO-6 


0.0774 -0.1295i 




30.102,433.12,477.24 


/12/IO-6 


-2.0249 - 0.0687i 


F[l,l,2] 


6.43,6.43 


/13/IO-5 


0.0636 + 0.0485i 




26.27,346.91 


/22/IO-5 


6.5620 -4.8133i 


^[1,1,0] 


0.063,0.48,0.60 


/23/10-^ 


2.0322 + 2.3187i 




0.600,26.21,26.37 


/33/10-=^ 


-1.0357 + 0.4496i 


/i[l,2,l] 


1.019,21.54,35.95 


^12/10-^ 


0.0555 + 0.1219i 




600.38,620.06 


^13/10-^ 


-0.0665 + 1.8279i 


/[3, 1,10/3] 


0.85 


^23/10-^ 


6.3795 + 5.6374i 


J[3, 1,4/3] 


0.959, 14.94, 14.94 


A/10-2 


-0.4433 -1.0332i 




46.51,403.62 


V 


-10.5295 + 2.6182i 


K[3, 1,8/3] 


53.75,492.39 


P 


0.6481 - 2.2964i 


L[6, 1,2/3] 


25.26,794.20 


k 


0.0149 - 0.0810i 


M[6, 1,8/3] 


802.14 


C 


1.6389 + 0.5058i 


Ar[6, 1,4/3] 


802.65 


c 


1.0205 + 0.5627i 


0[1,3,2] 


1583.99 


m/W^GeV 


0.03 


P[3,3,2/3] 


15.91,1232.28 


nio/lQ^^GeV 


-21.872e-^^''»W 


g[8,3,0] 


0.711 


7 


3.97 


i?[8,l,0] 


0.27,1.07 


7 


-2.9091 


5[1,3,0] 


1.2049 


X 


0.9564 + 0.6718i 


t[3, 1,2/3] 


0.78,20.40,50.58,87.07 


Ax 


1.21 




275.21,368.92,7729.25 


Ag 


1.734 


f/[3,3,4/3] 


1.018 


Aa3{Mz) 


-0.016 


\/[l,2,3] 


0.714 


{M'^yW^GeV} 


0.01,13.27,525.55 


Vr[6,3,2/3] 


948.60 


{M^j/10-^^eV} 


0.6555,1070.44,42380.30 


X[3,2,5/3] 


0.241,29.690,29.690 


M^{meV) 


2.42,7.81,43.04 


F[6,2, 1/3] 


0.30 


{Evals[f]}/10-^ 


0.19,306.30,12126.94 


Z[8,l,2] 


1.06 


Soft parameters 


mi = -143.002 

2 


mo = 3347.253 


Ao = -1.0946 X 10^ 


at Mx 


/i = 8.8928 X 10^ 


B = -5.4965 X 10^ 


tan/3 = 51.5000 




M| = -5.8300 X 10^ 


M2j = -5.5566 X 10^ 


Rb. = 2.6676 


Max{\LABCD\, \Rabcd ) 


5.5297 X 10-23 GeV^i 




S/J. 



Table 59: III-3-a : Column 1 contains values of the NMSGUT-SUGRY-NUHM parameters at 
Mx derived from an accurate fit to all 18 fermion data and compatible with RG constraints. See 
caption to Table 2 for explanation. 
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Parameter 


Target = Oi 


Uncert. = 6i 


Achieved = Oi 


Pull = (0, - 0,)/6, 


?/./io-6 


2.121322 


0.810345 


2.121299 


-0.000028 


l/c/10-' 


1.034062 


0.170620 


1.033951 


-0.000648 


Vt 


0.392382 


0.015695 


0.392383 


0.000076 


l/d/10-^ 


7.238271 


4.219912 


7.243490 


0.001237 


l/s/10-3 


1.372396 


0.647771 


1.371099 


-0.002002 


Vb 


0.512784 


0.266135 


0.513218 


0.001628 


1/e/lO-^ 


1.291572 


0.193736 


1.291652 


0.000409 


1/m/10-^ 


2.645649 


0.396847 


2.645665 


0.000038 


yr 


0.580184 


0.110235 


0.580135 


-0.000444 


sin ^f 2 


0.2210 


0.001600 


0.2210 


-0.0001 


sin^?3/10-^ 


28.7475 


5.000000 


28.7473 


0.0000 


sin^lg/lO-^ 


33.8250 


1.300000 


33.8252 


0.0002 


5« 


60.0207 


14.000000 


60.0163 


-0.0003 


{m\,)/lQ-\eVf 


5.5085 


0.583902 


5.5085 


0.0000 


{ml,)/lQ-\eVf 


1.7916 


0.358315 


1.7916 


0.0000 


sin^ ^^2 


0.2874 


0.057476 


0.2874 


0.0002 


sin^ ^^3 


0.4571 


0.137144 


0.4571 


-0.0004 


^f3(rfe^rees) 


3.7 


3.7 


5.41 




Eigenvalues{Ay) 


0.050027 


0.050031 


0.050037 




Eigenvalues{Ag) 


0.047655 


0.047660 


0.047665 




Eigenvalues{Ap) 


0.053435 


0.053439 


0.053445 




Eigenvalues{Ae) 


0.060550 


0.060554 


0.060560 




Eigenvalues{AQ) 


0.047387 


0.047391 


0.047396 




Eigenvalues{AL) 


0.055538 


0.055543 


0.055547 




An, Ah 


78.220842 


64.106181 






tti 


0.7970 + O.OOOOi 


ai 


0.8786 + O.OOOOi 




^2 


0.0641 + 0.0154i 


(52 


0.0540 + 0.0624i 




as 


-0.0400 - 0.0369i 


as 


-0.0557 -0.0165i 




q;4 


-0.4642 + 0.1746i 


(54 


0.3473 - 0.002H 




as 


0.1479 + 0.0832i 


"5 


0.0692 + 0.0216i 




"6 


0.1139 -0.2641i 


"6 


0.1278 -0.2750i 





Table 60: III-3-b: Fit with xx = \/Y.i=iiOi - Oi)^/5f = 0.0031. Target values, at Mx of 
the fermion yukawa couplings and mixing parameters, together with the estimated uncertainties, 
achieved values and pulls. See caption to Table 3 for explanation. 
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Parameter 


SM{Mz) 


m'^^'^iMz) 


^MSSM _ (^ ^ Amf^'^iMz) 


md/lO"^ 


2.90000 


0.62377 


2.90229 


m,/10^3 


55.00000 


11.80718 


54.95706 


rrib 


2.90000 


3.08261 


2.90002 


me/10-3 


0.48657 


0.47238 


0.48647 


rrif. 


0.10272 


0.09671 


0.10272 


m-r 


1.74624 


1.74579 


1.74575 


m„/10-3 


1.27000 


1.10616 


1.26997 


rric 


0.61900 


0.53916 


0.61901 


rat 


172.50000 


149.41400 


172.44284 



Table 61: III-3-c: Values of standard model fermion masses in GeV at Mz compared with the 
masses obtained from values of GUT derived yukawa couplings run down from M^ to Mz both 

before and after threshold corrections. Fit with xz = \/J2i=i{''^i 
0.0012. 



,MSSM 



m 



SM\2 



/{m. 



MSSM\2 



Parameter 


Value 


Parameter 


Value 


Ml 


98.22 


M^, 


4257.81 


M2 


273.31 


M-n, 


4256.95 


M3 


386.79 


Mn, 


15648.62 


^\ 


1055.61 


Ml) 
^11 


-66907.20 


^% 


149.79 


40(0 

^22 


-66822.57 


^3 


12367.95 


40(0 


-39462.25 


Mu 


5350.98 


40W 
^11 


-79457.41 


Mu 


5300.00 


4OH 
^22 


-79456.87 


Ml, 


10339.21 


^33 


-39415.43 


% 


2146.46 


.0(d) 

^11 


-67055.77 


M-- 

d2 


2145.58 


^22 


-67055.16 


^Ja 


25986.72 


40(d) 
^33 


-25007.87 


Mq, 


4117.71 


tan/3 


51.50 


Mq, 


4117.04 


KMz) 


69115.80 


Mq, 


21502.97 


B{Mz) 


7.5656 X 10^ 


M% 


-4.7132 X 10^ 


Ml 


-5.2004 X 10^ 



Table 62: III-3-d: Values (GeV) in of the soft Susy parameters at Mz (evolved from the soft 
SUGRY-NUHM parameters at Mx)- See caption to Table 5 for explanation. 



Field 


MassiGeV) 


Ma 


386.79 


M,± 


273.31,69115.90 


M^o 


98.22, 273.31, 69115.86, 69115.87 


M, 


5350.564, 5299.581, 10338.994 


Mg 


1056.57, 5351.20, 141.75, 5300.63, 10298.54, 12402.01 


Mu 


4117.34, 4257.65, 4116.48, 4256.98, 15643.84, 21507.38 


M;i, 


2146.61, 4118.16, 2145.70, 4117.51, 21489.76, 25997.73 


Ma 


197426.75 


Mh± 


197426.76 


Mho 


197426.73 


M,o 


117.60 



Table 63: III-3-e: Spectra of supersymmetric partners calculated ignoring generation mixing 
effects. See caption to Table 6 for explanation. 



Field 


Mass{GeV) 


Ma 


386.90 


M,± 


273.35,69113.11 


M,o 


98.23, 273.35, 69113.07, 69113.08 


M, 


5299.63,5350.60,10339.851 


Me 


118.46, 1056.41, 5300.46, 5351.24, 10300.80, 12402.58 


Mu 


4113.23, 4117.53, 4257.15, 4257.85, 15641.39, 21506.49 


Ma 


2145.97, 2146.92, 4114.23, 4118.35, 21488.89, 25996.33 


Ma 


197503.48 


Mh± 


197503.50 


Mho 


197503.46 


Mho 


117.80 



Table 64: III-3-f : Spectra of supersymmetric partners calculated including generation mixing 
effects. See caption to Table 7 for explanation. 
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Parameter 


Value 


Field 


Masses 






[SU{'i),SU{2),Y] 


{Units ofW^GeV) 


Xx 


0.0469 


Al,l,4 


153.03 


Xz 


0.0139 


5[6,2,5/3] 


0.3445 


/ill/10-6 


-0.0162 


C[8,2, 1] 


5.06,191.08,192.38 


/i22/10-^ 


-0.0156 


Dp, 2, 7/3] 


4.50,104.56,129.69 


/i33 


0.0033 


^[3,2,1/3] 


0.41,3.63,5.83 


/11/IO-6 


0.0145 + 0.1575i 




5.828,154.26,178.97 


/12/IO-6 


-1.2719 + 0.4383i 


F[l,l,2] 


1.06,1.06 


/13/IO-5 


0.3213 + 0.1092i 




2.94,77.81 


/22/IO-5 


5.9014 - 6.6842i 


G[1,1,0] 


0.097,0.28,0.28 


/23/10-^ 


2.0445 + 3. 1164i 




0.783,16.65,17.07 


/33/10-=^ 


-1.0245 + 0.7289i 


/i[l,2, 1] 


1.382,3.46,4.74 


^12/10-^ 


0.0428 - 0.095H 




134.98, 176.94 


^13/10-^ 


-1.5145 - 2.2227i 


/[3, 1,10/3] 


1.41 


^23/10-^ 


2.9856 + 6. 1462i 


J[3, 1,4/3] 


1.370,2.90,2.90 


A/10-2 


-7.7575 - 1.6636i 




7.29,119.06 


V 


-14.6064 + 1.2176i 


i^[3, 1,8/3] 


7.99, 140.68 


P 


0.4221 - 1.5473i 


L[6, 1,2/3] 


4.16,230.75 


k 


0.0267 -0.073H 


M[6, 1,8/3] 


236.13 


C 


1.0476 + 0.4885i 


Ar[6, 1,4/3] 


228.10 


c 


1.5631 -0.4286i 


0[1,3,2] 


421.09 


m/W^GeV 


0.03 


P[3,3,2/3] 


1.17,372.51 


nio/lQ^^GeV 


-3.437e-*^'^3(^) 


g[8,3,0] 


1.242 


7 


4.58 


i?[8, 1, 0] 


0.37,1.45 


7 


-2.2542 


5[1,3,0] 


1.6823 


X 


0.8706 + 0.4654i 


t[3, 1,2/3] 


1.02,2.78,7.81,9.33 


Ax 


0.43 




152.15,160.68,3371.92 


Ag 


20.355 


f/[3,3,4/3] 


1.485 


Aa3{Mz) 


-0.009 


\/[l,2,3] 


0.944 


{M'^yW^GeV} 


0.01,7.01,280.82 


1^(6,3,2/3] 


356.14 


{M^j/10-^^eV} 


6.3962,5020.83,201135.16 


X[3,2,5/3] 


0.316,6.340,6.340 


M^{meV) 


2.32,7.54,41.49 


F[6,2,l/3] 


0.39 


{Evals[f]}/10-^ 


0.43,340.95,13658.39 


Z[8,l,2] 


1.45 


Soft parameters 


mi = -258.670 

2 


mo = 4903.144 


Ao = -1.4665 X 10^ 


at Mx 


/i = 1.2975 X 10^ 


B = -9.5887 X 10^ 


tan/3 = 51.0000 




M| = -1.4950 X 10^° 


M2j = -1.4315 X 10^° 


Rbr = 3.7388 


Max{\LABCD\, \Rabcd ) 


4.7911 X 10-22 GeV-i 




SM 



Table 65: III-4-a : Column 1 contains values of the NMSGUT-SUGRY-NUHM parameters at 
Mx derived from an accurate fit to all 18 fermion data and compatible with RG constraints. See 
caption to Table 2 for explanation. 
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Parameter 


Target = Oi 


Uncert. = 6i 


Achieved = Oi 


Pull = (0, - 0^)/6^ 


1/./10-6 


2.106535 


0.804696 


2.106685 


0.000186 


l/c/10-3 


1.026841 


0.169429 


1.027073 


0.001369 


yt 


0.368035 


0.014721 


0.368036 


0.000048 


Vd/io-' 


6.294474 


3.669678 


6.323679 


0.007958 


ys/io~' 


1.194492 


0.563800 


1.201918 


0.013170 


Vb 


0.431150 


0.223767 


0.439556 


0.037566 


l/e/lO-4 


1.180231 


0.177035 


1.179590 


-0.003620 


y,/io-^ 


2.495809 


0.374371 


2.496108 


0.000799 


Vr 


0.530719 


0.100837 


0.528387 


-0.023127 


sin ^f 2 


0.2210 


0.001600 


0.2210 


0.0001 


sin ^^3/ 10-4 


29.3942 


5.000000 


29.3933 


-0.0002 


sin ^13/10-3 


34.5852 


1.300000 


34.5835 


-0.0013 


5« 


60.0209 


14.000000 


60.0050 


-0.0011 


iml,)/WieVr 


5.1427 


0.545127 


5.1428 


0.0002 


{ml,)/W-'{eVr 


1.6651 


0.333025 


1.6650 


-0.0003 


sin^ ^^2 


0.2885 


0.057694 


0.2885 


-0.0003 


sin^ 023 


0.4611 


0.138326 


0.4610 


-0.0005 


Ofsidegrees) 


3.7 


3.7 


4.54 




Eigenvalues{Au) 


0.010307 


0.010309 


0.010310 




Eigenvalues{Ai) 


0.006885 


0.006887 


0.006888 




Eigenvalues{Ap) 


0.000492 


0.000494 


0.000494 




Eigenvalues{Ag) 


0.010759 


0.010761 


0.010762 




Eigenvalues{AQ) 


0.013199 


0.013201 


0.013201 




Eigenvalues{AL) 


0.010228 


0.010230 


0.010230 




Ah, Ah 


58.562855 


50.761700 






ai 


0.7564 + O.OOOOi 


ai 


0.8100 + 0.0000i 




a2 


0.0475 + 0.0440i 


(52 


-0.0078 + 0.0789i 




as 


0.0250 - 0.0486i 


"3 


-0.0002 - 0.0292i 




a4 


-0.5034 - 0.0920i 


(54 


0.2831 + 0.1536i 




as 


0.2065 + 0.0095i 


«5 


0.1215 - 0.0343i 




«6 


0.0602 - 0.3353i 


"6 


0.0906 - 0.4570i 





Table 66: III-4-b: Fit with xx = ^/ElliiOi 
explanation. 



Oi)'^/6f = 0.0469. See caption to Table 3 for 
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Parameter 


SM{Mz) 


m'^^'^iMz) 


^MSSM _ (^ ^ Amf^'^iMz) 


md/lO"^ 


2.90000 


0.58241 


2.89785 


m,/10^3 


55.00000 


11.06967 


55.04827 


rrib 


2.90000 


2.97388 


2.91949 


me/10-3 


0.48657 


0.46121 


0.48338 


rrif. 


0.10272 


0.09755 


0.10210 


m-r 


1.74624 


1.73358 


1.73257 


m„/10-3 


1.27000 


1.11823 


1.26917 


rric 


0.61900 


0.54517 


0.61876 


rat 


172.50000 


147.05975 


172.17090 



Table 67: III-4-c: Values of standard model fermion masses in GeV at Mz compared with the 
masses obtained from values of GUT derived yukawa couplings run down from M^ to Mz both 

before and after threshold corrections. Fit with xz = \/J2i=i{''^i 
0.0139. 



,MSSM 



m 



SM\2 



/{m. 



MSSM\2 



Parameter 


Value 


Parameter 


Value 


Ml 


94.02 


M-^, 


5778.49 


M2 


284.49 


M^, 


5778.19 


M3 


269.37 


M-n, 


32851.51 


^H 


100.27 


.0(0 
^11 


-94185.07 


^h 


122.47 


40(0 

^22 


-94074.32 


^3 


28577.12 


.0(0 
^33 


-59361.23 


Ml, 


8090.35 


^11 


-107534.77 


Ml, 


8090.50 


^22 


-107534.06 


Ml, 


21787.08 


^33 


-56550.41 


^. 


1167.18 


.0(d) 
^11 


-94095.80 


M-,, 


1168.84 


^22 


-94095.08 


^2, 


40997.04 


40(d) 
^33 


-42029.00 


Mq, 


5546.60 


tan/3 


51.00 


Mq, 


5546.62 


^Ji{Mz) 


108204.37 


Mq, 


37188.99 


B{Mz) 


1.9013 X 10^ 


Ml 


-1.1922 X 10^° 


Ml 


-1.2689 X 10^° 



Table 68: III-4-d: Values (GeV) in of the soft Susy parameters at Mz (evolved from the soft 
SUGRY-NUHM parameters at Mx)- See caption to Table 5 for explanation. 
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Field 


MassiGeV) 


Ma 


269.37 


M,± 


284.48, 108204.44 


M^o 


94.02, 284.48, 108204.41, 108204.42 


M, 


8090.077, 8090.228, 21786.977 


Mg 


110.00, 8090.50, 111.64, 8090.93, 21780.86, 28581.94 


Mu 


5546.33, 5778.37, 5546.23, 5778.19, 32848.10, 37192.53 


M;i, 


1167.48, 5546.94, 1169.07, 5546.97, 37176.73, 41008.21 


Ma 


311453.73 


Mh± 


311453.74 


Mho 


311453.72 


M,o 


128.46 



Table 69: III-4-e: Spectra of supersymmetric partners calculated ignoring generation mixing 
effects. See caption to Table 6 for explanation. 



Field 


Mass{GeV) 


Ma 


269.50 


M,± 


284.54, 108200.79 


M,o 


94.04, 284.54, 108200.77, 108200.77 


M, 


8090.11,8090.27,21787.520 


Me 


88.50, 108.38, 8090.53, 8090.99, 21781.65, 28582.63 


Mu 


5546.48, 5569.26, 5778.35, 5778.53, 32847.38, 37188.53 


M,: 


1168.04, 1169.61, 5547.10, 5569.98, 37172.75, 41007.52 


Ma 


311542.91 


Mh± 


311542.92 


Mho 


311542.91 


Mho 


128.51 



Table 70: Spectra of supersymmetric partners calculated including generation mixing effects. See 
caption to Table 7 for explanation. 



83 



References 

[1] C.S. Aulakh and R.N. Mohapatra, CCNY-HEP-82-4 April 1982, CCNY-HEP-82-4- 
REV, Jun 1982 , Phys. Rev. D28, 217 (1983). 



[2] 
[3] 

[4; 

[5 
[6 
[7 



T.E. Clark, T.K.Kuo, and N.Nakagawa, Phys. lett. 115B, 26(1982). 

M. S. Carena, M. Olechowski, S. Pokorski and C. E. M. Wagner, Nucl. Phys. B 426, 
269 (1994) [arXiv:hep-ph/9402253 1 . B. Ananthanarayan, Q. Shafi and X. M. Wang, 
|arXiv:hep-ph/9311225 |; R. Rattazzi, U. Sarid and L. J. Hall, |arXiv:hep-ph/9405313( 

L. J. Hall, R. Rat tazzi and U. Sarid, Phys. Rev. D 50 (1994) 7048 
|arXiv:hep-ph/ 9306309]; R. Hempfling, Phys. Rev. D 49, 6168 (1994). 

C. S. Aulakh, B. Bajc, A. Melfo, A. Rasin and G. Senjanovic, Nucl. Phys. B 597, 89 
(2001) || arXiv:hep-ph/ 0004031|. 

C S. Aulakh, B. Bajc, A. Melfo, G. Senjanovic and F. Vissani, Phys. Lett. B 588, 
196 (2004) |arXiv:hep-ph/0306242| . 

C. S. Aulakh and S. K. Garg, Nucl. Phys. B 757, 47 (2006) |arXiv:hep-ph/05 12224] . 

L. J. Hall, R. Rattazzi and U. Sarid, Phys. Rev. D 50, 7048 (1994) 
|arXiv:hep-ph/9306309 1 ; G. Anderson, S. Raby, S. Dimopoulos, L. J. Hall and 
G. D. Starkman, Phys. Rev. D 49, 3660 (1994) |arXiv:hep-ph/9308333|; L. J. Hall and 
S. Raby, Phys. Rev. D 51, 6524 (1995) larXiv:hep-ph/9501298|; K.S.Babu and S.M. 



Barr, Phys. Rev. D51, 2463 (1995), hep-ph/9409285i R. Rattazzi and U. Sarid, Phys. 



Rev. D 53, 1553 (1996) |'arXiv:hep-ph/9505428|; S.M. Barr and S.Raby, Phys. Rev. 
Lett. 79, 4748 (1997), hep-ph/9705366i C. H. Albright, K. S. Babu and S. M. Barr, 
Phys. Rev. Lett. 81, 1167 (1998) |arXiv:hep-ph/98023l4|; H. Baer, M. A. Diaz, J. Fer- 
randis, X. Tata, Phys.Rev.D61:111701,2000, [hep /ph 9907211] ; K. S. B abu, J. C. Pati 
and F. Wilczek, Nucl. Phys. B 566, 33 (2000) |ar Xiv:hep-ph/9812538j ; C. H. Albright 
and S. M. Barr, Phys. Lett. B 461, 218 (1999) [arXiv:hep-ph/9906297]; Phys. Rev. 
Lett. 85, 244 (2000) |arXiv:hep-ph/0002155]; T. Blazek, R. Dermisek and S. Raby, 
Phys. Rev. Lett. 88, 111804 (2002) |arXiv:hep-ph/0107097|; T. Blazek, R. Dermisek 
and S. Raby, Phys. Rev. D 65, 115004 (2002) [arXiv:hep-ph/0201081]. 



[9] C.S.Aulakh and A. Girdhar, ■hep-ph/0204"097l v2 August 2003; v4, 9 February, 2004 



Int. J. Mod. Phys. A 20, 865 (2005) 

[10] B. Baj c, A. Melfo, G. Senjanovic and F. Vissani, Phys. Rev. D 70, 035007 (2004) 
[arXiv:hep-ph/0402122| . 

[11] T. Fukuyama, A. Ilakovac, T. Kikuchi, S. Meljanac and N. Okada, Eur. Phys. J. C 
42, 191 (2005) arXiv:hep-ph/0401213 vl.,v2. 

[12] C. S. Aulakh and A. Girdhar, Nucl. Phys. B 711 (2005) 275 |arXiv:hep-ph/0405074l . 



[13] C. S. Aulakh, From germ to bloom, arXiv:hep-ph/0506291 



84 



[14] C S. Aulakh and S. K. Garg, The new minimal supersymmetric GUT, 
|arXiv:hep-ph/0612021vl; The new minimal supersymmetric GUT: Spectra, RG anal- 
ysis and fitting formulae, jarXiv:hep-ph/061202r| V2. 

[15] B. Bajc, A. Melfo, G. Senjanovic and F. Vissani, Phys.Lett.B634:272-277,2006 

[16] S. Bertolini, T. Schwetz and M. Malinsky, Phys. Rev. D 73 (2006) 115012 
|arXiv:hep-ph7060 5006| . 

[17] C. S. Aulakh and S. K. Garg, "NMSGUT II: Pinning the Nmsgut@LHC," 
larXiv:0807.09T7V l [hep-ph]; "NMSGUT II: Fermion Fits and soft Spectra" 
larXiv:0807.09T7V 2 [hep-ph] . 

[18] P. Langacker and N. Polonsky, Phys. Rev. D 47, 4028 (1993) 
|arXiv:hep-ph/921023 5 1 ;P. Langacker and N. Polonsky, Phys. Rev. D 52, 3081 
(1995) |arXiv:hep-ph795032l4|. 



[19] C. S. Aulakh and S. K. Garg, Mod. Phys. Lett. A 24 (2009) 1711 |arXiv:0710.40T8l 
[hep-ph]]. 

[20] W. J. Marciano and G. Senjanovic, Phys. Rev. D 25 (1982) 3092. 

[21] U. Amaldi, W. de Boer and H. Furstenau, Phys. Lett. B 260, 447 (1991). 

[22] K. Tobe and J. D. Wells, Nucl. Phys. B 663 (2003) 123 |arXiv:hep-ph/0301015| . 



[23] M. S. Carena, M. Olechowski, S. Pokorski and C. E. M. Wagner, Nucl. Phys. B 
426, 269 (1994) |a rXiv:hep-ph/9 402253|; M. S. Carena, D. Garcia, U . Nierste and 
C. E. M. Wagner, Nucl. Phys. B 577, 88 (2000) | arXiv:hep-ph/9912516 ); A. J. Buras, 
P. H. Chankowski, J. Rosiek and L. Slawianowska, Nucl. Phys. B 659, 3 (2003) 
farXiv:hep-ph/0210l45] ;A. Freitas, E. Gasser and U. Haisch, Phys. Rev. D 76 (2007) 
014016 [ar Xiv:hep-ph/0702267] . 

[24] G. Anderson, S. Raby, S. Dimopoulos, L. J. Hall and G. D. Starkman, Phys. Rev. D 
49 (1994) 3660 [arXiv:hep-ph/9308333| . 



[25] V. Lucas and S. Raby, Phys. Rev. D 55 (1997) 6986 arXiv:hep-ph/9610293]. 



[26] T. Blazek , R. Dermisek and S. Raby, Phys. Rev. D 65 (2002) 115004 
[arXiv:hep-ph/0201081]; T. Blazek, R . Dermisek and S. Raby, Phys. Rev. Lett. 88 
(2002) 111804 |arXiv:hep-ph/0107097l . 



[27] H. Baer and J. Ferrandis, Phys. Rev. Lett. 87 (2001) 211803 |arXiv:hep-ph/0106352| . 

[28] D. G. Lee, Phys. Rev. D 49 (1994) 1417. 

[29] K.S.Babu and R.N.Mohapatra, Phys. Rev. Lett. 70(1993)2845. 

[30] C. S. Aulakh, Phys. Lett. B 661, 196 (2008) [arXiv:0710.3"945l [hep-ph]]. 



85 



[31] K. Y. Oda, E. Takasugi, M. Tanaka and M. Yoshimura, Phys. Rev. D 59, 
055001 (1999) |arXiv:hep-ph/9808241|; K. Matsuda, Y. Koide, T. Fukuyama and 
H. Nishiura, Phys. Rev. D 65, 033008 (2002) [Erratum-ibid. D 65, 079904 (2002)] 
jarXiv:h ep-ph/ 0108202j ; K. Matsuda, Y. Koide and T. Fukuyama, Phys. Rev. 
D 64, 053015 (2001) |arXiv:hep-ph/001002"6) . N. Oshimo, ; Phys. Rev. D 66, 
095010 (2002) |arXiv:hep-ph/0206239|; N. Oshimo, NucL Phys. B 668, 258 (2003) 



tarXiv:hep-ph/0305166 ];K. S. Babu and C. Macesanu, [arXiv:hep-ph/0505200^ 



[32] B. Bajc, G. Senjanovic and F. Vissani, Phys. Rev. Lett. 90 (2003) 051802 
|arXiv:hep-ph/0210207]; H. S. Goh, R. N. Mohapatra and S. P. Ng, Phys. Lett. B 
570, 215 (2003) [arXiv:hep-ph/0303055|. H. S. Goh, R. N. Mohapatra and S. P. Ng, 
Phys. Rev. D 68, 115008 (2003) |arXiv:hep-ph/030819 7]. H. S . Goh, R. N. Mohap- 
atra and S. Nasri, Phys. Rev. D 70 (2004) 075022 [ arXiv:hep-ph/0408139]; B. Bajc , 
G. Senjanovic and F. Vissani, Phys. Rev. D 70, 093002 (2004) |arXiv:hep-ph/0402140| ; 
B. Bajc, G. Senjanovic and F. Vissani, arXiv:hep-ph/0110310;K. S. Babu and C. Mace- 
sanu, Phys. Rev. D 72, 115003 (2005) [arXiv:hep-ph/0505200j . 

[33] S. Antusch and M. Spinrath, Phys. Rev. D 78, 075020 (2008) [ arXiv: 0804.07X71 [hep- 

ph]]. 

[34] G. Ross and M. Serna, Phys. Lett. B 664, 97 (2008) [arXiv: 0704. 1 2481 [hep-ph]]. 



[35 

[36 

[37 
[38; 
[39 

[40 

[41 

[42; 
[43; 
[44; 



L. Lavoura, H. Kuhbock and W. Grimus, NucL Phys. B 754 (2006) 1 
I ar Xiv:hep-ph/0603259j . 

V.V. Dixit and M. Sher, Phys. Rev. D40,3765(1989). 



B. D. Wright, jhij-ph/9404217| (1994). 



S. Weinberg, Phys. Rev. Lett. 43 (1979) 1566. L.J.Hall, Nucl. Phys. B178,75(1981). 



T2K Collaboration announcement, June 15, 2011.http://www.kek.jp/intra-e/press/2011/J-PARC_T2Kne 



MINOS Collaboration announcement, see|http://www.symmetrymagazine.org/breaking/2011/06/24/fern 
C. S. Aulakh, Fermion mass hierarchy in the Nu MS GUT. I: The real core, 



arXiv:hep-ph/0602132 



C. S. Aulakh, MSGUT Reborn ? arXiv:hep-ph/0607252 



Y.Fukuda et a/.,[Super-Kamiokande collaboration], Phys. Rev. Lett. 81,1562(1998) 

D. M. Pierce, J. A. Bagger, K . T. Matchev and R. j. Zhang, Nucl. Phys. B 491, 3 
(1997) |arXiv:hep-ph/9606211| . 



[45] K. S. Babu, C. N. L eung and J. T. Pantaleone, Phys. Lett. B 319, 191 (1993) 
|arXiv:hep-ph/9309223 1 . 



86 



[46] S. Antusch, M . Drees, J. K ersten, M. Lindner and M. Ratz, Phys. Lett. B 519, 238 
(2001) |arXiv:hep-ph/0108005|; S. Antusch, M. Drees, J. Kersten, M. Lindner and 
M. Ratz, |arXiv:hep-ph/0110366|; S. Antusch, J. Kersten, M. Lindner, M. Ratz and 
M. A. Schmidt, JHEP 0503, 024 (2005) | arXiv:hep-ph/0501272 ]. 

[47] Y.Hayato et al. ,[Super-Kaniiokande collaboration], Phys. Rev. Lett. 83,1529(1999); 
K.Kobayashi et al. , [Super-Kaniiokande collaboration] |hep-ex/0502026j . 



[48; 
[49; 
[5o; 

[51 

[52; 
[53; 

[54 

[55 
[56; 

[57 

[58; 

[59 
[60 

[61 
[62 
[63 

[64 



CLEO: S.Chen et al, Phys. Rev. Lett 87,251807(2001). 

BaBar: B. Aubert et a/., Phys. Rev. D72,052004(2005),Phys. Rev. D73,012005(2006). 

Belle: A.Limosani et al, Phys. Rev. Lett 103,241801(2009). 

K. Nakamura et al. (Particle Data Group), J. Phys. G 37, 075021 (2010) 

Numerical Recipes in Fortran 90, Second Edition(1996), Cambridge Univ. Press, 
W.H. Press, S.A.Teukolsky, W.T.Vetterling, and B.P.Flannery. 



T. Goto and T. Nihei, Phys. Rev. D 59, 115009 (1999) [arXiv:hep-ph/9808255 |. 

A. Kusenko, P. Langacker and G. Segre, Phys. Rev. D 54 (1996) 5824 
[arXiv:hep-ph/9602414 | . 

G. Dvali, ^''Safety of Minkowski Vacuum" larXiv: 1107. 09561 [hep-th]. 



W. Porod, Comput. Phys. Commun. 153 (2003) 275 |arXiv:hep-ph/0301101] ; W. 
Porod and F. Staub larXiv: 1104. 15731 



C. S. Aulakh, ^^ Truly minimal unification: Asymptotically strong panacea'^ 



?" 



arXiv:hep-ph/0207150 



C. S. Aulakh, ^''Taming asymptotic strength" arXiv:hep-ph/02 10337 

M. Fukugita and T. Yanagida, Phys. Lett. B 174 (1986) 45. 

PAMELA collab, O. Adriani et al.. Nature 458, 607 (2009). 

ATIC collab, J. Chang et al. Nature (London) 456, 362 (2008). 

FERMI/LAT collab, A. A. Abdo et al., Phys. Rev. Lett. 102, 181101 (2009). 

HESS collab, F. Aharonian et al., Phys. Rev. Lett. 101, 261104, (2008). 

S. F. King and M. Oliveira, king Phys. Rev. D 63 (2001) 015010 
[arXiv:hep-ph/0008183]. 



[65] K. Kowalska, Acta Phys. Polon. B 33 (2002) 1823 tarXiv:hep-ph/0203168| . 

S. Antusch, M. Drees, J. Kersten, M. Lindner and M. Ratz, Phys. Lett. B 519, 238 
(2001) [arXiv:hep-ph/0108005|; S. Antusch, M. Drees, J. Kersten, M. Lindner and 
M. Ratz, |arXiv:hep-ph/0110366|; S. Antusch, J. Kersten, M. Lindner, M. Ratz and 
M. A. Schmidt, JHEP 0503, 024 (2005) |arXiv:hep-ph/0501272|. 



[66] B. Baj c, P. F ileviez Perez and G. Senjanovic, Phys. Rev. D 66 (2002) 075005 
|arXiv:hep-ph/0204311j . 

[67] J. M. Frere, D. R. T. Jones and S. Raby, Nucl. Phys. B 222, 11 (1983); L. Alvarez- 
Gaume, J. Polchinski and M. B. Wise, Nucl. Phys. B 221, 495 (1983); J. P. Derendinger 
and C. A. Savoy, Nud. Phys. B 237, 307 (1984); C. Kounnas, A. B. Lahanas, D. V. 
Nanopoulos and M. Quiros, Nucl. Phys. B 236, 438 (1984). 

[68] J. A. Casas and S. Dimopoulos, Phys. Lett. B 387, 107 (1996) |arXiv:hep-ph/9606237| . 

[69] A. Riotto and E. Roulet, Phys. Lett. B 377, 60 (1996) | arXiv:hep-ph/9512401| . 

[70] G. Gamberini, G. Ridolfi, F. Zwirner, Nucl. Phys. B331(1990) 331. 

[71] J.Elhs, G.Ridolfi, F.Zwirner, Phys. Lett. B262(1991) 477. 

[72] U. Ellwanger, Phys. Lett. B141(1984) 435. 

[73] A. D. Box and X. Tata, Phys. Rev. D 77 (2008) 055007 [Erratum-ibid. D 82 
(2010) 119904] [arXiv:0712.2858 [hep-ph]]; A. D. Box and X. Tata, Phys. Rev. D 
79 (2009) 035004 [Erratum-ibid. D 82 (2010) 119905] jarXiv:0810.5765] [hep-ph]]; 

A. D. Box, "Renormalisation Group Analysis of Supersymmetric Particle Interactions, 
larXiv:0811.2444l [hep-ph]. (PhD Thesis , Univ. of Hawaii). 

[74] H. Arason, D. J. Castano, B. Keszthelyi, S. Mikaelian, E. J. Piard, P. Ramond and 

B. D. Wright, Ph ys. Rev. D 4 6 (1992) 3945; H. Arason, D. J. Castano, E. J. Piard 
and P. Ramond, |arXiv:hep-ph/9204225l . 

[75] S. P. Martin and M. T. Vaughn, Phys. Rev. D 50 (1994) 2282 |arXiv:hep-ph/9311340| ; 

[76] M. Ciuchini, A. Masiero, P. Paradisi, L. Silvestrini, S. K. Vempati and O. Vives, Nucl. 
Phys. B 783 (2007) 112 |a rXiv:hep-ph/0702144] . 

[77] V. A. Novikov, M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 
229 (1983) 381 ; M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Sov. J. Nucl. 
Phys. 43 (1986) 1028 . 

[78] J. Rosiek, P. Chankowski, A. Dedes, S. Jager and P. Tanedo, Comput. Phys. Comniun. 
181 (2010) 2180 [a rXiv: 1003.4260 [hep-ph]]. 

[79] F. Gabbiani, E. Gabrielh, A. Masiero and L. Silvestrini, Nucl. Phys. B 477 (1996) 321 
[arXiv:hep-ph/9604387j . 



